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Summary 

An  investigation  of  soot  formation  in  laminar  diffusion  flames  showed  that  soot  particle  surface 
growth  under  laminar  diffusion  flame  conditions  ceases  because  of  the  depletion  of  hydrocarbon  species 
and  not  soot  particle  reactivity  loss  due  to  thermal  aging  of  the  particles.  This  result  was  obtained  through 
direct  species  concentration  measurements  under  well-controlled  conditions,  while  the  particle  reactivity 
effects  were  calculated  based  on  premixed  flame  results  along  with  particle  temperature/time  information 
available  from  earlier  laminar  diffusion  flame  studies.  Comparisons  with  a  soot  formation  model  which 
incorporated  detailed  chemistry  effects  showed  good  agreement  in  terms  of  predicted  and  measured  species 
concentration  and  soot  particle  field  evolution.  In  addition,  a  novel  technique  for  measuring  soot  volume 
fraction  was  developed  based  on  laser-induced  incandescence  and  was  successfully  applied  to  similar 
laminar  diffusion  flame  studies.  This  technique  was  extended  to  droplet  and  turbulent  diffusion  flame 
conditions  where  a  two-dimensional  imaging  approach  was  employed  to  measure  soot  volume  fraction. 

The  development  of  laser-induced  incandescence  as  a  diagnostic  for  soot  represents  a  major  development  in 
terms  of  present  capability  to  make  soot  volume  fraction  measurements  in  unsteady  inhomogeneous 
combusting  flows.  Finally,  the  complete  data  set  from  these  studies  was  assembled  in  a  form  suitable  for 
dissemination  on  computer  diskettes  throughout  the  research  community  for  comparison  with  modeling 
efforts. 
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1.0  Research  Objectives 

The  objective  of  this  research  effort  was  to  investigate  the  fundamental  phenomena  controlling  soot 
particle  formation  and  destruction  in  combustion  systems.  The  emphasis  and  approach  chosen  for  the 
present  studies  was  based  on  previous  results  developed  under  AFOSR  support  which  have  proven  the 
soundness  of  the  methodology  being  utilized  [1-10].  The  present  studies  were  conducted  in  a  series  of  well- 
characterized  laminar  diffusion  flames  in  which  fuel  constituents,  temperature,  concentration  and  transport 
processes  are  systematically  varied  and/or  measured.  A  laminar  diffusion  flame  environment  was  selected 
because  it  provides  a  relatively  simple  environment  in  which  mixing  processes  are  important  while  allowing 
for  a  wide  range  of  measurement  techniques  to  be  employed.  The  importance  of  mixing  in  practical 
combustor  situations  is  well  recognized  and,  thus,  it  is  critical  to  study  soot  formation  processes  under 
conditions  where  mixing  has  a  significant  role. 

Specifically,  the  effort  emphasized  detailed  measurements  of  the  preparticle  chemistry,  particle 
inception  and  surface  growth  phenomena  important  in  soot  particle  formation.  These  studies  yielded 
quantitative  measurements  of  the  precursor  and  surface  growth  species  which  contribute  to  soot  particle 
formation  and  growth.  Data,  which  is  temporally  and  spatially  resolved,  are  of  great  value  to  modeling 
studies  as  past  studies  have  demonstrated  [1 1-14].  In  fact,  the  availability  of  such  results  has  allowed 
significant  progress  to  be  achieved  over  the  past  decade  in  soot  modeling  capabilities. 

Efforts  emphasized  measurements  of  species  concentration,  especially  acetylene  and  benzene 
concentrations,  throughout  the  flame  using  both  mass  spectrometry  and  gas  chromatography.  Analysis  has 
concentrated  on  explaining  the  soot  surface  growth  process  with  respect  to  particle  aging,  as  well  as 
collaborating  in  a  soot  modeling  effort  with  Prof.  Ian  Kennedy  of  the  University  of  California. 

Additionally,  a  new  approach  for  obtaining  quantitative  soot  volume  fraction  measurements  using  laser- 
induced  incandescence  was  developed  and  demonstrated.  In  the  following  sections,  the  techniques  used  and 
results  obtained  are  summarized.  More  detailed  discussions  are  included  in  the  appendices  for  most  of  the 
results.  For  this  reason,  the  main  body  of  the  report  will  strive  for  concise  presentations  of  the  most 
important  results. 

2.0  Research  Approach 

For  the  present  studies,  a  coannular  laminar  diffusion  flame  was  selected  for  study.  This  flame 
configuration  has  been  widely  studied  and  has  been  the  focus  of  extensive  fuel  molecular  structure  and 
pressure  sensitivity  studies  in  our  laboratory  [1-10].  Since  practical  combustors  typically  involve  the 
separate  injection  of  fuel  and  oxidizer,  laminar  diffusion  flames  provide  a  simple,  tractable  environment  in 
which  to  study  soot  formation  while  maintaining  conditions  where  mixing  effects  are  important.  Soot 
formation  can  be  viewed  as  a  series  of  sequential  processes  which  occur  to  some  degree  in  all  combustion 
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environments  where  soot  is  observed.  It  is  the  relative  importance  of  these  processes  and  their  coupling 
which  determines  the  amount  of  soot  formed,  its  distribution  throughout  the  combustion  zone  and  its 
eventual  emission  from  a  specific  combustor.  Essential  to  the  study  of  these  processes  is  the  selection  of  an 
experimental  configuration  which  can  provide  the  range  of  conditions  typical  of  fuel-rich  combustion  zones. 
Laminar  coannular  diffusion  flame  studies  have  clearly  demonstrated  that  such  conditions  can  be  achieved 
and  studied  [5].  Furthermore,  such  studies  are  proving  to  be  increasingly  valuable  for  the  validation  of 
modeling  studies.  With  these  points  in  mind,  the  coannular  laminar  diffusion  flame  represented  a  highly 
attractive  experimental  configuration. 

The  soot  formation  studies  have  been  conducted  in  an  atmospheric  diffusion  flame  facility 
consisting  of  a  coannular  diffusion  flame  burner,  burner  chimney,  positioning  system  and  gas  metering 
system.  The  burner  has  a  configuration  consisting  of  a  1.1  cm  fuel  tube  surrounded  by  a  10  cm  air 
annulus.  The  air  passage  is  partially  filled  with  glass  beads  followed  by  a  series  of  fine  screens  to  provide 
flow  conditioning.  A  ceramic  honeycomb  2.54  cm  in  thickness  is  used  at  the  exit  to  provide  a  uniform  flow 
field.  The  fuel  tube  which  extends  4.8  mm  above  the  ceramic  honeycomb  also  is  partly  filled  with  glass 
beads  to  condition  the  flow.  The  fuel  flow  can  consist  of  up  to  three  gases,  each  metered  with  a  separate 
rotameter.  This  allows  for  mixtures  of  fuels  as  well  as  nitrogen  dilution  of  the  fuel  for  temperature  and 
concentration  variation.  The  air  flow  is  metered  using  a  mass  flow  meter  or  rotameter  which  can  monitor 
flows  up  to  5  SCFM  of  air.  To  protect  the  flame  from  room  disturbances,  a  metal  chimney  has  been 
incorporated  into  the  burner  facility.  The  chimney  translates  horizontally  with  the  burner  while  the  burner 
slides  vertically  within  the  chimney.  Slots  machined  in  the  chimney  provide  access  for  the  optical  and 
intrusive  probing  of  the  flame. 

A  variety  of  diagnostic  techniques  were  applied  in  these  studies  to  measure  the  soot  particle, 
temperature,  species  concentration  and  velocity  fields.  Both  non-intrusive  laser-based  techniques,  such  as 
laser  light  scattering,  laser  velocimetry,  laser-induced  incandescence  and  laser-induced  fluorescence,  and 
intrusive  probing  followed  by  mass  spectrometric  and  gas  chromatographic  analysis  were  utilized.  A  brief 
description  of  the  apparatus  used  for  these  measurements  follows  for  the  reader's  information. 

2.1  Pulse  Laser  Systems 

An  Nd-YAG  pulsed-dye  laser  system  was  utilized  for  fluorescence  detection  of  OH  radicals  as  well 
as  for  PAH  species  measurements  while  the  Nd-YAG  laser  alone  provides  the  light  source  for  the  laser- 
induced  incandescence  measurements.  The  major  components  of  the  system  involve  the  dye  laser,  which 
consists  of  an  Nd-YAG  pulsed  laser  operating  at  532  nm.  a  dye  laser  module  and  a  doubling  crystal.  The 
present  system  has  the  capability  to  tune  over  a  wavelength  range  extending  from  220  to  600  nm.  This 
wavelength  range  is  quite  suitable  for  OH  fluorescence  as  well  as  PAH  measurements.  Detection  is 
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achieved  using  a  photomultiplier  or  an  intensified  CCD  array  camera.  In  the  case  of  the  photomultiplier 
(PMT),  the  PMT  signal  is  input  to  a  box  car  integrator  to  provide  for  signal  averaging  and  background 
noise  rejection.  For  the  intensified  CCD  array  camera,  direct  storage  of  the  images  on  a  personal  computer 
is  used  to  record  the  data.  This  unit  has  a  50  ns  gate  capability  and  can  be  used  for  signal  averaging  of 
multiple  laser  shots.  Using  this  camera,  planar  laser  imaging  is  also  possible,  which  allows  spatially 
extensive  measurements  to  be  obtained. 

2.2  Intrusive  Probe  Sampling  for  Species  Concentration  Measurements 

In  order  to  provide  more  quantitative  information  on  the  species  important  in  the  formation  and 
growth  of  soot  particles,  a  mass  spectrometer  system  was  assembled  for  these  studies.  Gas  samples  from 
the  flames  under  study  were  obtained  using  a  quartz  microprobe  technique  similar  to  the  approach  taken  by 
Smyth  et  al.  [15].  Measurements  of  gas  phase  species  in  the  soot  precursor  and  particle  inception  regions 
can  be  obtained  using  standard  sampling  probe  configurations.  Measurements  in  the  present  studies  were 
also  extended  into  the  surface  growth  region  in  order  to  ascertain  the  important  species  responsible  for 
mass  addition  through  surface  growth.  As  part  of  this  research  program,  approaches  allowing  sampling  of 
fuel-rich  particle  laden  regions  had  to  be  developed  to  allow  measurements  in  the  post-particle  inception 
region  of  the  flame.  Sampling  probe  techniques  usually  are  limited  by  orifice  clogging  problems  which  can 
only  be  overcome  by  increasing  the  orifice  dimensions  and  sacrificing  the  spatial  resolution  of  the 
measurements.  To  overcome  the  aforementioned  orifice  clogging  problem,  a  novel  sonic  sampling  probe 
has  been  developed  and  applied  to  a  series  of  laminar  diffusion  flames  containing  various  amounts  of  soot. 
Since  the  details  of  the  design  and  operation  of  this  probe  were  previously  presented  [16],  only  a  brief 
description  is  provided  below. 

The  approach,  employed  to  prevent  clogging  of  the  probe  orifice,  involves  mechanically  oscillating 
a  wire  through  the  orifice  region  using  a  spring-loaded  solenoid  plunger  whose  driving  circuit  is 
electronically  interrupted  periodically  (see  Figures  1  and  2).  The  constant  motion  of  the  wire  relative  to  the 
quartz  tube  is  responsible  for  keeping  the  orifice  open.  For  the  present  probe,  the  effective  orifice  is  the 
annulus  formed  by  an  oscillating  130  |i.m  diameter  Sapphire  fiber  of  uniform  cross  section  and  the  160  (am 
diameter  orifice  in  the  quartz  tube.  Therefore,  the  effective  orifice  is  equivalent  to  approximately  a  100  |im 
circular  hole.  The  Sapphire  fiber  is  attached  to  a  1.59  mm  tube  which  is  connected  to  a  spring-loaded  soft 
iron  solenoid  plunger  which  oscillates  in  the  orifice.  This  probe  is  referred  to  as  the  electromechanical 
sonic  probe  (EMS)  since  it  employs  a  sonic  sampling  orifice  with  an  electromagnetic  mechanism  to 
maintain  the  integrity  of  the  orifice  opening. 
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Figure  1.  Schematic  representation  of  the 


Figure  2.  Schematic  representation  of  the  quadrapole  mass  spectrometer  and  probe  assembly 


The  mass  spectrometer  used  for  these  studies  utilizes  a  quadrupole  mass  filter  with  a  mass  range  of 
1  to  500  amu.  This  unit  (an  Extrel  Model  C-50)  has  excellent  mass  resolution  (0.1  amu)  and  high 
sensitivity.  A  gas  chromatograph  (Hewlett-Packard  5860)  is  also  used  for  gas  analysis  to  complement  the 
mass  spectrometer  measurements.  This  system  utilizes  a  specially  designed  sampling  and  pressurization 
approach  to  allow  the  EMS  probe  to  be  used.  Radial  profiles  of  the  species  concentration  as  a  function  of 
axial  position  are  obtained  by  traversing  the  flame  in  the  radial  direction.  When  these  measurements  are 
combined  with  the  soot  particle  and  velocity  measurements,  the  temporal  evolution  of  the  species 
concentration  and  its  relationship  to  soot  formation  can  be  investigated. 

3.0  Research  Accomplishments 

The  research  accomplishments  of  the  grant  are  summarized  in  the  following  sections.  Specifically, 
three  areas  of  activities  are  described: 

(1)  soot  particle  surface  growth  mechanisms,  (2)  soot  modeling  studies  and  (3)  laser-induced  incandescence 
measurements  of  soot  volume  fraction.  Detailed  discussions  of  these  accomplishments  are  also  provided  in 
the  appendices,  as  noted.  In  addition,  an  extensive  review  of  the  current  understanding  of  soot  formation  is 
presented  in  Appendix  A  for  the  interested  reader. 

3.1  Soot  Particle  Surface  Growth  Mechanisms 

Extensive  spatially  resolved  measurements  of  the  species  concentration  field  have  been  obtained 
with  both  mass  spectrometry  and  gas  chromatography  analysis  approaches  for  a  well-studied  ethene/air 
laminar  diffusion  flame.  The  complete  set  of  mass  spectrometry  measurements  obtained  previously  has 
been  repeated  to  allow  for  an  estimation  of  the  accuracy  and  precision  of  the  measured  data.  This 
information  is  critical  for  modeling  comparisons  in  order  to  accurately  assess  when  real  differences  exist 
between  measurement  and  prediction.  Additionally,  similar  spatially  resolved  concentration  measurements 
were  obtained  with  a  gas  chromatography  approach  to  provide  an  independent  measurement  of  the 
concentrations  for  comparison  with  the  mass  spectrometer  data.  The  gas  chromatograph  also  allowed 
determinations  of  the  concentrations  of  some  key  species  not  easily  measured  with  the  mass  spectrometer, 
in  particular  carbon  monoxide. 

The  major  conclusions  of  these  studies  were: 

1 .  Soot  particle  surface  growth  terminates  because  of  the  depletion  of  the  hydrocarbon  species 

responsible  for  surface  reactions  which  add  mass  to  the  particle.  In  particular,  the  measurements 
of  acetylene  and  benzene  clearly  show  these  species  decrease  to  trace  concentration  levels  at  the 
peak  of  the  soot  volume  fraction  location. 
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2.  Comparisons  with  results  describing  the  aging  process,  which  soot  particles  have  been  observed  to 
experience  under  premixed  flame  conditions,  show  that  this  mechanism  has  only  a  small 
contribution  to  the  reduction  in  the  surface  growth  rate  for  these  laminar  flame  conditions. 

These  conclusions  are  consistent  with  recent  work  by  Honnery  and  Kent  [17,18],  who  have  shown  that  in 
long  laminar  diffusion  flames,  soot  growth  can  be  extended  to  very  long  residence  times.  A  more  detailed 
discussion  of  these  results  is  contained  in  Appendix  B,  where  the  data  and  analysis  to  support  the  above 
conclusions  are  presented. 

3.2  Soot  Modeling  Studies 

A  collaborative  effort  was  initiated  with  Prof.  Ian  Kennedy  of  the  University  of  California  at  Davis 
to  compare  the  measured  soot  and  species  concentrations  fields  with  that  predicted  by  a  model  developed 
by  Prof.  Kennedy  and  his  co-workers.  A  system  of  elementary  reactions  were  used  to  describe  the  gas 
phase  Ci  and  C2  chemistry.  The  model  incorporates  a  simple  description  of  the  growth  of  soot  particles 
which  performs  a  good  job  of  reproducing  the  amount  of  soot  in  the  flame  and  the  transition  from  non¬ 
sooting  to  sooting  conditions.  Stable  gas  phase  species  measured  in  the  flame  along  with  OH 
concentrations  previously  obtained  [10]  were  compared  with  model  predictions.  The  model  reproduced 
most  of  the  specified  concentration  profiles  quite  well.  It  was  found  that  the  assumption  of  equilibrium 
concentrations  for  the  minor  radical  species  was  grossly  in  error,  particularly  near  the  end  of  the  flame 
where  temperatures  were  low.  Calculated  radical  species  were  found  in  concentrations  that  were  orders  of 
magnitude  greater  than  equilibrium.  Comparisons  with  measured  OH  concentrations  supported  this 
observation.  A  more  detailed  discussion  of  these  results  is  contained  in  Appendix  C,  which  has  recently 
been  submitted  to  Combustion  and  Flame  for  publication. 
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3.3  Laser-Induced  Incandescence  Measurements  of  Soot  Volume  Fraction 

Laser-induced  incandescence  has  been  used  to  obtain  spatially-resolved  measurements  of  soot 
volume  fraction  in  a  laminar  diffusion  flame,  in  which  comparisons  with  laser  scattering/extinction  data 
yield  excellent  agreement.  In  addition,  the  laser-induced  incandescence  signal  is  observed  to  involve  a  rapid 
rise  in  intensity  followed  by  a  relatively  long  (ca.  600  ns)  decay  period  subsequent  to  the  laser  pulse,  while 
the  effect  of  laser  fluence  is  manifest  in  non-linear  and  near-saturated  response  of  the  laser-induced 
incandescence  signal  with  the  transition  occurring  at  a  laser  fluence  of  approximately  1.2  x  108  W/cm2. 
Spectral  response  of  the  laser-induced  incandescence  involves  a  continuous  spectrum  in  the  visible 
wavelength  range  due  to  the  blackbody  nature  of  the  emission.  Simultaneous  measurements  of  laser- 
induced  incandescence  and  light-scattering  yield  encouraging  results  concerning  the  mean  soot  particle 
diameter  and  number  concentration.  Thus,  laser-induced  incandescence  can  be  used  as  an  instantaneous, 
spatially-resolved  diagnostic  of  soot  volume  fraction  without  the  need  for  the  conventional  line-of-sight 
laser  extinction  method,  while  potential  applications  in  two-dimensional  imaging  and  simultaneous 
measurements  of  laser-induced  incandescence  and  light-scattering  to  generate  a  complete  soot  property 
characterization  are  significant.  These  results  are  presented  in  detail  in  Appendix  D. 

Recently,  the  laser-induced  incandescence  technique  was  used  to  make  novel  planar  measurements 
of  soot  volume  fraction  within  turbulent  diffusion  flames  and  droplet  flames.  The  two-dimensional  imaging 
technique  was  developed  and  assessed  by  systematic  experiments  in  a  coannular  laminar  diffusion  flame,  in 
which  the  soot  characteristics  have  been  well  established.  With  a  single-point  calibration  procedure, 
agreement  to  within  10%  was  found  between  the  values  of  soot  volume  fraction  measured  by  this  technique 
and  those  determined  by  conventional  laser  scattering-extinction  methods  in  the  flame.  As  a  demonstration 
of  the  wide  range  of  applicability  of  the  technique,  soot  volume  fraction  images  were  also  obtained  from 
both  turbulent  ethene  diffusion  flames  and  from  a  freely  falling  droplet  flame  that  bums  the  mixture  of  75% 
benzene  and  25%  methanol.  For  the  turbulent  diffusion  flames,  approximately  an  80%  reduction  in  soot 
volume  fraction  was  found  when  the  Reynolds  number  of  the  fuel  jet  increased  from  4000  to  8000.  In  the 
droplet  flame  case,  the  distribution  of  soot  field  was  found  to  be  similar  to  that  observed  in  coannular 
laminar  diffusion  flames.  The  development  of  the  laser-induced  technique,  therefore,  is  a  significant  new 
contribution  to  the  array  of  laser-based  diagnostics  available  for  combustion  measurements.  The  results  of 
studies  involving  extensions  of  the  technique  to  planar  measurements  are  given  in  Appendix  E. 
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4.0  Conclusions 

The  major  conclusions  of  studies  undertaken  as  part  of  this  contract  are: 

1 .  Soot  particle  surface  growth  terminates  because  of  the  depletion  of  the  hydrocarbon  species 
responsible  for  surface  reactions  which  add  mass  to  the  particle.  In  particular,  the  measurements 
of  acetylene  and  benzene  clearly  show  these  species  decrease  to  trace  concentration  levels  at  the 
peak  of  the  soot  volume  fraction  location.  Additional  measurements  using  gas  chromatography 
analysis  have  provided  additional  species  concentration  information  and  allowed  an  estimate  of  the 
accuracy  and  precision  of  the  measurements.  The  availability  of  such  an  estimate  was  viewed  as 
critical  to  meaningful  comparison  with  soot  chemistry  model  predictions. 

2.  Comparisons  with  results  describing  the  aging  process,  which  soot  particles  have  been  observed  to 
experience  under  premixed  flame  conditions,  show  that  this  mechanism  has  only  a  small 
contribution  to  the  reduction  in  the  surface  growth  rate  for  these  laminar  flame  conditions. 

3.  A  current  detailed  model  of  soot  formation  in  laminar  diffusion  flames  was  shown  to  be  capable  of 
reproducing  the  major  gas  phase  species  and  soot  particle  fields  reasonably  well  (i.e.  to  within  a 
factor  of  two).  Better  agreement  appears  possible  with  further  refinements  in  the  model.  Of 
particular  note  is  that  both  measurement  and  prediction  show  that  radical  species  in  these  soot 
laden  diffusion  flames  are  far  from  equilibrium.  This  observation  is  critical  to  the  development  of 
appropriate  combustion  models  for  soot  laden  flows. 

4.  A  novel  quantitative  laser  diagnostic  approach  has  been  developed  and  applied  to  the  measurement 
of  soot  volume  fraction  in  laminar  diffusion  flames  based  laser-induced  incandescence.  This 
diagnostic  approach  has  been  validated  using  well  known  soot  particle  measurements  previously 
obtained  in  laminar  diffusion  flames.  The  technique  has  recently  been  extended  to  turbulent  flame, 
droplet  burning  and  Diesel  engine  experiments. 
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SOOT  FORMATION:  A  VIEW  OF  OUR  CURRENT  UNDERSTANDING 
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Propulsion  Engineering  Research  Center 
and 

Department  of  Mechanical  Engineering 
The  Pennsylvania  State  University 
University  Park,  PA  16802-2320 


INTRODUCTION 

Over  the  past  three  decades  a  significant  amount  of  research  has  been  directed  to  understanding,  from  both  a 
fundamental  and  applied  viewpoint,  the  processes  responsible  for  soot  formation  and  destruction  in  combustion  systems. 
Several  excellent  reviews  have  been  published  which  describe  much  of  this  previous  work  [1-6].  The  basic  sequence  of 
events  which  compose  the  formation  and  destruction  processes  for  soot  particles  in  combustion  systems  is  now  well 
established.  This  sequence  includes:  1)  a  chemically  kinetic  controlled  reaction  sequence  which  results  in  the  formation  of 
the  precursor  species  needed  to  form  the  first  particles,  2)  a  particle  inception  stage  which  results  in  the  formation  of  large 
numbers  of  small  primary  particles,  3)  a  particle  growth  region  in  which  surface  growth  and  particle  coagulation  processes 
contribute  to  the  increase  in  particle  size,  and,  finally,  4)  a  stage  in  which  material  is  no  longer  added  to  the  soot  particles, 
and  size  and  number  concentration  are  controlled  by  agglomeration  or  even  are  reduced  through  oxidative  attack  [7] .  This 
sequence  of  events  describes  soot  formation  in  a  wide  variety  of  combustion  systems  including  laboratory-scale  premixed 
and  diffusion  flames  as  well  as  practical  combustors  such  Diesel  and  gas  turbine  engines.  Differences  observed  in  various 
combustion  environments  are  a  result  of  variations  in  the  importance  of  the  individual  steps  under  the  actual  operating 
conditions.  Thus,  an  understanding  of  soot  formation  and  destruction  not  only  requires  a  knowledge  of  the  individual  steps 
involved,  but  also  how  they  are  coupled. 

Below  a  review  of  the  soot  formation  and  destruction  is  presented  in  four  sections  dealing  with  soot  particle 
precursor  chemistry,  soot  particle  inception,  surface  growth  and  soot  oxidation.  The  studies  reviewed  emphasize  simple 
combustion  systems  in  which  suitable  conditions  can  be  achieved  which  allow  differentiation  of  the  various  aspects  of  the 
problem.  A  variety  of  experimental  apparatus  have  been  employed  for  such  studies  including  shock  tubes  [9,10],  laminar 
premixed  [11-20]  and  diffusion  [21-46]  flames,  and  turbulent  flames  [47-53].  Key  to  these  laboratory-scale  studies  has  been 
the  development  of  a  variety  of  diagnostic  capabilities,  both  for  gas  phase  species  and  soot  particle  measurements.  Most 
prominent  among  these  techniques  have  been  mass  spectrometric  sampling  methods  for  species  concentration  measurements 
[12]  and  laser  light  scattering  techniques  for  soot  particle  characterization  [15].  Of  particular  importance  has  been  the 
capability  to  obtain  spatially  resolved  measurements  of  key  species  and  soot  particle  concentrations  [12,  23-25]  to  provide 
mechanistic  information  on  the  evolution  of  the  soot  precursor  and  particle  fields. 

An  additional  development  which  has  had  an  important  impact  is  the  application  of  large  detailed  numerical  models 
representing  the  entire  soot  formation  and  destruction  process  [54-71].  The  development  of  these  models  has  helped  to  focus 
research  in  several  areas  including  chemical  precursor  mechanisms,  particle  nucleation,  surface  growth  phenomena,  and 
oxidation  processes.  Although  differences  in  interpretation  and  understanding  remain,  a  coherent  modeling  process  is 
emerging  from  which  to  examine  likely  routes  to  soot  particle  formation. 

SOOT  PARTICLE  PRECURSOR  CHEMISTRY 

Any  understanding  of  the  soot  formation  process  must  begin  with  the  reactions  involving  the  parent  hydrocarbon 
fuel.  These  typically  are  relatively  small  hydrocarbon  species  which  undergo  pyrolysis  and  oxidative  reaction  steps  to  form 
intermediate  radical  species  that  react  with  other  hydrocarbons  present  to  form  larger  species.  It  is  useful  in  this  endeavor 
to  concentrate  on  two  distinct  groups  of  hydrocarbon  reactions,  following  the  work  of  Frenklach  and  co-workers  [54-59]. 
The  first  group  involves  reactions  which  lead  to  formation  of  the  first  and  second  aromatic  ring  species.  While  the  second 
group  involves  the  reactions  which  lead  to  subsequent  multiple  ring  species,  i.e.,  polycyclic  aromatic  hydrocarbons  (PAH). 
The  emphasis  of  aromatic  hydrocarbons  is  a  reflection  of  the  established  fact  that  soot  is  composed  of  species  with  an 
aromatic  ring  structure  [1].  Consequently,  at  some  point  in  the  soot  formation  process  these  aromatic  species  must  form. 
Currently,  data  supports  the  view  that  the  gas  phase  precursors  to  soot  particles  are  PAH  species.  Although  the  description 
here  emphasizes  a  neutral  species  mechanism,  other  studies  have  examined  an  ion  molecule  reaction  pathway  to  soot  [5, 18,72] . 


A  key  feature  of  any  mechanism  explaining  soot  formation  must  account  for  the  thermodynamic  stability  of  the 
intermediate  hydrocarbon  species.  Since  soot  formation  occurs  under  high  temperature  conditions,  one  would  expect  larger 
molecular  species  to  dissociate  unless  they  are  thermodynamically  stable.  Additionally,  the  reactions  leading  to  these  larger 
hydrocarbon  species  should  be  fast,  in  order  to  account  for  the  observed  rapid  formation  of  soot  particles  in  combustion 
systems.  In  some  sense,  the  conditions  of  rapid  reaction  rates  and  thermodynamically  stable  molecules  pose  contradictory 
requirements  which  must  be  resolved  within  the  model. 

Frenklach  et  al.  [54,55]  proposed  a  mechanism  in  which  reactions  between  acetylene  (C2H2)  and  key  radicals,  such 
as  vinyl  (C2H3),  vinylacetylene  (C4H4)  and  butadienyl  (C4H5),  lead  to  the  formation  of  the  first  aromatic  ring  through  the 
cyclization  of  QHj.  Cole  et  al.  [20]  proposed  a  slightly  different  reaction  route  involving  QH?  to  form  benzene  (QH6) 
which  Colket  [60]  used  to  explain  the  temperature  sensitivity  observed  in  shock  tube  studies  of  C2H2  at  lower  temperatures 
than  the  original  work  of  Frenklach  et  al.  [54].  Frenklach  et  al.  [55]  later  showed  that  the  difference  in  sooting  tendency 
of  various  fuels  was  due  to  the  sensitivity  of  these  radical  species  concentrations  to  fuel  molecular  structure  and  temperature. 
Other  routes  to  the  formation  of  the  first  ring  have  also  been  examined  including  the  reaction  involving  the  propargyl  radical 
(C3H3)  to  form  1,5  hexadiyne  which  rapidly  is  converted  to  benzene  [73].  Although  details  remain  to  be  resolved,  significant 
progress  has  been  made  in  achieving  a  basic  framework  for  understanding  this  important  first  step  in  the  soot  formation 
process. 

The  formation  of  the  second  ring,  begins  with  the  QH5  radicals  reacting  with  C2H2,  followed  by  a  subsequent 
reaction  with  H-atom  to  reform  a  radical  species.  This  radical  again  reacts  with  C2H2  and  undergoes  ring  closure  to  form 
the  naphthalene  radical  [54].  Acetylene,  which  is  present  in  large  concentrations,  thus,  represents  a  primary  building  block 
in  the  growth  of  large  hydrocarbon  species.  Hydrogen  atom  reaction  steps  are  key  to  recreating  the  radical  species,  and  are 
continually  necessary  in  the  subsequent  steps  to  the  formation  of  multiple  ring  PAH  species.  Since  all  these  reactions  are 
reversible,  dissociation  of  the  products  to  reactants  is  likely  to  some  degree,  thus  both  the  forward  and  reverse  reactions  are 
important. 

A  point  to  emphasize  is  that  fuel  molecular  structure  is  likely  important  in  establishing  the  concentrations  of 
acetylene  and  important  QHX  and  C4Hy  hydrocarbon  radicals.  Furthermore,  modeling  studies  have  shown  that  these  reactions 
are  very  sensitive  to  temperature  in  the  range  around  where  soot  is  observed  to  rapidly  form  (1300K<T<  1800K).  For 
aromatic  fuels  some  support  has  been  given  for  a  more  direct  route  to  a  two  ring  aromatic  species  through  biphenyl  formation 
[55] ,  again  emphasizing  that  fuel  molecular  structure  likely  impacts  the  early  steps  in  the  soot  formation  process  significantly. 

Once  one  and  two  ring  compounds  are  formed  growth  to  much  larger  species  is  required.  Typically  incipient  soot 
particles  are  viewed  to  achieve  molecular  weights  between  3000  and  10,000  amu.  To  achieve  the  continued  molecular 
growth  needed,  Frenklach  et  al.  [54,55]  proposed  a  series  of  reactions  where  an  aromatic  radical  undergoes  sequential 
reactions  with  C2H2  and  H-atoms  to  form  a  subsequent  ring.  Interspersed  along  the  reaction  path  is  the  formation  of  highly 
stable  species,  such  as  acenaphthylene,  which  contain  a  five  membered  ring.  The  formation  of  these  species  is  substantially 
more  favorable  in  the  forward  direction  which  "pushes”  the  species  growth  to  larger  molecules.  Because  of  their  stability 
these  species  have  been  termed  "islands  of  stability”  and  are  key  to  providing  a  thermodynamically  allowable  path  to  soot. 
Before  leaving  this  point  it  should  be  noted  that  although  such  species  are  stable,  they  do  not  represent  termination  points 
in  the  reaction  pathway.  Subsequent  attack  by  H-atoms  provides  a  means  to  again  form  a  radical  species  which  will  react 
with  C2H2  to  eventually  lead  to  another  "island  of  stability".  The  above  described  process  answers  the  earlier  mentioned 
requirement  for  highly  stable  intermediate  species  which  can  also  react  rapidly  to  form  larger  hydrocarbons. 

Although  developments  related  to  the  details  of  individual  steps  involved  in  the  mechanism  continue,  the  breadth 
of  application  of  this  approach  in  both  full  and  reduced  reaction  schemes  is  impressive  [62-63 ,67,7 1 ,74].  The  advances  made 
in  the  last  decade  in  the  area  of  reaction  mechanisms  describing  soot  precursor  formation  may  be  the  most  impressive 
achievement  in  the  field  over  that  time  period. 

SOOT  PARTICLE  INCEPTION 

Although  the  basis  the  soot  particle  formation  is  closely  linked  to  the  formation  and  growth  of  PAH  species,  little 
direct  measurements  exists  of  the  rate  of  inception  of  soot  particles  under  conditions  where  PAH  concentration  fields  are 
known.  Furthermore,  light  scattering  techniques  use  to  measure  soot  particles  have  difficulties  yielding  quantitative  results 
due  broadband  background  fluorescence  induced  by  the  laser  itself.  Consequently,  in  the  author's  opinion,  soot  inception 
is  an  area  where  further  quantitative  results  are  needed  to  further  clarify  the  phenomena  resulting  in  soot  particles.  However, 
much  has  been  learned  over  the  past  decade  in  this  area  as  well.  It  is  now  clear  that  the  formation  of  soot  particles  is  not 
due  to  a  classical  nucleation  event,  but  is  a  result  of  continual  growth  of  hydrocarbon  precursors  to  soot.  Recent  work  has 
more  clearly  linked  the  decrease  in  larger  PAH  species  with  the  appearance  of  soot  particles  [75-77],  qualitatively 


demonstrating  the  relationship  between  PAH  production  and  soot  particle  formation.  These  studies  of  particle  inception 
conducted  under  premixed  flame  conditions  are  being  complemented  by  laminar  diffusion  flame  studies  where  soot  particle 
inception  is  directly  measured  [781.  Previous  qualitative  measurements  indicated  that  the  first  detected  soot  particles  in 
laminar  ethene  flames  occur  in  the  region  of  maximum  observed  PAH  fluorescence,  but  no  quantitative  inception  rate  data 
was  obtained  [24].  Recently  based  on  the  limited  experimental  measurements  related  to  inception,  simple  nucleation  rate 
expressions  are  being  developed  [68]  and  compared  with  measurements  where  possible  [78].  These  initial  developments  show 
some  potential  for  success  and  should  be  followed  closely  in  the  future  as  more  data  become  available. 

Despite  the  lack  of  detailed  experimentally  measured  inception  rates,  modeling  efforts  have  met  with  considerable 
success  in  representing  inception  as  a  coagulation  process  with  a  critical  molecular  size  identified  with  the  incipient  soot 
particles  [57,58,62,66,71,74,75,79].  McKinnon  and  Howard  have  used  the  term  "reactive  coagulation"  to  describe  their 
work  in  this  area  to  indicate  that  the  "sticking  collisions  between  PAH  are  stabilized  by  the  formation  of  a  chemical  bond" 
[75].  Recent  modeling  work  has  also  considered  the  relative  importance  of  C2H2  and  PAH  collision  partners  on  the  growth 
of  soot  precursors  and  incipient  particles  and  identified  under  what  conditions  these  species  contribute  to  precursor  formation, 
as  well  as  surface  growth  process  [57,58,62,66,71,74,75]. 

SURFACE  GROWTH 

Although  particle  inception  is  one  of  the  key  processes  in  the  formation  of  soot,  it  actually  accounts  for  only  a  small 
fraction  of  the  mass  addition.  The  major  mass  addition  process  involves  the  surface  growth  stage  which  follows  inception. 
A  great  deal  of  attention  has  been  given  to  this  process  in  both  premixed  and  diffusion  flames  [14,19,80,81,82].  Major 
contributions  to  the  present  understanding  of  surface  growth  have  been  made  by  Wagner  [13]  and  Harris  [14].  The 
mechanisms  developed  in  these  studies  differ  significantly  in  concept  and  work  is  continuing  to  further  understand  the  surface 
growth  process  [83-86]. 

Most  approaches  to  soot  surface  growth  phenomena  argue  that  acetylene,  or  a  species  closely  associated  with 
acetylene,  provides  the  major  surface  reactant.  Simple  mass  balance  considerations  provide  the  basis  for  this  observation, 
since  only  acetylene  appears  to  be  present  in  sufficient  concentration  to  account  for  the  amount  of  mass  added  to  the  soot 
particles  [14].  However,  other  workers  have  pointed  out  that  under  suitable  conditions  large  PAH  species  may  directly  add 
to  the  growing  soot  particles  [57,61,66,75]. 

Recently  there  has  been  increased  interest  in  the  specific  phenomena  controlling  the  surface  reactivity  in  soot  surface 
growth  processes  [57,61,66,83-86].  Each  of  these  studies  has  focussed  on  the  concept  of  active  surface  sites  as  reaction 
centers  for  the  addition  of  carbon  through  reactions  with  gas  phase  species  such  as  acetylene.  The  objectives  of  these  studies 
differ,  but  they  do  raise  a  set  of  common  features.  First,  soot  particles  possess  highly  reactive  surfaces,  exceeding  by  two 
orders  of  magnitude  the  reactivity  of  pyrolytically  formed  carbon  particles  (pyrocarbon)  [85].  The  regeneration  of  active 
sites  appears  to  be  related  to  atomic  hydrogen  reactions  with  the  surface  to  form  new  reactive  sites  [57,61,66,79],  which  have 
a  relativity  constant  concentration  in  many  situations  [79,85].  This  approach  has  been  used  to  reconcile  the  disparate 
observations  in  the  literature  regarding  the  sensitivity  of  the  surface  growth  process  to  the  available  soot  particle  surface  area. 
The  results  of  Harris  and  Weiner  [14,19]  indicate  a  direct  proportionality  to  available  surface  areas  as  indicated  by: 


dm 

dt 


k  S[CJiy] 


where  m  is  the  soot  mass,  k  is  the  growth  rate  constant,  S  is  the  surface  area,  and  CxHy  represents  the  growth  species,  e.g. 
C2H2. 

Results  from  the  work  of  Haynes  et  al.  [86]  and  Weischnowsky  et  al.  [16]  indicate  that  the  surface  growth  process 
is  independent  of  the  available  surface  area.  In  the  recent  work,  Harris  [79],  Haynes  [85]  and  Frenklach  [57,58]  have 
replaced  or  augmented  the  surface  area  arguments  with  the  addition  of  a  surface  active  site  concentration.  Although  the 
details  of  the  representations  are  different,  the  basic  approach  involves  similar  concepts  with  regard  to  the  importance  of  these 
surface  sites.  For  example,  Haynes  [85]  represents  the  increase  in  soot  volume  fraction  with  time  as: 


d£v 

dt 


where  fv  is  the  soot  volume  fraction,  a  is  the  growth  reaction  probability  on  collision,  Ng  is  the  number  density  of  gas  phase 
growth  species,  vg  is  the  molecular  volume  of  the  growth  species  when  incorporated  in  the  soot  particle,  c  is  the  mean  free 
speed  of  the  gas  phase  growth  species  (SkT/Tmig)172,  Ns  is  the  total  concentration  of  active  sites  and  A.  is  the  active  area  of 
a  single  site.  In  this  representation  the  difference  in  the  experimental  results  is  determined  by  the  behavior  of  the  active  sites 
in  the  various  flames.  In  the  results  of  Harris  and  Weiner  [14, 17],  the  active  site  concentration  remains  nearly  constant  and, 
thus,  scales  with  particle  surface  area.  In  contrast,  for  the  higher  temperature  flame  studies  by  Haynes  et  al.  [86]  and 
Weischnowsky  [16],  the  surface  site  concentration  is  argued  to  decrease  with  time.  The  basis  for  the  decay  in  the  active 
surface  site  is  not  known,  although  the  higher  temperature  flames  can  be  argued  to  provide  an  environment  in  which  more 
rapid  annealing  of  the  soot  particles  may  be  present.  Frenklach  [57]  has  argued  in  a  somewhat  similar  manner  that  the  state 
of  the  initial  surfaces  of  the  soot  particles  may  differ  in  the  two  flames  with  the  higher  temperature  flames  having  fewer 
active  sites. 

Much  of  the  above  work  emphasizes  detailed  treatment  of  the  chemistry  and  physics  important  in  surface  growth. 
Other  workers  have  explored  the  potential  to  relate  surface  growth  to  more  global  properties  of  the  flames  which  may  allow 
more  universal,  simpler  representations  of  soot  formation  that  can  be  incorporated  into  more  complicated  flow  field 
calculations,  such  as  require  for  turbulent  combustion  conditions  [63,81,82,87].  In  these  approaches,  the  surface  growth 
process,  or  in  some  cases  the  actual  soot  volume  fraction  [49],  are  determined  from  measured  mixture  fraction  and 
temperature  conditions  in  diffusion  flames.  These  results  are  then  used  as  a  basis  for  predicting  surface  growth  over  a  wide 
range  of  flame  conditions.  A  major  advantage  of  this  approach  is  the  reduced  detail  which  must  be  given  to  representing 
the  chemical  kinetic  processes  involved  in  the  surface  reaction  mechanism. 

The  above  work  illustrates  the  detail  to  which  present  soot  models  and  experiments  can  consider  the  fundamental 
mechanisms  controlling  individual  processes.  Despite  the  questions  concerning  the  specifics  of  the  surface  growth 
mechanism,  recent  efforts  to  model  soot  particle  growth  in  diffusion  flames  have  met  with  some  success  in  reproducing  the 
general  features  of  the  soot  formation  process  [57-68,81,82,88].  However,  the  model  representation  for  the  surface  growth 
process  shows  a  strong  sensitivity  to  the  specific  surface  growth  rate  constant  [65].  In  this  case,  small  differences  in  the 
growth  rate  constant  could  be  argued  to  be  responsible  for  variations  in  the  sooting  propensity  of  various  fuels.  To  resolve 
such  questions  further  detailed  measurements  of  surface  growth  rates,  particularly  under  diffusion  flame  conditions  are 
needed. 

Before  leaving  the  topic  of  surface  growth,  it  is  worth  mentioning  the  importance  that  developments  in  the 
measurement  of  soot  size  and  morphology  are  having  in  this  area.  The  application  of  fractal  concepts  to  laser  light  scattering 
measurements  is  providing  more  accurate  determination  of  the  actual  soot  particle  surface  area  [7,33,34,41,89-92].  It  has 
been  well  established  that  soot  particle  size  is  affected  by  coagulation  and  classical  coagulation  models  have  been  well  suited 
to  describing  the  evolution  of  the  number  concentration  of  particles  [9],  However,  the  specific  nature  of  the  resulting 
aggregate  was  typically  argued  to  be  spherical  by  necessity  of  the  optical  analysis  employed.  Although  this  approach  yields 
reasonable  relative  trends  with  regard  to  particle  size  and  surface  area,  it  does  have  an  inherent  error  since  the  particles  are 
not  truly  spherical.  The  fractal  approaches  introduced  recently  are  providing  a  better  measurement  capability  for  the 
determination  of  soot  particle  surface  area  and  allowing  a  more  detailed  understanding  of  the  structure  of  soot  formed  in 
flames. 


SOOT  PARTICLE  OXIDATION 

Following  the  termination  of  surface  growth,  oxidation  of  the  soot  particles  can  occur  resulting  in  either  total 
destruction  or  emission  of  some  fraction  of  the  soot  particles.  To  effectively  predict  the  oxidative  destruction  of  soot 
particles,  an  understanding  of  the  reaction  rates  as  a  function  of  temperature  for  various  oxidizing  species  is  required. 
Limited  attention  has  been  given  to  this  area  recently.  An  appropriate  review  of  the  present  understanding  of  soot  particle 
oxidation  is  provided  by  Neoh  et  al.  [93].  Although  discrepancies  between  various  studies  of  soot  oxidation  have  been 
described  in  the  literature,  there  is  a  general  consensus  that  the  reactions  with  02  and  OH  are  the  dominant  reaction  pathways. 
Under  lean  conditions,  02  is  expected  to  be  the  dominant  oxidizing  species  and  the  Nagle  and  Strickland-Constable  expression 
has  been  traditionally  used  [94].  For  fuel  rich  conditions,  where  the  02  concentration  is  small,  OH  is  the  major  oxidizing 


species  whose  reactivity  has  been  quantified  by  Neoh  et  al.  [93], 

Recent  work  has  examined  reaction  rate  data  for  02  and  indicates  that  the  tradition  use  of  the  Nagle  and  Strickland- 
Constable  expression  may  lead  to  errors  in  assessing  the  importance  of  this  oxidation  route  [95].  Studies  of  the  importance 
of  OH  as  oxidation  species  have  also  raised  questions  about  the  rate  of  02  reaction  with  soot  particles  [96a,b],  Additional 
work  has  also  been  reported  on  O  and  OH  radicals  [97].  These  developments  point  to  the  need  for  continued  examination 
of  the  oxidation  stage  of  the  soot  formation  and  destruction  process. 

SUMMARY 

The  above  discussion  provides  an  overview  of  the  key  developments  which  have  occured  over  the  last  decade  in 
terms  of  increasing  our  understanding  of  soot  formation  and  destruction.  Not  all  the  important  work  could  be  reviewed  with 
the  recent  studies  of  soot  formation  in  high  pressure  flames  the  most  glaring  omission.  As  we  look  into  the  future,  we  can 
expect  to  see  increase  attention  given  to  turbulent  flame  systems  and  high  pressure  environments.  More  emphasis  is  likely 
to  focussed  on  practical  systems  such  as  gas  turbine  and  Diesel  engines  where  soot  impacts  both  performance  and  emissions 
issues. 
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ABSTRACT 


Newly  added  measurements  of  the  spatial  distribution  of  species  concentrations  in  a  ethene/air, 
coannular,  laminar  diffusion  flame  have  been  integrated  with  previous  measurements  of  soot 
and  flowfield  (parameters,  quantities)  to  establish  the  mechanism  controlling  the  cessation  of  soot 
surface  growth.  In  particular,  the  measurements  of  C2H2  and  other  hydrocarbons  have  revealed 
that  soot  surface  growth  ceases  due  to  the  depletion  of  hydrocarbons  and  not  due  to  the  loss  of  soot 
particle  reactivity.  This  observation  is  in  agreement  with  recent  work  by  Honnery  and  Kent,  who 
have  shown  that  in  long  laminar  diffusion  flames,  soot  growth  can  be  extended  to  very  long 
residence  times.  Values  for  the  specific  surface  growth  rates  in  this  diffusion  flame  are  in 
agreement  with  previous  measurements  in  both  premixed  and  diffusion  flames.  Now  that 
concentrations  of  many  intermediate  hydrocarbons,  including  C2H2  and  C6H6  are  available  in  a 
diffusion  flame,  the  rate  coefficient  for  surface  growth  by  C2H2  can  be  calculated  and  is  in 
agreement  with  values  found  in  premixed  flame  studies. 


INTRODUCTION 

The  evolution  of  soot  particles  has  been  the  object  of  study  in  a  wide  variety  of  flame  environments. 
Although  studies  conducted  in  practical  devices,  such  as  gas  turbine  combustors  or  Diesel  engines, 
have  yielded  much  useful  engineering  insight,  fundamental  understanding  of  soot  formation 
and  destruction  has  been  typically  gleaned  from  simpler  flame  studies.  Much  of  the  basic 
understanding  of  soot  formation  is  a  result  of  experimental  work  in  shock  tubes,  low  pressure  and 
atmospheric  pressure  premixed  flames,  and  diffusion  flames  [1-19].  However,  knowledge  of  the 
chemical  environment  as  it  relates  to  soot  formation  has  come  primarily  from  experiments  in 
premixed  flat  flames  [3-11]  and  been  supplemented  by  computer  models  [20-23].  These 
experimental  conditions  differ  from  many  practical  applications  which  involve  non-premixed 
conditions.  Therefore,  it  is  important  to  study  soot  formation  processes  under  conditions  where 
diffusive  and  mixing  processes  control  the  reaction  rate  (high  Damkohler  Number). 

Furthermore,  when  focusing  on  fundamental  soot  phenomena  it  is  important  to  select  as  simple  an 
environment  as  possible  while  retaining  the  major  characteristics  of  diffusion  flames.  Though 
all  practical  combustion  devices  can  have  turbulence  to  very  high  levels,  it  is  more  tractable  to 
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study  soot  formation  under  laminar  conditions  where  the  complexities  of  turbulence  do  not  obscure 
the  primary  object  of  study. 

Among  the  simplest  laminar  diffusion  flames  which  contains  many  characteristics  of  practical 
combustors  is  an  axisymmetric  system  of  gaseous  fuel  issuing  from  a  tube  with  a  larger  flow  of  air 
surrounding  it.  The  experimental  setup  is  straightforward  and  the  flames  are  quite  stable.  In  the 
present  experiment,  the  burner  configuration  and  dimensions  are  identical  to  that  used  by 
previous  researchers  to  acquire  a  significant  set  of  data  on  flow  variables  (temperature  and 
velocity)  [14],  and  on  the  soot  aerosol  (volume  fraction,  number  density  and  particle  diameter) 
[121.  Recently,  hydroxyl  radical  (OH)  mole  fractions  have  been  determined  using  laser-induced 
fluorescence  [24]  and  fluorescence  from  polycyclic  aromatic  hydrocarbons  (PAH)  have  been 
examined  [25].  Due  to  this  extensive  database,  researchers  involved  in  the  numerical  modeling  of 
the  soot  formation  and/or  chemical  kinetics  in  diffusion  flames  have  endeavored  to  model  this 
system  [26,  27] . 

The  present  work  expands  on  this  set  of  data  for  a  laminar  diffusion  flame  burning  C2H4  in  air  at 
atmospheric  pressure.  The  information  to  be  added  to  the  existing  measurements  involve  the 
spatial  distribution  of  the  molecular  species  in  the  flame.  When  combined  with  the  velocity 
measurements,  knowledge  of  the  environment  surrounding  the  soot  particles  on  a  temporal  basis 
is  gained.  The  following  discussion  will  pursue  a  more  sophisticated  analysis  of  the  data  than 
previously  possible  in  the  surface  growth  region. 

In  order  to  determine  the  distribution  of  molecular  species  throughout  the  flame,  a  method  of 
measuring  the  concentrations  is  necessary.  Although  several  non-intrusive  laser-based 
techniques  exist  for  species  measurements  (e.g.  LIF,  absorption,  CARS,  Raman)  many  are 
limited  in  the  number  of  species  that  can  be  quantitatively  measured.  Furthermore,  when  soot 
particles  are  present  interference  with  the  measurements  can  make  these  techniques  difficult  if 
not  impossible  to  implement  [28].  Therefore,  intrusive  probes,  which  have  been  used  for  many 
years,  will  be  used  to  extract  gas  samples  for  on-line  or  later  batch  analysis. 

EXPERIMENTAL 


A  laminar  ethene/air  diffusion  flame  was  established  on  a  coannular  burner  similar  in  design  to 
that  used  in  previous  studies  [12]  and  only  a  brief  description  will  be  given  here.  This  burner 
consists  of  a  1.11  cm  diameter  brass  fuel  tube  surrounded  by  a  10.2  cm  diameter  air  annulus.  A 
brass  cylindrical  chimney  was  used  to  shield  the  flame  from  laboratory  air  currents  with  access 
for  sampling  provided  by  slots  machined  in  the  chimney  wall.  The  ethene  and  air  flow  rates  were 
3.85  cm3/s  and  713  cm3/s  respectively,  which  were  measured  with  calibrated  rotameters.  Ethene 
with  a  stated  purity  of  99.5%  was  used,  while  filtered  and  dried  air  was  supplied  by  an  in-house 
compressor.  These  conditions  replicate  those  of  Santoro  and  coworkers  [12-14,  24,  25,  29]  so  that  the 
current  and  previous  measurements  can  be  brought  together  synergistically. 

The  extraction  of  gas  samples  is  accomplished  throughout  the  flame  using  an  intrusive  quartz 
microprobe  sampling  approach.  Quartz  microprobes  have  been  used  in  many  studies  to  perform  the 
extraction  of  gas  samples  from  premixed  flames  [4,  11,  30,  31]  and  diffusion  flames  [32-36].  Until 
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recently  their  use  has  been  limited  to  regions  where  the  soot  volume  fraction  is  low  to  avoid 
clogging  the  small  probe  orifice  (—100  um).  A  variation  of  the  standard  quartz  microprobe  design, 
referred  to  as  an  Electro-Mechanical  Sonic  (EMS)  probe  developed  by  Puri  [36],  is  shown  in  Figure 
1  and  enables  sampling  in  flame  regions  where  soot  volume  fractions  are  on  the  order  of  10*5. 

The  EMS  probe  utilizes  a  fiber  that  extends  through  the  orifice  and  is  attached  to  a  solenoid  core 
with  springs  on  either  side  of  the  core.  A  solenoid  is  placed  around  the  outside  of  the  quartz  tube  and 
the  current  to  the  solenoid  is  interrupted  at  20  -  30  Hz,  providing  for  constant  oscillatory  motion  of 
the  fiber  within  the  orifice,  thus  keeping  an  annular  region  open  for  gas  sampling.  The  fiber  for 
this  study  is  =125  pm  in  diameter  and  is  made  from  sapphire,  while  the  orifice  diameter  is  =200 
pm.  The  annular  sampling  region  has  an  area  equivalent  to  a  =150  pm  diameter  orifice. 

At  the  soot  levels  encountered  in  this  flame  (volume  fraction  up  to  10"5),  it  is  required  to  supplement 
the  EMS  probe  with  a  rapid  insertion  procedure  similar  to  that  originally  described  by  Kent  and 
Wagner  [19]  for  use  with  thermocouples  in  flames  with  high  soot  concentrations.  This  procedure 
involves  alternately  removing  the  probe  tip  from  the  sample  location  in  the  flame  to  a  lean,  high 
temperature  region  where  soot  will  burn  from  the  tip  leaving  a  clean  orifice  before  moving  to  the 
next  sample  location. 

SPECIES  ANALYSIS 

Analysis  of  the  collected  gas  samples  was  performed  using  a  gas  chromatograph  (GC)  and  a  mass 
spectrometer  (MS).  The  use  of  these  two  techniques  allowed  an  independent  verification  of  the 
concentration  profile  measurements  as  well  as  providing  specific  species  measurement  capability 
as  described  below. 

Gas  Chromatograph:  A  Hewlett-Packard  5890  GC  employed  Porapak  Q  and  Porapak  S  columns 
coupled  with  a  molecular  sieve  (5A)  for  species  separation  using  a  column  sequence  reversal 
technique  and  temperature  programming.  A  dry  ice  bath  was  used  to  extract  H20  from  the  sample 
to  simplify  storage  requirements.  A  Thermal  Conductivity  Detector  (TCD)  was  used  for  detection 
of  all  non-hydrocarbons  (N2,  02/Ar,  CO,  and  C02).  A  Flame  Ionization  Detector  (FID)  was  used 
for  detection  of  hydrocarbons  in  small  concentrations.  For  concentrations  of  hydrocarbons  which 
were  beyond  the  linear  range  of  the  FID,  the  TCD  was  used.  A  signal  integrator  is  attached  to  each 
detector  for  the  recording  of  the  chromatogram. 

Since  the  gas  chromatograph  requires  the  sample  to  be  at  atmospheric  pressure  for  analysis  and  the 
EMS  probe  collects  sample  gas  at  low  pressure,  a  method  was  devised  to  collect  a  sample  in  a 
cylinder  and  compress  it  with  a  piston  for  storage  in  a  multiloop  sample  storage  system  [36].  The 
schematic  of  this  compression  apparatus  is  shown  in  Figure  2.  Each  sample  is  stored  in  a  10  cm3 
loop  at  a  pressure  of  two  atmospheres.  After  fifteen  samples  are  collected  in  separate  storage  loops, 
the  sample  storage  system  is  connected  to  the  GC  and  each  sample  is  analyzed  sequentially.  One  of 
the  GC  signal  integrators  controls  the  GC  and  the  multiposition  valve  connected  to  the  sample 
storage  system. 
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Calibration  of  the  GC  was  primarily  done  using  commercially  available  gas  mixtures  with 
accuracy  of  ±2%.  The  GC  was  also  calibrated  for  fuel  species  concentrations  much  greater  than 
that  available  commercially  using  laboratory  prepared  mixtures.  The  GC  signal  integrators  used 
with  the  TCD  and  FID  were  able  to  perform  multi-level  calibrations  for  species  which  spanned 
large  ranges  in  concentrations.  The  molecular  species  calibrated  and  measured  directly  with  the 
GC  are  No,  02/Ar,  CO2,  CO,  CH4,  C2H2,  C2H4,  C2H6,  CgHg,  C3H8,  C4H6  and  C6H6.  The  columns 
cannot  separate  02  and  Ar.  Therefore,  a  procedure  which  assumes  the  mole  fractions  of  N2  and  Ar 
remains  in  the  same  ratio  as  in  air  is  used  to  determine  the  mole  fraction  of  02  and  Ar 
individually. 

For  optimal  detection  of  species  with  the  TCD,  helium  was  used  as  the  carrier  gas.  Since  detection 
depends  on  the  difference  in  thermal  conductivities,  measurements  of  molecular  hydrogen  are 
problematic.  In  addition,  the  samples  were  dried  prior  to  storage,  thus,  the  H20  concentration  is 
obtained  by  a  hydrogen  atom  balance  approach  developed  by  Puri  [36]  which  assumes  all  species 
diffuse  similarly.  This  method  calculates  the  total  amount  of  H  atom  present  in  the  sample  from 
the  hydrocarbon  species  present,  equation  la.  The  total  amount  of  carbon  atom  is  then  calculated 
from  the  carbon  containing  species,  equation  lb.  Then  the  total  amount  of  H  atom  which  should  be 
present  in  the  sample  based  on  the  amount  of  carbon  present  and  the  carbon  to  hydrogen  ratio  in  the 
fuel,  and  the  actual  H  atom  present  are  used  to  calculate  the  amount  of  H20  present  in  the  flame  at 
that  position,  equation  lc.  Finally  the  dry  measurements  were  normalized  with  the  calculated  H20 
concentration  assuming  the  measured  species  mole  fractions  sum  to  one.  This  method 
overestimates  the  concentration  of  H20  since  it  ignores  molecular  hydrogen,  but  the  mole 
fractions  are  unaffected  since  both  species  contain  two  hydrogen  atoms. 


«  =  S-vXc,», 

(la) 

C  ~  ^  HvO. 

(lb) 

X„,o=<7c„yC-H)l2 

(10 

Since  not  all  hydrocarbons  in  the  flame  sample  are  quantitatively  determined  and  the  hydrogen  to 
carbon  ratio  in  these  unsaturated  hydrocarbons  and  soot  is  less  than  in  the  fuel,  the  amount  of 
H20/H2  is  underpredicted  as  a  result.  This  factor  along  with  the  higher  diffusivity  of  H2  could  be  a 
factor  in  bringing  the  H20  measurement  back  in  line  with  what  is  actually  present  in  the  flame. 

Mass  Spectrometer:  The  mass  spectrometer  employed  an  Extrel  model  7-324-9  quadrupole  mass 
filter  with  mass  range  of  1-500  u  for  separation  of  ions  and  an  electron  multiplier  for  detection.  A 
computer  controlled  the  ionizer,  the  mass  filter  and  the  burner  position  while  the  signal  was 
recorded  by  the  computer  through  an  A/D  converter. 

As  shown  in  Figure  3,  a  staged  vacuum  pumping  system  was  employed.  The  EMS  probe  was  used  to 
continuously  extract  sample  gas  from  the  flame.  The  pressure  in  the  small  vacuum  chamber  was 
on  the  order  of  1  Torr.  The  orifice  of  a  second  quartz  probe  was  situated  in  the  exit  stream  of  the 
EMS  probe  and  allowed  a  small  amount  of  the  flame  sample  gas  to  pass  into  the  high  vacuum 
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chamber  where  the  pressure  was  on  the  order  of  10'6  Torr.  The  secondary  probe  discharges  the 
sample  into  the  electron  impact  ionizer  which  has  tungsten  filaments.  The  resultant  ions  were 
extracted,  focused  and  accelerated  into  the  mass  filter  for  separation.  The  mass  filter  was  operated 
in  single  ion  mode  and  the  ionizer  was  operated  at  20  eV  or  less  to  reduce  molecular 
fragmentation. 

The  calibration  method  for  the  MS  follows  the  procedure  given  by  Ermolin  et  al.  [37]  as  shown  in 
equation  2  where  Ft  and  St  are  the  calibration  factor  and  signal  intensity  respectively  for  species  i. 
This  method  is  mathematically  similar  to  Bittner  [31]  and  Crowhurst  and  Simmons  [33].  The 
method  relies  on  the  signal  response  relative  to  other  species  remaining  constant.  That  is,  for  a 
given  mixture  the  ratio  of  the  signals  due  to  two  gases  will  always  be  the  same  and,  therefore,  is 
independent  of  variations  in  mass  flow  through  the  orifice  due  to  orifice  clogging  or  variations  in 
density.  Secondly,  this  method  assumes  that  the  species  measured  constitute  all  of  the  gases 
present  which  are  normalized  to  a  mole  fraction  of  one. 


X..= 


F;S; 


If ,s, 


(2) 


Determination  of  the  calibration  factors  were  performed  with  a  gas  of  known  mole  fraction  of  a 
given  species  flowing  through  the  burner  fuel  tube  so  that  any  biases  relating  to  the  sampling 
procedure  or  MS  characteristics  are  included  in  the  calibration  factor,  Fr  Calibrations  for  N2,  02, 
and  Ar  used  air  as  the  standard  gas.  Calibration  of  H20  employed  a  N2  stream  saturated  with  H20 
vapor  at  a  known  temperature.  Other  species  used  commercial  calibration  mixtures  or  calibrated 
flowmeters  to  mix  pure  gases  in  the  appropriate  ratios.  The  molecular  species  quantitatively 
measured  by  the  MS  are  N2/CO,  02,  H20,  C02,  H2,  CH4,  C2H2,  C2H4,  and  Ar. 


One  limitation  of  the  use  of  a  MS  in  the  present  study  is  the  coincidence  of  C2H4,  N2  and  CO  at  mass 
28.  For  C2H4,  the  fragment  at  mass  27  was  used  as  an  independent  measurement.  Calibration  of 
C2H4  at  both  masses  27  and  28  was  necessary  so  that  the  signal  contribution  of  C2H4  at  mass  28 
could  be  subtracted  from  the  mass  28  signal.  C2H4  was  similarly  calibrated  and  subtracted  at  mass 
26  where  C2H2  was  measured. 

Furthermore,  it  is  not  possible  to  determine  the  signal  of  N2  or  CO  independently  with  the  MS. 
Therefore,  the  analysis  is  performed  as  though  there  was  no  CO  present.  Since  CO,  as  determined 
by  the  GC  measurements,  can  yield  mole  fractions  of  close  to  0.10  at  some  locations  in  the  flame, 
this  approach  will  overpredict  the  N2  by  somewhat  more  than  the  measured  CO  concentration  due  to 
the  greater  sensitivity  of  the  MS  to  CO  and  reduce  the  mole  fractions  of  other  species  slightly  due  to 
the  normalization  procedure  shown  in  equation  2. 

To  determine  the  centerline  position,  the  MS  data  was  used  since  these  data  were  obtained  across 
the  entire  diameter  of  the  flame  and  were  generally  symmetrical.  The  primary  method  was  to 
assign  the  two  locations  where  fuel  and  oxygen  disappear  to  the  same  radius.  At  higher  axial 
positions  where  the  fuel  was  totally  consumed,  either  the  C2H2  or  H2  profiles  were  used  in  a  similar 
manner  or  directly  to  center  the  data.  Due  to  sample  storage  limitations,  samples  for  GC  analysis 
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were  obtained  for  only  one-half  of  the  flame  diameter.  Therefore,  the  radial  position  of  the  GC  data 
was  determined  by  setting  the  position  where  fuel  and  oxygen  disappear  to  match  the  MS 
measurements. 

Measurements  using  the  above  instrumentation  will  be  combined  with  the  temperature  (fine  wire 
thermocouples),  velocity  (Laser  Doppler  Anemometer)  and  soot  aerosol  information  (laser 
extinction  and  scattering)  which  have  been  published  earlier  for  this  same  flame  [12-14],  The 
experimental  details  of  those  measurements  are  not  discussed  here. 

RESULTS 


Spatially  resolved  measurements  of  molecular  species  have  been  determined  throughout  the 
region  of  the  flame  where  soot  undergoes  growth  processes.  Measurements  at  a  given  axial 
position  were  taken  at  0.5  mm  radial  steps  with  the  MS,  while  the  GC  measurements  were  limited 
to  15  samples  per  experiment,  so  samples  were  taken  judiciously  to  resolve  the  important  features 
in  the  flame.  For  the  species  measured  by  both  instruments  the  results  generally  agreed  to  ±20%  or 
better  with  the  exception  of  acetylene  which  displayed  a  somewhat  larger  difference  for  the  two 
measurement  techniques.  However,  due  to  the  subtraction  procedure  for  calculating  the  mole 
fraction  of  C2H2  with  the  MS  the  measurements  of  acetylene  agreed  within  the  measurement 
errors.  As  an  illustration  of  the  spatial  distribution  of  the  mole  fractions  of  C2H4  and  C2H2,  Figure 
4  shows  the  radial  measurements  of  these  species  for  several  heights  in  the  flame  for  both  the  GC 
and  MS.  Note  that  the  rapid  consumption  of  the  fuel  results  in  a  maximum  C2H4  mole  fraction  of 
less  than  0.01  at  30  mm  above  the  burner  outlet  while  the  luminous  flame  height  is  approximately 
88  mm. 

Examination  of  the  radial  profile  for  the  C2H2  mole  fraction  at  the  7  and  10  mm  heights  indicates 
the  maximum  value  occurs  at  a  radial  position  different  from  the  centerline.  This  result  implies 
that  the  highest  production  rates  for  C2H2  occurs  in  an  annular  region  approximately  4-5  mm  in 
radius.  The  temperature  and  C2H4  mole  fraction  in  this  region  of  the  flame  ranges  from 
approximately  1300  -  1700  K  and  0.0345  to  0.185  respectively.  The  quantity  of  C2H2  at  the  centerline 
indicates  some  production  of  C2H2  throughout  the  central  region  of  the  flame  where  the  C2H4  mole 
fractions  are  much  greater.  In  fact,  the  increase  in  the  maximum  mole  fraction,  indicating  net 
production  of  C2H2,  continues  through  to  the  20  mm  axial  measurement  location  where  the 
temperature  is  1421  K.  After  the  30  mm  height,  the  peak  C2H2  mole  fraction  decreases  by  a  factor  of 
about  four  by  40  mm  to  0.01  and  further  decreases  by  a  factor  of  about  ten  by  50  mm  to  0.001  as 
reactions  proceed.  The  peak  in  the  integrated  soot  volume  fraction  in  this  flame  occurs  at  the  40 
mm  height  and  the  concentrations  measured  show  only  a  small  amount  of  measured 
hydrocarbons  remaining  to  be  oxidized.  At  this  height,  soot  represents  14%  of  the  total  carbon  flow 
rate  in  the  flame  [14].  The  major  gaseous  species  remaining  to  be  oxidized  at  this  height  are  CO, 
H2,  and  C2H2  with  peak  mole  fractions  of  approximately  0.07.  0.02,  and  0.01  respectively. 

Several  other  intermediate  hydrocarbon  species  were  measured  quantitatively  by  the  GC  and  may 
be  considered  important  in  soot  formation  in  this  flame.  These  include  C4H6  and  C6H6  shown  as 
radial  profiles  at  several  axial  positions  in  Figure  5.  C4H6  is  closely  related  to  C4H5  which  has 
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been  proposed  to  play  an  important  role  in  the  formation  of  the  first  aromatic  ring  [6]  and  C6H6  is 
the  first  aromatic.  Note  that  the  mole  fractions  for  these  species  have  a  spatial  development  which 
is  characteristic  of  intermediate  species  in  this  flame,  that  is  similar  to  C2H2.  For  these 
intermediate  species,  the  mole  fraction  is  low  initially,  sometimes  with  a  peak  concentration  in 
the  annular  region  or  a  broad  region  where  the  mole  fraction  is  relatively  constant.  The 
concentration  increases  with  height  and  the  maximum  moves  towards  the  centerline,  and  then  the 
concentration  decreases  as  consumption  becomes  dominant.  This  behavior  has  been  seen 
previously  in  similar  flames  burning  methane  [38].  The  species  CH4,  C2H6  and  C3H6  also  show 
similar  characteristics  during  their  development  through  the  flame.  These  species  have 
maximum  mole  fractions  in  the  flame  of  0.0038,  0.00027  and  0.00020  respectively. 

Finally,  several  peaks  on  the  GC  chromatogram  were  present  throughout  much  of  the  fuel  rich 
region.  However,  only  the  elution  times  and  integrated  areas  are  available  since  these  species  are 
not  identified.  A  double  peak  elutes  very  close  in  time  to  that  of  C3H6  and  C3H8.  Analysis  of  a 
flame  sample  on  a  GC  with  an  HP-1  capillary  column  shows  a  large  peak  with  a  retention  index 
very  close  to  both  propadiene  (allene)  and  propyne  (methyl- acetylene),  both  are  C3H4  structures. 
Propadiene  is  related  to  the  propargyl  radical  (C3H3)  which  has  been  proposed  to  play  an  important 
role  in  the  formation  of  the  first  aromatic  ring  directly  from  C3  species  [23].  Estimating  the 
calibration  factor  for  these  unidentified  species  as  similar  to  that  of  the  C3H6  and  C3H8  would  give 
these  species  maximum  mole  fractions  greater  than  100  ppm. 

Measurements  were  also  taken  with  the  MS  at  masses  50  and  78.  The  molecular  formula  for  these 
masses  are  C4H2  and  C6H6  respectively.  These  species  have  no  calibrations  for  them  so  only 
relative  signal  values  are  available.  These  masses  could  also  be  fragments  of  species  at  larger 
masses.  For  instance,  at  mass  50,  the  signal  could  be  a  combination  of  this  species  that  is  present 
in  the  flame  and  other  species  from  the  flame  with  larger  masses  (e.g.  C4H4,  C4H6)  which  have 
been  fragmented  in  the  ionizer  of  the  MS.  The  measurements  of  these  intermediate  species  with  the 
MS  indicate  a  similar  development  through  the  flame  as  found  in  the  GC  measurements  for  other 
intermediate  species. 

Most  of  the  discussion  which  follows  will  focus  on  the  measurements  from  the  GC.  This  is  due  to 
the  lack  of  interference  between  species  at  different  masses  from  fragmentation  during  ionization 
in  the  MS.  Moreover,  identification  and  quantification  of  several  additional  species  is  possible 
with  the  GC. 

DISCUSSION 


One  question  which  arises  in  both  premixed  and  diffusion  flames  with  respect  to  soot  formation 
pertains  to  the  cessation  of  growth  of  soot  particles  in  the  flame.  In  rich  premixed  flames,  it  has 
been  found  that  significant  quantities  of  C2H2  and  other  hydrocarbons  are  present  in  the  post  flame 
region  [5,  8] .  It  is  in  this  region  that  particle  inception  occurs  and  the  surface  growth  reactions 
proceed.  Furthermore,  the  rate  of  these  surface  reactions  decreases  with  time  while  the  smaller 
hydrocarbon  species  concentrations  do  not  change  appreciably.  This  experimental  result  has  been 
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used  to  argue  that  the  reactivity  of  the  particles  decreases  with  time  and  has  been  discussed  by 
others  [15,  39.  40]. 


In  diffusion  flames  with  significant  soot  levels,  no  experimental  information  has  been  available 
about  the  chemical  environment  surrounding  the  soot  particles  as  they  traverse  the  flame.  The 
combination  of  the  previous  measurements  of  Santoro  and  coworkers  [12-14]  with  the  spatial 
species  measurements  enables  analysis  which  can  help  to  elucidate  greater  understanding  about 
some  of  the  soot  chemistry  which  is  occurring. 

The  following  analysis  integrates  the  previous  measurements  of  flowfield  and  soot  with  the 
molecular  species  information  now  available.  From  the  velocity  and  temperature  measurements, 
soot  particle  paths  have  been  determined  and,  thus,  the  spatial  measurements  can  be  transformed 
into  temporal  measurements  along  a  particle  path  [14].  Consequently,  the  detailed  temporal 
evolution  of  soot  particle  growth  with  respect  to  its  thermochemical  environment  can  be  obtained. 

In  the  present  case,  the  soot  growth  process  will  be  examined  in  relation  to  the  concentration  field  of 
several  hydrocarbon  species.  In  particular,  the  C2H2  mole  fraction  will  be  examined  in  greater 
detail  since  this  species  has  been  argued  to  be  the  primary  growth  species  in  premixed  flames  [8] . 

In  examining  this  process,  a  total  of  seven  particle  paths  will  be  considered  which  are  numbered  1  - 
7  with  the  centerline  being  particle  path  1  and  proceeding  radially  outwards  through  particle  path 
7 .  Three  particle  paths  will  be  discussed  in  greater  detail.  One  particle  path  is  the  centerline  of  the 
flame  (PP  1).  A  second  considers  the  particle  path  which  passes  through  the  maximum  observed 
local  soot  volume  fraction  for  the  entire  flame  (PP  6).  A  third  particle  path  is  partway  between  these 
two  and  is  designated  the  middle  particle  path  (PP  3).  For  the  paths  away  from  the  centerline,  the 
appropriate  data  from  both  sides  of  the  flame  have  been  averaged  and  then  the  analysis  performed. 
This  is  due  to  the  symmetric  nature  of  the  flame  and  the  arbitrary  choice  which  would  be  necessary 
when  combining  the  two  sets  of  data. 

Figure  6  shows  the  local  soot  volume  fraction,  fv,  and  the  C2H2  mole  fraction  as  a  function  of  time 
along  particle  paths  1,  3  and  6.  The  trends  on  these  three  particle  paths  are  typical  of  the  evolution 
observed  along  other  particle  paths  in  the  flame.  The  important  characteristic  is  that  the  soot 
volume  fraction  achieves  a  maximum  at  the  same  axial  location  where  the  C2H2  concentration 
falls  to  near  or  below  detectability  (approximately  2xl0'5  for  the  GC).  The  measurements  of  the 
other  hydrocarbons  from  both  the  GC  and  MS  indicate  a  similar  depletion  at  the  same  height  or 
even  earlier  in  the  flame.  Thus,  it  is  argued  that  the  depletion  of  hydrocarbons  available  for  soot 
particle  growth  is  the  reason  for  particle  growth  ceasing.  The  present  measurements  for  this 
ethene/air  diffusion  flame  show  a  different  mechanism  for  the  cessation  of  soot  growth  is 
operative  as  compared  to  pre mixed  flames. 

The  underlying  phenomenon  of  this  growth  cessation  presented  above  has  two  parts.  These  are  best 
explained  by  looking  at  a  general  rate  equation  for  soot  formation  in  equation  3.  It  is  first 
important  to  notice  that  when  the  changes  in  the  local  soot  volume  fraction  are  plotted  on  a  time 
coordinate  along  a  particle  path,  it  is  the  net  result  of  the  growth  and  oxidation  rates  which  is 
determined.  So  obviously,  a  transition  from  net  soot  formation  to  net  soot  destruction  occurs  when 
the  soot  concentration  reaches  its  maximum  along  a  given  particle  path. 
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d)xoot  =  Growth  -  Oxidation 


(3) 


When  one  proceeds  to  examine  why  such  a  transition  occurs,  it  is  necessary  to  consider  the 
structure  of  the  flame  and  the  specifics  of  the  soot  growth  process.  In  a  coannular,  laminar 
diffusion  flame,  fuel  is  transformed  into  products  as  one  proceeds  higher  in  the  flame.  The  data 
from  the  present  study  shows  that  beyond  the  30  mm  height  above  the  burner  exit,  all  measured 
hydrocarbon  concentrations  decrease  dramatically.  This  decrease  in  hydrocarbons  available  for 
soot  growth  will  slow  the  growth  reactions  since  this  rate  would  generally  be  dependent  upon  the 
concentration  of  the  hydrocarbon  species  participating  in  surface  growth  reactions. 

Secondly,  this  decrease  in  hydrocarbons  is  followed  by  an  increase  in  species  (especially  02) 
available  to  participate  in  the  oxidative  reactions.  The  concentrations  of  H20  and  C02  which  could 
also  participate  in  soot  oxidation  is  relatively  high  everywhere  in  the  flame.  Coincident  with  this 
transition  from  hydrocarbon  rich  to  oxygen  rich  is  a  reaction  zone  with  an  (associated) 
temperature  rise.  Therefore,  the  growth  reactions  in  equation  3  decrease  and  the  oxidation 
reactions  increase,  which  yields  the  position  of  maximum  local  soot  volume  fraction  along  a 
given  particle  path.  This  corresponds  to  the  position  where  the  local  equivalence  ratio  based  on 
measured  gaseous  species  as  calculated  with  the  method  discussed  by  Mitchell  et  al.  [41]  has  fallen 
below  a  value  of  one  as  shown  in  figure  7. 

Additional  evidence  that  supports  a  mechanism  of  soot  growth  cessation  in  laminar  diffusion 
flames  which  differs  from  that  dominant  in  premixed  flames  can  be  gleaned  by  examining  the 
temperatures  along  these  particle  paths  and  using  the  analysis  results  from  premixed  flames  of 
Dasch  [39]  and  Woods  and  Haynes  [40]  to  examine  the  decrease  in  soot  particle  reactivity. 
Equation  4  expresses  the  decay  in  the  reactivity  of  the  soot  surface  as  taken  from  Dasch  [39]  and 
Woods  and  Haynes  [40].  The  term  kSG ,  which  was  determined  from  Woods  and  Haynes  [40],  is  the 
same  term  as  a  used  by  Dasch  [39].  Figure  8  shows  the  temperatures  along  the  three  particle  paths 
with  the  time  a  particle  path  crosses  the  3  mm  axial  height  set  to  zero. 


exp(-ar) 


(4a) 


a  -  3. 5 X107  exp(-7^  /  T) 


(4b) 


While  it  may  not  be  good  practice  to  apply  the  above  formula  for  reactivity  decay  to  non-isothermal 
flames  [15,  42]  there  is  insight  to  be  gained  by  proceeding  along  this  line  of  analysis.  If  the 
maximum  temperature  along  the  particle  path  for  the  soot  growth  region  (i.e.  from  initial  detection 
of  soot  to  the  maximum  soot  volume  fraction  along  a  given  particle  path)  is  chosen,  one  can 
calculate  a  conservative  estimate  of  the  decrease  in  the  soot  particle  reactivity  along  that  particle 
path  based  on  the  results  from  premixed  flames.  For  particle  paths  1,  3  and  6  the  maximum 
temperatures  and  soot  residence  times  are  shown  in  Table  1.  Data  for  two  additional  flames  are 
also  shown  in  Table  1.  The  first  is  an  atmospheric  pressure  premixed  flame  burning  C2H4/air  and 
the  second  is  the  centerline  of  a  C2H4/air  diffusion  flame  discussed  by  Honnery  and  Kent  [43]. 

The  latter  flame  has  a  residence  time  that  is  about  five  times  longer  than  any  observed  for  the 
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TABLE  1  Data  for  and  calculation  of  particle  reactivity  decrease  along 
several  particle  paths 


Particle 

Path 

Maximum 
Temperature  (K) 

Soot  Growth 
Residence 

Time  (s) 

Characteristic 
Time  (s) 

Fraction 
of  Initial 
Reactivity 

PP  1 

1439 

0.017 

0.196 

0.92 

PP  3 

1567 

0.035 

0.054 

0.69 

PP  6 

1613 

0.026 

0.036 

0.49 

Harris  and 

1750 

0.02 

0.012 

0.19 

Weiner  [9] 

Honnery 

=1250 

>0.2 

2.11 

0.91 

and  Kent 

[43] 

present  study.  The  characteristic  time,  I/a,  for  soot  particle  reactivity  decrease  is  calculated  using 
the  activation  temperature  and  expression  of  Dasch  [39]  (see  equation  3)  which  uses  an  activation 
temperature  of  22650  K.  The  fraction  of  initial  reactivity  that  would  remain  when  growth  ceases  in 
this  flame  is  also  shown  in  Table  1.  The  Arrhenius  nature  of  the  expression  for  the  reactivity 
decrease  can  be  expected  to  differ  somewhat  in  the  temperature  range  observed  in  the  present 
diffusion  flame.  However,  Woods  and  Haynes  [40]  did  not  expect  this  effect  to  be  significant. 

Particle  path  7  (not  shown)  has  the  highest  maximum  temperature  (1786  K)  and  the  fraction  of 
initial  reactivity  by  this  conservative  estimate  is  0.03.  However,  the  temperature  during  much  of 
the  growth  along  this  and  all  the  other  particle  paths  is  typically  1300  -  1600  K.  So  the  values  of  the 
fraction  of  initial  reactivity  from  Table  1  are  more  typical  of  that  in  the  bulk  of  the  flame. 
Furthermore,  when  the  data  from  Honnery  and  Kent  [43]  is  used  along  the  centerline  of  their 
flame,  this  analysis  yields  a  fraction  of  initial  reactivity  available  at  the  end  of  growth  of  0.91. 
These  values  for  the  surface  reactivity  indicate  that  the  continued  growth  of  soot  in  diffusion 
flames  to  very  long  residence  times  is  not  in  conflict  with  the  results  obtained  by  others  in 
premixed  flames.  The  premixed  flame  has  higher  temperatures  and  comparable  residence  time  to 
the  current  flame,  resulting  in  a  lower  reactivity.  These  calculations  further  support  the 
interpretation  that  soot  surface  growth  in  laminar  diffusion  flames  ceases  due  to  the  depletion  of 
gas  phase  reactants. 

Rate  Analysis :  Previous  analysis  of  soot  surface  growth  in  laminar  diffusion  flames  has  been 
limited  to  analysis  of  data  from  the  soot  aerosol  itself  [14,  25,  44]  due  to  a  lack  of  species 
measurements  in  flames  with  significant  soot  levels.  To  perform  an  analysis  of  the  chemistry  of 
the  growth  process  requires  the  determination  of  a  specific  chemical  mechanism.  Surface  growth 
of  soot  particles  can  result  from  a  number  of  processes  including  heterogeneous  reactions  of  a 
gaseous  molecule  with  the  surface  and  condensation  of  large  hydrocarbons,  generally  PAH,  on  the 
surface.  The  mechanism  which  has  received  much  attention  is  a  reaction  between  the  soot  surface 
and  a  gas  phase  hydrocarbon.  The  work  of  Harris  and  Weiner  [8]  has  been  instrumental  in 
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determining  constraints  on  the  potential  growth  species  which  could  be  important  in  the  surface 
growth  process  in  premixed  flames.  Their  conclusion  was  that  the  rate  of  growth  of  the  volume 
fraction  of  soot  is  a  result  of  the  heterogeneous  reaction  of  C2H2  with  the  soot  surface  and  can  be 
expressed  as  equation  5.  With  the  data  available,  the  rate  coefficient  in  this  expression  can  be 
determined  and  compared  to  premixed  flames. 

^  soot  ~  kc2  h2  SXCi  n2  ( 5 ) 

This  equation  assumes  that  the  oxidation  rates  under  the  fuel  rich,  premixed  flame  conditions  they 
studied  are  much  smaller  than  the  growth  rates,  and  the  species  terms  in  the  equation  describing 
general  soot  growth,  equation  3,  are  small  except  for  the  term  involving  C2H2.  Harris  and  Weiner 
[10]  estimated  that  only  10-15%  of  the  growth  of  soot  comes  from  C4H2.  It  was  also  found  that  C2H2 
explained  the  soot  growth  rate  in  an  ethene/toluene  premixed  flame  [10],  where  levels  of  aromatic 
species  are  expected  to  be  much  greater  than  in  a  pure  C2H4  flame.  Moreover,  recent  modeling 
work  [22]  has  shown  that  C2H2  can  explain  the  growth  of  soot  in  high  pressure,  premixed  ethene 
flames. 

In  previous  studies  of  both  premixed  and  diffusion  flames  [8,  17,  25,  29,  44,  45],  the  formation  rate  of 
soot,  CQsn()[ ,  along  a  particle  path  was  calculated  from  equation  6.  This  expression  is  valid  when  the 
gas  density  (temperature  and  molecular  weight)  remains  constant.  However,  because  the  volume 
fraction  is  tied  directly  to  the  gas  density,  the  mass  growth  of  soot  becomes  equation  7  when  the  soot 
conservation  equation  is  solved.  This  expression  assumes  an  ideal  gas  and  a  constant  molecular 
weight.  The  latter  assumption  is  well  satisfied  in  this  flame  since  the  fuel  molecular  weight  is 
28.05  and  the  air  has  an  average  molecular  weight  of  28.96.  Measurements  of  the  average 
molecular  weight  also  bear  this  out.  As  can  be  seen  from  figure  8  the  rate  of  change  of  temperature 
is  relatively  low  except  as  the  particle  path  crosses  a  reaction  zone.  In  the  case  of  a  particle  path 
exiting  the  fuel  rich  region,  this  increased  density  correction  term  coincides  with  the  soot  volume 
fraction  reaching  its  peak  and  the  corresponding  mass  growth  rate  of  soot  being  small.  Therefore, 
measurements  in  this  region  can  have  large  relative  corrections  due  to  changes  in  gas  density. 

For  this  reason,  locations  where  this  correction  is  greater  than  20%  of  df\/  ,  are  not  included  in  the 

i  *  /dt 

analysis. 


<»=Psoot 


dt 


(6) 


®  =  P« 


(4L+L*L) 

k  dt  T  dt  ) 


(7) 


As  previously  mentioned,  the  changes  in  the  soot  mass  along  a  particle  path  is  a  reflection  of  the 
difference  of  the  growth  rate  and  the  oxidation  rate  of  soot,  equation  3.  Therefore,  it  could  be 
important  to  account  for  oxidation  reactions  when  determining  the  rate  of  reactions  contributing  to 
the  growth  of  soot.  The  major  species  available  to  participate  in  soot  oxidation  are  02,  CO2  and 
H2O.  Rate  expressions  for  the  reaction  rates  of  these  species  with  soot  are  given  in  references  [46- 
48],  As  seen  in  figure  9  there  is  generally  about  0.003  mole  fraction  of  02  in  the  fuel  rich  region  of 
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the  flame.  This  has  been  noted  previously  in  coannular  diffusion  flames  [34,  38.  41,  49]  and  is 
generally  attributed  to  leakage  of  02  at  the  burner  lip  where  reactions  are  quenched  [41].  In 
addition,  the  mole  fractions  of  C02  and  H20  in  the  soot  growth  region  range  from  0.06  to  0.14  and 
0.14  to  0.19  respectively. 

The  rate  expression  for  02  from  Bradley,  Dixon-Lewis  et  al.  [48]  is  the  mechanism  from  Nagle 
and  Strickland-Constable  [50] .  This  mechanism  for  soot  oxidation  has  been  used  quite  extensively 
in  recent  studies  [26,  51,  52].  This  expression  was  validated  by  Park  and  Appleton  [53]  in  a  shock 
tube  study  of  soot  oxidation.  However,  the  range  of  temperature  and  02  concentrations  in  that  study 
are  above  those  found  in  this  experiment.  In  addition,  the  Nagle  and  Strickland-Constable  rate 
has  been  found  to  overpredict  the  oxidation  rate  of  soot  at  lower  temperatures  [54,  55],  In 
comparison,  the  02  reaction  mechanism  from  Libby  and  Blake  [46,  47]  indicates  a  reaction  rate 
from  five  to  ten  times  faster  than  the  Nagle  and  Strickland-Constable  rate  [50]  at  the  temperatures 
and  02  mole  fractions  encountered  in  this  flame.  Additional  evidence  that  the  Nagle  and 
Strickland-Constable  [50]  rate  could  be  too  fast  is  the  calculated  oxidation  rate  in  the  lean  part  of 
several  C2H2  flames  exceeding  the  measured  rate  of  soot  oxidation  by  a  factor  of  roughly  7  [52]. 

Further  examination  of  the  rate  for  C02  oxidation  by  the  Libby  and  Blake  [46,  47]  mechanism  is 
generally  within  a  factor  of  two  of  the  rate  given  by  the  Nagle  and  Strickland-Constable  [50] 
mechanism  for  02  at  their  respective  concentrations.  Therefore,  these  oxidation  rates  will  not  be 
used  in  further  analysis.  Though  the  evidence  cited  above  indicates  that  the  Nagle  and 
Strickland-Constable  [50]  soot  oxidation  rate  is  not  accurate  under  the  conditions  present  in  this 
flame,  this  is  still  the  most  popular  mechanism  for  calculating  soot  oxidation  rates  by  02  by  those 
mentioned  earlier.  Therefore,  these  02  rates  and  the  rates  for  C02  and  H20  from  Bradley  et  al.  [48] 
will  be  used.  In  the  region  of  large  net  soot  growth,  the  correction  due  to  oxidation  is  on  the  order  of 
10%  of  >  primarily  from  02.  However,  near  the  end  of  soot  growth,  this  correction  becomes  a 
much  greater  fraction  of  the  growth  rate.  The  following  analysis  will  exclude  locations  where  the 
oxidation  correction  is  greater  than  50%  of  the  df  ./  . 

Two  procedures  were  performed  to  find  4/^/  in  this  two  dimensional  flame.  For  the  first 
procedure,  the  first  step  is  to  transform  the  volume  fraction  data  along  a  given  particle  path  into 
rates  of  growth  )  at  the  measurement  points.  To  accomplish  this,  a  weighted  polynomial  curve 
fit  (typically  fifth  order)  was  performed  on  the  volume  fraction  data  and  the  derivative  taken.  This 
approach  was  chosen  since  all  of  the  data  (fv,  species,  etc.)  were  taken  at  the  same  axial  positions. 
The  soot  volume  fraction  data  and  the  particle  paths  coordinates  for  this  analysis  were  taken  from 
Santoro  et  al.  [14].  Deciding  the  best  fit  to  represent  the  soot  growth  when  varying  the  polynomial 
order  and  capturing  the  transition  from  growth  to  oxidation  was  somewhat  subjective. 

Therefore,  a  second  approach  which  used  data  available  in  our  lab  through  the  technique  of  2-D 
Laser  Induced  Incandescence  (LII)  was  performed  [56] .  The  steps  were  to  take  the  2-D  image  of  the 
soot  volume  fraction  field  and  extract  a  smaller  section  with  dimensions  6.02  mm  axially  and 
1.35  mm  radially  (this  corresponds  to  a  section  41  by  10  pixels  respectively)  that  was  centered  on  the 
location  of  interest.  A  third  order,  polynomial  least  squares  fit  of  this  surface  (ten  coefficients) 
was  performed  on  this  smaller  section  of  volume  fraction  data.  The  particle  path  locations  and  the 
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velocities  were  taken  from  [14].  The  measured  velocities  are  important  for  the  directional 
derivative  along  the  surface,  equation  8. 

The  two  methods  of  determining  L^/dt  give  comparable  growth  rates.  Typically,  the  LII  method 
resulted  in  slightly  greater  values.  On  average,  the  LII  results  were  approximately  30%  greater 
than  that  determined  by  curve  fitting  and  most  were  within  a  factor  of  two.  Many  of  the  instances 
where  the  two  growth  rates  differed  by  the  largest  amounts  were  at  the  point  where  soot  is  first 
detected  or  near  the  peak  in  soot  volume  fraction  where  the  transition  from  growth  to  oxidation 
occurs.  Both  of  these  locations  are  difficult  to  capture  accurately  with  the  limited  data  available 
along  a  particle  path  from  light  scattering  and  extinction  since  data  is  only  taken  every  10  mm  in 
the  axial  direction.  Other  factors  also  contribute  to  the  systematic  uncertainty  in  estimating  the 
soot  mass  growth  rate  [9] . 

The  values  of  the  species  mole  fractions  at  a  given  axial  position  along  a  particle  path  were 
interpolated  from  the  measured  radial  locations.  This  procedure  was  done  on  the  same  seven 
particle  paths  in  order  to  accumulate  27  points.  When  locations  with  corrections  due  to  density 
changes  larger  than  that  allowed  in  this  were  eliminated,  22  points  remained.  The  locations  that 
were  eliminated  for  having  large  density  corrections,  also  corresponded  to  the  same  locations  that 
had  large  oxidation  corrections.  It  is  important  to  note  that  C2H2  and  CH4  are  the  only  hydrocarbon 
species  which  are  detected  throughout  the  soot  growth  region. 

Before  comparing  soot  growth  results  from  premixed  and  diffusion  flames,  several  things  should 
be  pointed  out.  The  soot  growth  region  in  the  premixed  flames  listed  in  Table  2  can  be  at 
significantly  greater  temperatures  than  the  bulk  of  the  soot  growth  region  in  the  diffusion  flame  of 
this  study.  The  temperature  range  of  some  of  the  premixed  flames  is  1600-1700K  [9].  In  the  laminar 
diffusion  flame  studied  here,  the  temperatures  where  soot  growth  occurs  is  lower,  generally  1300- 
1500K  [14],  though  several  points  near  the  end  of  growth  have  temperatures  above  1600  K.  The 
higher  temperatures  in  the  premixed  flames  could  lead  to  oxidation  corrections  that  are  important. 
However,  in  most  previous  studies,  the  calculation  of  the  Specific  Surface  Growth  Rate  (SSGR), 
equation  9,  and  the  rate  coefficient  for  the  premixed  flames,  kaHn ,  no  correction  for  oxidation  was 
made.  Therefore,  the  SSGR  and  rate  coefficients  of  this  study  presented  in  Table  2  have  only  been 
corrected  for  the  density  changes  which  occur  along  a  particle  path.  The  rate  coefficients  with 
oxidation  correction  included  are  not  appreciably  different  than  those  calculated  with  only  the 
density  correction.  In  addition,  no  density  correction  is  applied  to  the  premixed  flame  data  since 
the  temperature  in  the  postflame  zone  of  the  1-D  flames  is  relatively  constant. 

SSGR  =  kCzHi  [C2H2  ]  =  (9) 
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Table  2 


Specific  Surface  Growth  Rate  (SSGR),  Rate  Coefficient  (k)  and 
Equivalence  Ratio  (O)  from  selected  flames  burning  ethene  except 
as  noted 


Source 

SSGR 

(g-soot  cm"2-soot  s'1) 

Rate  Coefficient,  k 
(g-soot  cm_2-soot  s_1 
atm'1-C2H2) 

Equivalence 

Ratio, 

Harris  and  Weiner  [8] 

8xl0-6  to  6xl0"5 

NA 

2. 1-2.4 

Harris  and  Weiner  [9] 

7xl0'5  to  lxlO'4 

4xl0-4  to  4xl0-3 

2.28-2.82 

Harris  and  Weiner  [10]  a 

6xl0'6  to  3xl0-5 

7xl0"4  to  3xl0  3 

1.78-1.92 

Kent  and  Wagner  [16] 

5xl0"5  to  3X10"4 

NA 

b 

Vandsburger  et  al.  [44] 

lxlO'5  to  lxlO"4 

NA 

c 

Puri  et  al.  [29]  d 

3xl0"5  to  9xl0"5 

lxlO"3  to  3xl0"3 

1.12  to  1.77 

Present  Study 

3x10-5  to  5xl0-4 

6xl0-4  to  2x10-4 

0.71  to  4.3 

Richardson  [25] 

2x10-5  to  2xl0-4 

NA 

e 

Sunderland  et  al.  [52] 

2xl0-6  to  2xl0'3 

NA 

f 

NA  -  Not  Available 
NOTES  FOR  TABLE: 

a)  A  mixture  of  ethene  and  toluene  was  the  fuel  for  this  study. 

b)  The  SSGR  was  estimated  by  Harris  and  Weiner  [9]  from  Kent  and  Wagner  [16]  from  laminar 
ethene  diffusion  flames  on  a  Wolfhard-Parker  burner. 

c)  This  was  a  counterflow  ethene  diffusion  flame. 

d)  Puri  et  al.  [29]  uses  Fractal  Aggregate  Theory  for  the  aerosol  surface  analysis.  The  rate 
coefficient  quoted  is  calculated  from  the  SSGR  as  published  [29]  and  the  C2H2  mole  fraction 
from  the  experimental  section  of  this  study. 

e)  The  flames  of  Richardson  [25]  are  CH4  and  mixture  flames  with  CH4  and  either  C4H10,  C4H8, 
or  C4H6. 

f)  These  were  C2H2  diffusion  flames  below  atmospheric  pressure  and  includes  density  and 
oxidation  corrections  during  growth. 
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In  general,  the  values  of  SSGR  from  both  premixed  and  diffusion  flames  overlap  as  noted  in  Table 
2.  The  SSGR  from  diffusion  flames  are  greater  than  those  found  in  premixed  flames.  In  premixed 
flames,  it  was  estimated  that  uncertainties  in  the  surface  area,  index  of  refraction  of  soot,  and  the 
differentiation  procedure  could  lead  to  systematic  errors  of  a  factor  of  two  to  three  in  the  calculated 
rate  coefficients  [9].  The  rate  coefficients  as  calculated  in  premixed  and  diffusion  flames  follow 
the  same  pattern  as  the  SSGR.  That  is.  the  ranges  of  calculated  values  overlap  with  the  maximum 
value  from  this  diffusion  flame  exceeding  the  maximum  value  found  in  premixed  flames. 

Recent  work  introducing  Fractal  Aggregate  Theory  (FAT)  to  the  measurement  of  the  soot  aerosol 
has  enabled  an  improvement  in  the  accuracy  of  the  determination  of  the  soot  surface  area  [291 .  The 
volume  fraction  along  PP  6  that  is  calculated  by  the  scattering/extinction  data  and  analyzed  by 
Mie  theory  differs  by  approximately  20%  from  that  analyzed  by  the  fractal  approach.  However,  the 
surface  density  (soot  area  per  unit  volume)  as  calculated  with  FAT  is  over  a  factor  of  3  greater  than 
the  Mie  calculation  as  a  result  of  the  more  detailed  treatment  of  shape  of  the  scattering  particles 
and  the  additional  information  from  the  TEM  measurements.  In  a  premixed  flame  [8]  the  surface 
density  from  optical  measurements  gave  a  soot  area  almost  twice  that  found  with  a  BET  analysis 
of  that  soot.  If  the  errors  in  measuring  the  areas  as  described  above  are  systematic  and  consistent, 
and  the  FAT  and  BET  areas  are  the  appropriate  areas  for  normalizing  the  soot  growth  rate,  then  the 
net  effect  would  be  to  make  the  SSGR  and  rate  coefficient  for  premixed  and  diffusion  flames  much 
more  similar. 

Another  possible  explanation  for  the  greater  growth  rate  in  diffusion  flames  could  be  a  greater 
contribution  to  the  growth  rate  due  to  condensation  of  PAH  species  [57-60].  The  two  dimensional 
nature  of  this  flame  could  even  enhance  this  effect.  That  is,  PAH  species  that  form  in  the  center 
region  of  the  flame  could  diffuse  radially  outwards  and  be  consumed  by  a  growth  reaction  with  the 
soot  surface  or  oxidized.  Therefore,  the  soot  along  a  given  particle  path  is  not  limited  by  the  carbon 
which  is  initially  present  on  that  particle  path. 

The  range  of  equivalence  ratios  where  measurements  of  growth  rate  are  available  is  broader  in 
this  diffusion  flame  than  in  premixed  flames.  This  is  due  to  the  nature  of  this  diffusion  flame 
where  the  inception  of  particles  can  occur  at  equivalence  ratios  higher  than  premixed  flame 
experiments  and  the  particles  are  swept  into  regions  with  progressively  decreasing  equivalence 
ratio.  Not  enough  measurements  are  available  to  determine  any  equivalence  ratio  effect  on  the 
growth  rate.  In  general,  the  study  of  soot  growth  in  premixed  flames  is  limited  to  a  range  of 
equivalence  ratios  which  will  incept  particles,  while  not  producing  so  many  particles  that  the 
burner  surface  becomes  contaminated  from  the  pyrolysis  products  and  soot  particles. 

SUMMARY 


Measurements  of  the  spatial  distribution  of  species  concentrations  in  a  laminar  diffusion  flame 
with  significant  soot  loadings  have  been  carried  out.  The  shape  of  the  radial  profiles  of 
intermediate  hydrocarbon  species  indicates  the  importance  of  temperature  in  their  formation. 
Q>H2,  C4H6  and  C6H6,  along  with  CH4,  are  the  most  prominent  pyrolysis  products  and,  except  for 
CH4,  are  considered  to  play  important  roles  in  the  formation  and  growth  of  soot. 
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These  spatial  measurements  have  been  integrated  with  previous  flowfield  and  soot  measurements 
to  establish  the  mechanism  controlling  the  cessation  of  soot  surface  growth  in  this  laminar 
diffusion  flame.  In  particular,  the  measurements  of  C2H2  and  other  hydrocarbons  have  revealed 
that  soot  surface  growth  ceases  due  to  the  depletion  of  hydrocarbons  and  not  due  to  the  loss  of  soot 
particle  reactivity.  This  observation  is  in  agreement  with  recent  work  by  Honnery  and  Kent  [43], 
who  have  shown  that  in  long  laminar  diffusion  flames,  soot  growth  can  be  extended  to  very  long 
residence  times.  In  addition,  C2H2  and  CH4  are  the  only  hydrocarbons  detected  throughout  the 
region  of  soot  growth. 

Calculation  of  a  rate  coefficient  for  soot  surface  growth  from  C2H2  can  now  be  done  for  this 
diffusion  flame.  The  observed  rate  coefficient  in  this  diffusion  flame  is  in  agreement  with  values 
found  in  premixed  flame  studies.  This  result  suggests  that  the  mechanism  of  soot  surface  growth 
dominant  in  premixed  flames  is  also  dominant  in  this  diffusion  flame. 
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NOMENCLATURE 


^  Particle  Reactivity  Decrease  Factor  is*1) 

F,  Mass  Spectrometer  Calibration  Factor  for  species  i  (V-1) 

fv  Soot  Volume  Fraction  (em3-soot  cm’3) 

k ,  kC  H ^  Surface  Growth  Rate  Coefficient  (g-soot  cm'2-soot  s’1  atm"1-C2H2) 

kSG  Particle  Reactivity  Decrease  Factor  (s’1) 

Ns  Concentration  of  Active  Sites  (sites  cm’2) 

Psoot  Soot  Density  (1.8  g  cm-3) 

S  Soot  Surface  Area  per  unit  volume  (cm2-soot  cm-3) 

Si  Mass  Spectrometer  Signal  for  Species  i  (V) 

t  Time  (s) 

Ta  Activation  Temperature  (K) 

XL  Mole  Fraction  of  species  i  (-) 

CO  Mass  growth  rate  of  soot  (g-soot  cm-3  s_1) 

MS()of  Mass  of  soot  per  unit  volume  (g-soot  cm*3) 
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Solenoid 


Figure  1  Electromechanical  Sonic  (EMS)  probe  Schematic 


FIGURE  2  Cylinder  system  for  compression  and  storage  of  gas  sample  in  preparation  of 
injection  into  gas  chromatograph  (after  Leonard  [61]).  (1)  Inlet  to  sample 
accumulation  cylinder,  (2)  Shut-off  valve  to  probe,  (3)  Exhaust  valve  to  vacuum 
pump,  (4)  and  (5)  three-way  valves.  A  and  B  indicate  piston  position  prior  to  and 
following  helium  pressurization. 
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FIGURE  3  Schematic  of  vacuum  chambers  with  mass  spectromenter. 


Radial  Position  (mm) 


FIGURE  4  Mole  fraction  as  a  function  of  radial  position  for  several  heights  above  the 
burner  outlet  for  (a)  C2H4  and  (b)  C2H2.  Comparisons  can  be  made  between 
data  from  the  gas  chromatograph  (open  symbols)  and  the  mass  spectrometer 
(filled  symbols). 


C  H  Mole  Fraction,  ppm 


Radial  Position  (mm) 


FIGURE  5  Mole  fraction  as  measured  by  the  gas  chromatograph  as  a  function  of  radial 
position  for  several  heights  above  the  burner  outlet  for  (a)  C4H6  and  (b)  C6H6 
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3  Soot  volume  fraction.  fv,  and  the  mole  fraction  of  C2H2  as  a  function  of  tin 
(a)  particle  path  6,  (b)  particle  path  3  and  (c)  particle  path  1.  (see  text  for 


Soot  Volume  Fraction,  f  Soot  Volume  Fraction,  f  Soot  Volume  Fraction,  f 
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FIGURE  8  Temperatures  along  particle  paths  1,  3  and  6.  The  open  symbols  are  locations 
with  no  soot  present.  Locations  with  filled  symbols  indicate  the  presence  of 
soot.  The  circled  symbol  on  each  particle  path  has  the  maximum  local  soot 
volume  fraction  along  that  particle  path. 
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Mole  Fraction 


Radial  Position  (mm) 


FIGURE  9  Mole  fraction  of  02  as  measured  by  the  gas  chromatograph  as  a  function  of 
radial  position  for  selected  heights  above  the  burner  outlet. 
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Abstract 


A  model  of  a  laminar,  ethene  diffusion  flame  has  been  developed  and  compared  with 
measurements.  A  system  of  elementary  reactions  is  used  to  describe  the  gas  phase  Q  and 
C2  chemistry.  The  model  incorporates  a  simple  description  of  the  growth  of  soot  which 
performs  a  good  job  of  reproducing  the  amount  of  soot  in  the  flame  and  the  transition  from 
non-sooting  to  sooting  conditions.  Stable  gas  phase  species  have  been  measured  in  the 
flame  with  a  combination  of  gas  chromatography  and  mass  spectrometry.  Laser-induced 
fluorescence  was  used  to  measure  the  OH  concentration.  The  model  predicted  most  of  the 
species  quite  well.  It  was  found  that  the  assumption  of  equilibrium  concentrations  for  the 
minor  radical  species  was  grossly  in  error,  particularly  near  the  end  of  the  flame  where 
temperatures  were  low.  Calculated  radical  species  were  found  in  concentrations  that  were 
orders  of  magnitude  greater  than  equilibrium.  Comparisons  with  measured  OH 
concentrations  supported  this  observation. 


Introduction 


Previous  studies  of  a  laminar  ethene  -  air  diffusion  flame  have  suggested  that  the  modeling 
of  soot  formation  with  the  assumption  of  equilibrium  chemistry  may  underestimate 
significantly  the  super-equilibrium  concentrations  of  important  radicals  such  as  O  and 
OH [1], [2] /These  species  are  important  oxidants  of  soot,  particularly  near  the  tip  of 
diffusion  flames  where  soot  can  escape.  The  original  code  that  has  been  developed  at  U.  C. 
Davis  made  use  of  the  equilibrium  assumption[2],[3]  in  order  to  determine  the  mole 
fractions  of  all  the  species  in  the  flame.  Finite  rate  chemistry  is  required  in  order  to  relax 
this  assumption.  An  entirely  new  code  was  written  to  accommodate  this  modification. 
Detailed  chemistry  has  been  incorporated  into  a  parabolic  code  that  solves  the  boundary 
layer  form  of  the  conservation  equations  for  mass,  momentum  and  energy  conservation.  A 
much  more  accurate  representation  of  radical  concentrations  has  been  achieved. 
Furthermore,  interesting  questions  that  are  related  to  the  interaction  of  soot  with  flame 
chemistry  can  now  be  investigated.  The  results  of  the  modeling  have  been  compared  with 
the  results  of  an  experimental  program  at  The  Pennsylvania  State  University. 

Experimental  Method 

A  laminar  ethene/air  diffusion  flame  was  established  on  a  co-annuiar  burner  similar  in 
design  to  that  used  in  previous  studies [4]  and  only  a  brief  description  will  be  given  here. 
This  burner  consists  of  a  11.1  mm  diameter  brass  fuel  tube  surrounded  by  a  102  mm 
diameter  air  annulus.  A  brass  cylindrical  chimney  was  used  to  shield  the  flame  from 
laboratory  air  currents  with  access  for  sampling  provided  by  slots  machined  in  the  chimney 
wall.  The  ethene  (C2H4)  and  air  flow  rates  were  3.85  cm3  s_1  and  713  cm3  s'1 
respectively,  which  were  measured  with  calibrated  rotameters.  Ethene  with  a  stated  purity 
of  99.5%  was  used:  dried,  filtered  air  was  supplied  by  an  in-house  compressor. 

Radial  profile  measurements  of  species  concentrations  have  been  obtained  throughout  the 
flame  using  an  intrusive  quartz  microprobe  sampling  approach.  Quartz  microprobes  have 
been  used  in  many  studies  to  perform  the  extraction  of  gas  samples  from  premixed  flames 
[5], [6]  and  diffusion  flames[7],  [8],  [9],  [10],  [11],  Until  recently  their  use  has  been 
limited  to  regions  where  the  soot  volume  fraction  is  low  to  avoid  clogging  the  probe  orifice. 
A  variation  of  the  standard  quartz  microprobe  design,  referred  to  as  an  Electro-Mechanical 
Sonic  (EMS)  probe,  was  developed  by  Puri  [11]  and  enabled  sampling  in  flame  regions 
where  volume  fractions  were  on  the  order  of  10“/  This  probe  has  a  fiber  which  extends 
throush  the  orifice  and  is  attached  to  a  solenoid  core  with  springs  on  either  side  of  the  core. 
A  solenoid  is  placed  around  the  outside  of  the  quartz  tube  and  the  current  to  the  solenoid  is 
interrupted  at  about  30  Hz.  providing  for  constant  motion  of  the  fiber  in  the  orifice,  thus 
keepins  an  annular  region  open  for  gas  sampling.  The  fiber  is  =125  |im  in  diameter  and  is 
made  from  sapphire,  while  the  orifice  diameter  is  =190  pm.  The  annular  sampling  region 
has  an  area  equivalent  to  a  =130  pm  diameter  orifice. 

A  2 as  chromatograph  and  a  mass  spectrometer  (MS)  were  used  to  analyze  the  samples.  The 
use  of  these  two  techniques  offered  an  independent  verification  of  the  concentration  profile 
measurements  as  well  as  providing  specific  species  measurement  capability  as  described 
below.  In  general,  the  two  analytical  approaches  yielded  agreement  to  within  ±20%. 
Because  of  difficulties  in  detecting  H2  with  the  GC.  the  MS  measurements  provide  the 
determinations  for  this  species.  Since  all  of  the  data  reported  in  this  paper  were  obtained 
usine  the  MS.  only  this  instrumentation  will  be  described  in  further  detail. 
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The  mass  spectrometer  employed  an  Extrel  model  7-324-9  quadrupole  mass  filter  with 
mass  range  of  1-500  amu  for  separation  of  ions  and  an  electron  multiplier  for  detection. 
Calibration  of  the  MS  employed  a  method  used  by  Korobeinichev  [12]  which  is 
mathematically  similar  to  that  of  Bittner  [13]  and  Crowhurst  and  Simmons[9].  The  method 
relies  on  the  signal  response  relative  to  other  species  remaining  constant.  That  is,  given  the 
same  mixture,  the  signals  of  two  gases  will  always  be  in  the  same  ratio  and  is  therefore 
independent  of  mass  flow  through  the  orifice  due  to  any  possible  orifice  clogging  or 
variations  in  density.  Secondly,  this  method  assumes  that  the  species  measured  constitute 
all  of  the  gases  present  and  they  are  normalized  to  a  mole  fraction  of  one. 

A  limitation  of  the  MS  in  this  flame  is  the  coincidence  of  C2H4>  N2  and  CO  at  mass  28. 
For  C2H4,  the  fragment  at  mass  27  was  used  as  an  independent  measurement.  Calibration 
of  C2H4  at  both  masses  27  and  28  was  necessary  so  that  the  signal  contribution  of  C2H4  at 
mass  28'could  be  subtracted  from  the  mass  28  signal. 

Hydroxyl  radical  (OH*)  concentration  measurements  were  obtained  in  a  previous  study  [1]. 
The  OH*  measurements  were  obtained  using  the  output  from  a  Nd:YAG  pumped  dye  laser 
which  was  frequency  doubled  to  produce  a  UV  beam.  The  fluorescence  signal  was 
collected  at  90°  with  an  intensified  cooled  CCD  camera.  The  laser  wavelength  was  tuned  to 
either  278.83  nm  or  283.55  nm.  corresponding  to  S  12(8)  and  Qi(8)  lines  of  the  A^X+  «- 
X2rii  (1,0)  band  of  OH*,  respectively.  Filters  attenuated  elastically  scattered  light  from 
soot  particles  and  transmitted  the  (0.0)  and  (1,1)  emission  bands  of  OH*.  Excitation  of  the 
Si2(8)  line  of  OH*  was  used  in  regions  of  low  soot  particle  concentration  while  the 
stronger  Qj(8)  line  was  utilized  when  soot  concentrations  were  high  (see  Puri  et  al.  [1] ). 

Numerical  Model 

The  most  general  formulation  of  the  governing  equations  consists  of  the  compressible 
Navier-Stokes  equations,  an  equation  to  describe  the  conservation  of  energy  and  species 
transport  equations  that  are  elliptic  in  nature.  The  computer  time  to  solve  these  fully  elliptic 
equations  is  large  [14] .  A  number  of  assumptions  are  introduced  in  the  present  model  that 
permit  a  significant  saving  in  computer  time  without  a  serious  impact  on  accuracy.  The 
boundary  layer  forms  of  the  governing  equations  are  solved  with  zero  pressure  gradient. 
Hence,  the  equations  can  be  solved  with  a  marching  method.  It  is  expected  that  the 
boundary  layer  assumption  will  introduce  the  greatest  errors  in  the  vicinity  of  the  nozzle 
exit.  The  flame  that  has  been  investigated  in  this  study  is  relatively  long,  of  the  order  of  7 
nozzle  diameters.  Therefore,  it  is  expected  that  the  errors  will  not  be  large  in  the  region  of 
most  interest  in  relation  to  soot  oxidation,  radiation  and  their  impact  on  gas  phase 
chemistry. 

The  equations  are  solved  in  primitive  variable  form.  Hence,  the  steady  state  continuity 
equation 
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is  solved  alone  with  the  momentum  equation 
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The  energy  equation  is  written  in  terms  of  the  temperature  as 
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in  which  the  complete  diffusion  velocity  of  species  n  is  Vn 

Yn  =  vn  +  Vj  +  ynCOT  [4] 

The  divergence  of  the  radiative  heat  flux  is  represented  by  qR.  The  normal  diffusion 
velocity  of  species  n  is  given  by 
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and  the  contribution  of  the  thermal  diffusion  phenomenon  to  the  total  diffusion  velocity  is 
given  by 
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The  thermal  diffusion  ratio,  Kn,  was  obtained  from  CHEMKIN  [15].  A  correction 
diffusion  velocity,  V^0TT 


yncorT  =  -  £  Yn  (vn  +  VT)  [7] 

ensures  that  net  diffusive  flux  of  all  gaseous  species  is  zero.  An  equation  is  solved  for  N-l 
species  mass  fractions 
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with  the  Nth  species  being  N2-  Finally,  an  equation  is  solved  for  the  soot  mass  fraction. 
Ysoot 
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In  this  equation,  the  thermophoretic  velocity  is  found  from  the  appropriate  expression  for  a 
free  molecular  aerosol  [16] 
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The  source  term  in  equation  (9)  includes  the  contributions  of  soot  nucieation  (wnucj),  soot 
surface  growth  (wg)  and  soot  oxidation  (w0) 

^(^SOOP  4-  f)  —  Wnucj  +  Wg  —  w0  [11] 

An  equation  for  soot  particle  number  density  is  not  included.  Rather,  an  average  number 
density  is  assumed.  The  justification  for  this  assumption  is  discussed  by  Kennedy  et  ai.[3] 
The  rates  of  reaction  for  the  surface  growth  and  soot  oxidation  terms  are  identical  to  those 
used  by  Kennedy  et  al.[3] . 

Method  of  Solution 

Initial  Conditions 

Initial  conditions  for  all  the  major  species  as  well  as  the  temperature  profiles  are  required  in 
order  to  start  the  calculation.  These  initial  conditions  are  obtained  by  first  starting  the 
calculations  with  an  equilibrium  chemistry  model  with  the  Shvab-Zeldovich  approximation 
[17]  so  that  the  mixture  fraction  can  be  used  to  describe  the  scalar  field.  This  is  achieved 
by  generating  a  table  of  mass  fractions  of  all  species,  density  and  temperature  as  functions 
of  mixture  fractions  with  the  chemical  equilibrium  code  STANJAN  [18].  The  momentum, 
continuity  and  mixture  fraction  equations  are  then  solved  simultaneously.  Once  the  mixture 
fraction  profile  is  known,  all  the  other  dependent  variables  can  be  obtained  by  using  the 
table.  Numerical  trials  have  established  that  reasonable  initial  conditions  for  the  finite  rate 
chemistry  calculation  can  be  obtained  by  using  the  mixture  fraction/equilibrium  calculation 
up  to  0. 15  nozzle  diameters. 


Chemical  Mechanism 

A  "skeletal"  mechanism  has  been  used  to  describe  concisely  the  Cl  and  the  C2  reactions  for 
ethene  oxidation.  The  methane  mechanism  of  Smooke  et  al  [19]  has  been  found  to  be 
numerically  efficient  and  it  has  been  implemented  in  the  current  code  with  the  rates  for  all 
the  Cj  chemistry  from  Reference  [19] .  It  contains  a  sufficient  degree  of  detail  to  describe 
methane  oxidation  for  the  purposes  of  this  project.  The  C2  chemistry  is  described  by  a 
series  of  reactions  that  convert  C2H4  to  C2H3  and  C2H2.  The  immediate  source  of  the 
rates  for  the  Ci  reactions  is  the  paper  by  Frenklach  et  al.  [20]  ;  primary  sources  of  the  rate 
data  may  be  found  in  that  reference.  The  entire  reaction  mechanism  is  shown  in  Table  1 
and  is  composed  of  62  reactions  involving  24  species.  The  important  oxidation  steps  by  O, 
O2  and  OH  are  included.  The  CO  oxidation  steps  and  the  H2  -  O2  steps  are  also  included. 


5 


Table  1 

Chemical  Reaction  Mechanism  in  the  form  k  =  A  TP  e  -E/RT  in  units  of  moles,  cubic 
centimetres,  seconds,  Kelvin  and  calories/mole. 

A  PH 

3.00E+13  0  .  0 


1.  CH2+OH  -4  CH20+H 

2.  CH2+02  -4  CO+OH+H 

3.  C2H+H2  -4  C2H2+H 

4.  C2H+02  -4  HCCO+O 

5.  HCCO+H  -4  CH2+CO- 

6.  C2H2+0  -»  CH2+CO 

7.  C2H2+OH  -4  HCCOH+H 

8.  CH2CO+H  -4  HCCO+H2 

9.  CH2CO+OH  -4  HCC0+H20 

10.  C2H3+0  -4  CH2CO+H 

11.  C2H3+02  -4  CH20+HC0 

12.  C2H3+M  -4  C2H2+H+M 

13.  C2H4+OH  -4  C2H3+H20 

14.  C2H4+M  -4  C2H2+H2+M 

15.  C2H4+H  -4  C2H3+H2 

16.  C2H4+0  -4  CH3+HCO 

17.  CH4+M  -4  CH3+H+M 

18.  CH4+02  -4  CH3+H02 

19.  CH4+H -4  CH3+H2 


3.10E+13 

0 

0 

1.10E+13 

0 

2866.0 

6.02E+1 1 

0 

0 

1.50E+14 

0 

0 

7.81E+03 

2.8 

502.0 

5.06E+05 

2.3 

13496.0 

3.00E+13 

0 

8599.0 

1.00E+13 

0 

2627.0 

3.00E+13 

0 

0 

4.00E+12 

0 

-239.0 

2.00E+38 

-7.2 

50567.0 

3.00E+13 

0 

2986.0 

2.60E+17 

0 

79230.0 

3.16E+1 1 

7 

8002.0 

1.60E+08 

1.4 

525.0 

1.00E+17 

0 

86000.0 

7.90E+13 

0 

56000.0 

2.20E+04 

3.0 

8750.0 

20.  CH4+0  -+  CH3+OH 

21.  CH4+OH  -+  CH3+H20 

22.  CH20+OH  HC0+H20 

23.  CH20+H  -»  HCO+H2 

24.  CH20+M  -+  HCO+H+M 

25.  CH20+O  -4  HCO+OH 

26.  HCO+OH  -+  C0+H20 

27.  HCO+M  -+  H+CO+M  - 

28.  HCO+H  -+  CO+H2 

29.  HCO+O  -+  OH+CO 

30.  HC0+O2  ->  H02+C0 

31.  CO+O+M  -+  C02+M 

32.  CO+OH  -+  C02+H 

33.  C0+02  -+  C02+0 

34.  CH3+02  -+  CH30+0 

35.  CH30+M  -+  CH20+H+M 

36.  CH30+H  -+  CH20+H2 

37.  CH30+0H  -+  CH20+H20 

38.  CH30+0  -+  CH20+0H 

39.  CH30+02  -+  CH20+H02 

40.  CH3+02  ->  CH20+0H 

41.  CH3+0  — »  CH20+H 


1.60E+06 

2.4 

7400.0 

1.60E+06 

2.1 

2460.0 

7.53E+12 

0 

167.0 

3.31E+14 

0 

10500.0 

3.31E+16 

0 

81000.0 

1.81E+13 

0 

3082.0 

5.00E+12 

0 

0 

1.60E+14 

0 

14700.0 

4.00E+13 

0 

0 

1.00E+13 

0 

0 

3.00E+12 

0 

0 

3.20E+13 

0 

-4200.0 

1.51E+07 

1.3 

-758.0 

1.60E+13 

0 

41000.0 

7.00E+12 

0 

25652.0 

2.40E+13 

0 

28812.0 

2.00E+13 

0 

0 

1.00E+13 

0 

0 

1.00E+13 

0 

0 

6.30E+10 

0 

600.0 

5.20E+13 

0 

34574.0 

6.80E+13 

0 

0 

42.  CH3+0H  ->  CH20+H2 

43.  H02+C0  C02+0H 

44.  H2+02  20H 

45.  0H+H2  H20+H 

46.  H+02  OH+O 

47.  0+H2  OH+H 

48.  H+02+M  H02+M 

49.  H+02+02  H02+02- 

50.  H+02+N2  H02+N2 

51.  0H+H02  ->  H20+02 

52.  H+H02  -¥  20H 

53.  0+H02  02+0H 

54.  20H  -4  0+H20 

55.  H2+M  H+H+M 

56.  02+M  -4  O+O+M 

57.  H+OH+M  -4  H20+M 

58.  H+H02  -4  H2+02 

59.  H02+H02  -4  H202+02 

60.  H202+M  -»  OH+OH+M 

61.  H202+H  — >  H02+H2 

62.  H202+0H  -»•  H20+H02 


7.50E+12 

0 

0 

5.80E+13 

0 

22934.0 

1.70E+13 

0 

47780.0 

1.17E+09 

1.3 

626.0 

2.20E+14 

0 

16800.0 

1.80E+10 

1.0 

8826.0 

2.10E+18 

-1.0 

0 

6.70E+19 

-1.4 

0 

6.70E+19 

-1.4 

0 

5.00E+13 

0 

1000.0 

2.50E+14 

0 

1900.0 

4.80E+13 

0 

1000.0 

6.00E+08 

1.3 

0 

2.23E+12 

.5 

92600.0 

1.85E+11 

.5 

95560.0 

7.50E+23 

-2.6 

0 

2.50E+13 

0 

700.0 

2.00E+12 

0 

0 

1.30E+17 

0 

45500.0 

1.60E+12 

0 

3800.0 

1.00E+13 

0 

1800.0 

Numerical  Scheme 


The  governing  differential  equations  are  written  in  a  discrete  form  by  using  a  second  order 
accurate  scheme  for  the  diffusion  terms  and  first  order  accurate  scheme  for  the  convective 
terms  with  the  delta  formulation.  The  non-linear  terms  (convective  and  rate  terms)  are 
linearized  by  using  Newton's  method.  The  resulting  Jacobian  matrices  are  determined 
numerically  in  the  manner  described  by  Smooke  [21] .  The  Jacobian  does  not  need  to  be 
updated  continuously  since  the  changes  from  one  step  to  the  next  step  are  small.  We  found 
that  the  evaluation  of  the  Jacobian  at  every  fifth  spatial  step  yielded  the  same  results  as  an 
evaluation  at  every  step. 

The  resulting  sets  of  coupled  algebraic  equations  are  solved  implicitly  and  simultaneously 
by  the  usage  of  a  block  solver  with  iterations  to  account  for  the  non-linear  term.  Typically, 
8  iterations  per  step  are  required  to  converge  the  equations  to  a  local  solution.  The  step  size 

(Ax)  is  determined  by  limiting  the  maximum  changes  of  all  the  dependent  variables  by  no 
more  than  10  percent.  If  the  maximum  changes  in  any  dependent  variable  exceed  this  limit, 

the  solution  is  rejected  and  a  smaller  Ax  is  used  and  the  calculation  is  redone.  The  equations 
are  assumed  to  be  solved  when  the  residual  approaches  machine  zero.  If  a  local 

convergence  cannot  be  reached,  a  smaller  Ax  is  used  and  the  solution  is  recalculated. 

Grid  adaptation  is  used  in  order  to  resolve  the  high  gradient  region  in  the  flow  field  with  a 
minimum  number  of  grid  points.  This  is  done  by  adapting  the  grid  distribution  to  the 
temperature  gradient.  The  results  shown  in  this  report  are  obtained  by  using  only  61  points 
in  the  transverse  direction.  All  the  properties  of  the  fluid,  the  diffusion  coefficients  of  all 
species  and  the  chemical  source  term  for  each  species  were  obtained  by  using  CHEMKIN 
[15]  . 

Results 

Solutions  for  a  non-sooting  and  a  sooting  flame  are  presented.  The  predicted  integrated 
soot  volume  fractions  are  shown  in  Fig.  1 ;  they  are  a  measure  of  the  total  amount  of  soot  at 
a  given  axial  plane.  Experimental  data  from  Santoro  et  al.  [22]  are  also  shown  in  Figure  1  . 
The  agreement  between  the  model  predictions  and  the  measurements  of  the  integrated  soot 
volume  fraction  is  considered  quite  good.  Bum  out  of  soot  at  the  end  of  the  flame  is 
evident.  Generally,  the  soot  model  tends  to  underpredict  the  measured  values  of  soot 
volume  fraction.  However,  the  spatial  distribution  of  soot  is  reproduced  very  well 
including  the  location  of  the  peak  integrated  soot  volume  fraction.  Although  it  is  not  shown 
in  this  paper,  the  model  manages  to  reproduce  the  transition  from  a  non-sooting  flame  to  a 
sooting  flame  quite  well  as  the  fuel  flow  rate  is  increased.  The  agreement  that  is  shown  in 
Fig.  1  ,  along  with  more  detailed  radial  profiles  of  soot  volume  fractions  which  are  not 
shown,  provide  the  basis  for  confidence  in  the  application  of  this  relatively  simple  soot 
model  to  a  calculation  of  detailed  chemistry  in  an  ethene  flame. 

In  addition  to  a  simple  soot  model,  the  code  includes  a  radiation  model,  in  which  an 
optically  thin  model  is  adopted  for  simplicity.  A  comparison  of  the  measured  and  predicted 
temperature  fields  offers  a  measure  of  its  accuracy.  Figure  2  presents  radial  profiles  of  the 
measured  and  predicted  temperatures  at  two  axial  locations  in  the  non-sooting  flame.  The 
drop  in  the  peak  temperature  at  the  more  downstream  location  is  evidence  of  the  heat  loss 
due  to  luminosity.  The  agreement  in  the  temperature  field  is  considered  good,  although  the 
radiation  model  leads  to  underprediction  of  the  radiative  heat  loss,  in  part  because  gas 
radiation  is  ignored  and  also  because  the  soot  volume  fractions  are  underpredicted.  Similar 
agreement  is  evident  in  the  velocity  field  although  these  results  are  not  shown. 


The  fuel  mole  fractions  in  the  non-sooting  flame  are  shown  in  Fig.  3  at  a  location  that  is  7 
mm  from  the  fuel  nozzle.  A  profile  of  C2&2  mole  fractions  at  the  same  location  is  shown  in 
Fig.  4  .  There  is  a  general  tendency  to  overpredict  the  measured  mole  fractions  of  C2H4 
and  C2H2  in  the  central  region  of  the  flame  while  the  model  results  show  steeper  gradients 
as  the  radial  location  increases.  The  finite  spatial  resolution  of  the  EMS  probe  contributes 
to  the  observed  differences  in  this  region  since  species  concentration  gradients  in  the  flame 
are  steep.  Further  improvement  in  the  comparisons  between  model  and  experiments  will 
require  higher  spatially  resolved  measurements  and  consequently  the  current  agreement  is 
considered  satisfactory  for  model  validation.  The  loss  of  C2H2  onto  soot  is  not  accounted 
for  at  present  in  the  balance  equation  for  acetylene.  This  sink  of  acetylene  could  reduce  the 
predicted  C2H2  concentrations,  particularly  "further  downstream  where  soot  mass  growth 
rates  are  more  intense.  The  hydrogen  mole  fraction  profile  at  7  mm  is  shown  in  Fig.  5  . 
The  calculated  values  exhibit  reasonably  good  agreement  with  the  experimental  data  in  view 
of  the  difficulty  in  measuring  this  species.  At  a  location  20  mm  from  the  fuel  nozzle,  the  H2 
distribution  has  moved  towards  the  centerline  as  the  flame  reaction  zone  approaches  the 
centerline.  The  results  indicate  that  the  general  structure  of  this  luminous  flame  is 
reproduced  quite  well  by  the  numerical  model. 

Hydroxyl  mole  fractions  have  been  measured  with  laser  induced  fluorescence.  The 
calculated  hydroxyl  mole  fractions  are  compared  with  the  measurements  at  7  mm  and  70 
mm  in  Fig.  6  .  The  model  tends  to  overpredict  the  OH  mole  fractions.  This  tendency  may 
be  due,  in  part,  to  the  absence  of  a  numerical  sink  for  OH  on  soot  particles.  The  model 
does  not  incorporate  this  aspect  of  the  physics  at  present.  A  full  coupling  of  the  flame 
chemistry  and  the  soot  chemistry  will  be  attempted  in  the  future  with  a  modified  soot  model 
that  will  incorporate  elementary  steps  for  soot  mass  growth  and  oxidation.  It  must  also  be 
bom  in  mind  that  T.TF  measurements  involve  uncertainties  that  are  related  to  quenching 
corrections  and  other  errors.  The  errors  were  estimated  to  be  ±50%  [1]. 


The  preceding  comparison  of  numerical  and  experimental  results  provides  grounds  for 
some  confidence  in  the  results  of  the  calculations.  The  use  of  a  detailed  kinetics  model  with 
a  model  for  soot  formation  and  radiation  permits  us  to  examine  questions  regarding  the 
impact  of  the  soot  process  on  flame  chemistry  in  a  maimer  that  experiments  alone  cannot 
attempt.  For  example,  calculations  of  turbulent  flames  require  a  model  for  the 
thermochemistry.  A  convenient  and  very  common  approach  is  to  use  chemical  equilibrium. 
It  is  instructive  to  examine  the  accuracy  of  this  assumption  when  there  is  substantial 
radiative  heat  loss  from  a  hydrocarbon  flame.  The  non-equilibrium  effect  is  greatest  for  the 
minor,  radical  species  such  as  O  and  OH.  These  species  are  of  particular  interest  in  the 
formation  and  oxidation  of  pollutants  such  as  NO  and  soot.  The  ratio  of  computed  OH 
mole  fractions  to  equilibrium  mole  fractions  at  the  local  computed  temperature  were 
calculated  at  40  mm  from  the  nozzle.  At  a  radius  of  r/D  =  0.745, 


^^numerical 


O 


=  22  while 


numerical 


=  478.  This  location  is  to  the  air  side  of  the  reaction 


OHeq  —  —  Ogq 

zone  at  a  mixture  fraction  of  around  0.03  (stoichiometric  mixture  fraction  is  0.064).  It  is 
apparent  that  large  super-equilibrium  concentrations  can  be  achieved  as  a  result  of  the 
depressed  temperatures  in  this  radiating  flame.  Puri  et  al.  [1]  found  that  the  OH  super¬ 
equilibrium  ratio  was  9.4  at  the  peak  OH  location  in  this  flame  at  a  height  of  70  mm.  The 
results  suggest  that  assumption  of  equilibrium  chemistry  is  grossly  inadequate  for  the 
prediction  of  the  minor  species. 


Conclusions 

The  interacuon  of  soot  and  the  flame  chemistry  has  been  modeled  with  a  relatively  simple 
fluid  and  chemical  scheme.  A  comparison  of  the  numerical  results  with  detailed 


experimental  results  indicates  that  the  flame  structure  is  reproduced  quite  well.  Radicals 
such  as  OH  and  O  are  predicted  to  be  present  in  concentrations  that  are  significantly  greater 
than  equilibrium.  The  degree  of  super-equilibrium  is  a  function  of  the  temperature  and, 
hence,  the  soot  loading  in  the  flame.  Future  work  requires  the  addition  of  sink  terms  in  the 
species  equations  for  OH  and  C2H2.  It  shoulcfbe  possible  to  treat  other  fuels  by  including  a 
small  number  of  additional  reactions  for  the  pyrolysis  of  the  parent  fuel  to  C2  species. 
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Nomenclature 

A  pre-exponential  factor  in  rate  coefficient 

Cp  specific  heat  at  constant  pressure 

D  nozzle  diameter 

Dn  effective  diffusivity  of  species  n  in  the  mixture 

E  activation  energy 

f  mixture  fraction 
a  gravitational  acceleration 

hn  specific  enthalpy  of  species  n 

k  rate  coefficient 

Mn  molecular  weight  of  species  n 

N  number  of  species 

qR  divergence  of  the  radiative  heat  flux  vector 

R  gas  constant 

r  radial  distance 

S  source  or  sink  term  for  soot  equation 

T  temperature 

u  velocity'  in  the  axial  direction 

v  velocity  in  the  radial  direction 

vn  normal  diffusion  velocity  of  species  n 

vT  thermal  diffusion  velocity  of  species  n 

n 

V™  correction  to  diffusion  velocities 
n 

Vj  thermophoretic  velocity  of  soot  particles 

wn  reaction  rate  of  species  n 
Xn  mole  fraction  of  species  n 


X 


distance  along  the  axis 
Yn  mass  fraction  of  species  n 

P  exponent  on  temperature  in  rate  coefficient 

Kn  thermal  diffusion  ratio 

A  thermal  conductivity 

fi  absolute  viscosity 

v  kinematic  viscosity 

p  density 

Subscripts 

eq  equilibrium 

g  soot  mass  growth 

n  soot  nucleation 

numerical  calculated  with  computer  code 

o  soot  oxidation 
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Figure  Captions 

Fig.  1  Measured  and  predicted  soot  volume  fractions  integrated  across  a  non- 

sooting  flame 

-  calculations 

■  experiments 

Fig.  2  Measured  and  predicted  temperature  profiles 

-  calculations 

■  experiments  at  7  mm  from  nozzle  exit 

A  experiments  at  70  mm  from  nozzle  exit 

Fig  3  Measured  and  predicted  C2H4  mole  fractions  at  7  mm  from  the  nozzle 
-  calculations 

■  experiments 

Pig  4  Measured  and  predicted  C2H2  mole  fractions  at  7  mm  from  the  nozzle 

-  calculations 

9  experiments 

Pig.  5  Measured  and  predicted  Kb  mole  fractions  at  7  mm  from  the  nozzle 

-  calculations 

■  experiments 

Fig.  6  Measured  and  predicted  OH  mole  fractions 

-  calculations 

■  experiments  at  7  mm  from  nozzle  exit 

A  experiments  at  70  mm  from  nozzle  exit 
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Spatially  Resolved  Measurements  of  Soot  Volume  Fraction 
Using  Laser-Induced  Incandescence 

B.  QUAY.  T.-W.  LEE,  T.  NI,  and  R.  J.  SANTORO* 

Department  of  Mechanical  Engineering,  The  Pennsylvania  State  University,  University  Park,  PA  16802 


Laser-induced  incandescence  is  used  to  obtain  spatially  resolved  measurements  of  soot  volume  fraction  in  a 
laminar  diffusion  flame,  in  which  comparisons  with  laser  scattering/extinction  data  yield  excellent  agreement. 
In  addition,  the  laser-induced  incandescence  signal  is  observed  to  involve  a  rapid  rise  in  intensity  followed  by 
a  relatively  long  (ca.  600  ns)  decay  period  subsequent  to  the  laser  pulse,  while  the  effect  of  laser  fluence  is 
manifest  in  nonlinear  and  near-saturated  response  of  the  laser-induced  incandescence  signal  with  the 
transition  occurring  at  a  laser  fluence  of  approximately  1.2  x  10s  W/cm2.  Spectral  response  of  the 
laser-induced  incandescence  involves  a  continuous  spectrum  in  the  visible  wavelength  range  due  to  the 
blackbody  nature  of  the  emission.  Simultaneous  measurements  of  laser-induced  incandescence  and  light 
scattering  yield  encouraging  results  concerning  the  mean  soot  particle  diameter  and  number  concentration. 
Thus,  laser-induced  incandescence  can  be  used  as  an  instantaneous,  spatially  resolved  diagnostic  of  soot 
volume  fraction  without  the  need  for  the  conventional  line-of-sight  laser  extinction  method,  while  potential 
applications  in  two-dimensional  imaging  and  simultaneous  measurements  of  laser-induced  incandescence  and 
light-scattering  to  generate  a  complete  soot  property  characterization  are  significant. 


INTRODUCTION 

Formation,  growth,  and  oxidation  of  soot  parti¬ 
cles  in  diffusion  flames  involve  a  complex  in¬ 
teraction  between  chemistry  and  fluid  mechan¬ 
ics.  The  understanding  of  these  chemical  and 
physical  processes  is  important  not  only  from  a 
fundamental  scientific  standpoint,  but  also  due 
to  their  applications  in  practical  combustion 
devices.  For  example,  soot  emission  from  auto¬ 
motive  and  gas  turbine  engines  constitutes  one 
of  the  major  pollutants  that  needs  to  be  mini¬ 
mized,  while  excessive  soot  formation  and  radi¬ 
ation  in  propulsion  devices  have  adverse  ef¬ 
fects  on  combustor  and  flow  components.  In 
this  regard,  sooting  characteristics  of  both  tur¬ 
bulent  and  laminar  flames  have  been  investi¬ 
gated  by  numerous  researchers,  while  in  this 
laboratory  attention  has  been  focused  on  ax- 
isymmetric  laminar  diffusion  flames.  The  soot 
property  measurements  made  in  this  flame, 
thus  far,  involve  the  laser  scattering/extinction 
method.  This  technique  yields  measurements 
of  soot  volume  fraction,  mean  soot  particle 
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diameter,  and  number  density  after  tomo¬ 
graphic  inversion  of  the  laser  extinction  data 
due  to  the  line-of-sight  nature  of  these  mea¬ 
surements. 

However,  recent  studies  of  a  process  involv¬ 
ing  laser-induced  incandescence  (LII),  in  which 
the  soot  particles  are  heated  up  by  the  laser 
energy  and  emit  blackbody  radiation  or  incan¬ 
descence  at  elevated  temperatures,  have  shown 
that  LII  can  be  used  as  a  nonintrusive  spatially 
resolved  soot  diagnostic  [1-6].  In  particular,  it 
has  been  pointed  out  by  Melton  [1]  that  the 
LII  signal  is  nearly  proportional  to  the  local 
soot  volume  fraction  for  sufficiently  large  laser 
fluence;  thus,  a  pointwise  measurement  of  soot 
volume  fraction  can  be  made  without  the  need 
for  the  line-of-sight  laser  extinction  and  time- 
consuming  tomographic  inversion  method. 
While  other  applications  of  LII  in.  soot  diag¬ 
nostics,  for  measurements  of  soot  particle  size 
distribution  for  example,  have  been  suggested 
[1],  the  most  direct  and  significant  application 
of  LII  may  be  in  obtaining  point  measure¬ 
ments  of  soot  volume  fraction  since  the  line-in¬ 
tegral  nature  of  laser  extinction  and  subse¬ 
quent  tomographic  inversion  technique  have 
deficiencies  in  some  laminar  flame  and  most 
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turbulent  flame  configurations.  For  example, 
instantaneous  measurements  of  local  soot  vol¬ 
ume  fraction  can  be  made  in  turbulent  diffu¬ 
sion  flames  using  L1I  without  being  limited  to 
time-averaged  data  or  an  axisymmetric  burner 
geometry.  Furthermore,  applications  of  LII  in 
investigations  of  soot  properties  include  two- 
dimensional  imaging  of  soot  volume  fraction 
distributions  and  simultaneous  LII  and  light¬ 
scattering  measurements  to  construct  a  com¬ 
plete  soot  property  characterization.  Recently, 
qualitative  information  on  soot  formation  un¬ 
der  Diesel  engine  conditions  have  been  re¬ 
ported  in  which  simultaneous  two-dimensional 
LII  and  light-scattering  images  were  obtained 
[5,  6], 

In  spite  of  these  potentially  significant  appli¬ 
cations  of  LII  in  soot  diagnostics,  no  experi¬ 
mental  verification  of  the  LII  technique  for 
determining  the  local  soot  volume  fraction  has 
been  reported  to  date.  The  objectives  of  this 
investigation,  therefore,  were  to  experimentally 
determine  the  applicability  of  the  LII  method 
to  obtain  spatially  resolved  measurements  of 
soot  volume  fraction,  to  study  the  feasibility  of 
making  simultaneous  LII  and  light-scattering 
measurements  to  obtain  a  complete  soot  prop¬ 
erty  characterization,  as  well  as  to  investigate 
the  detailed  characteristics  of  LII  in  laminar 
diffusion  flames. 

LASER-INDUCED  INCANDESCENCE 

Laser-induced  incandescence  involves  the 
heating  of  soot  particles  to  temperatures  above 
the  surrounding  gas  temperature  due  to  the 
absorption  of  laser  energy,  and  subsequent  de¬ 
tection  of  the  blackbody  radiation  correspond¬ 
ing  to  the  elevated  soot  particle  temperature. 
The  temperature  of  the  soot  particle  is  deter¬ 
mined  by  the  rate  of  laser  energy  absorption, 
conductive  heat  transfer  to  the  surrounding 
gas,  soot  vaporization,  and  radiative  heat  loss 
through  blackbody  radiation  [I,  2],  For  exam¬ 
ple,  a  Nd:YAG  pulsed  laser  beam  of  ca.  7  ns 
duration  used  in  the  present  LII  measure¬ 
ments  represents  an  energy  source  in  the  en¬ 
ergy  balance  equation,  and  the  soot  particle 
temperature  rapidly  rises  during  the  duration 
of  the  laser  pulse  as  the  soot  particles  absorb 
the  laser  energy.  The  heat  sink  terms  in  this 


phase  are  the  conductive  and  radiative  heat 
losses  to  the  surrounding  gas.  which  are  much 
smaller  than  the  laser  energy  absorption  rate 
for  laser  lluence  levels  relevant  to  LII.  Near  a 
soot  particle  temperature  of  ca.  4000  K,  which 
is  close  to  the  soot  vaporization  point,  the 
temperature  rise  is  severely  curtailed  by  the 
energy  expended  in  the  vaporization  of  soot 
particles  [1],  although  soot  surface  tempera¬ 
tures  as  high  as  5000  K  have  been  observed  for 
sufficiently  large  laser  fluence  [7],  Subsequent 
to  the  laser  pulse,  the  temperature  of  the  soot 
particles  gradually  decreases  due  to  conductive 
and  radiative  heat  losses. 

The  intensity  of  the  LII,  or  the  blackbody 
radiation  due  to  laser  heating,  for  a  single  soot 
particle  has  a  dependence  on  the  soot  particle 
temperature,  detection  wavelength,  and  the 
laser  lluence.  The  total  incandescence  emitted 
from  the  soot  particle  surface  has  a  fourth- 
order  dependence  on  the  soot  particle  temper¬ 
ature.  The  spectral  shape  of  the  incandescence 
is  determined  by  Planck’s  law  with  the  maxi¬ 
mum  in  the  blackbody  radiation  occurring  at  a 
wavelength  inversely  proportional  to  the  soot 
particle  temperature  according  to  Wien’s  dis¬ 
placement  law.  Thus,  the  temporal  variation  in 
the  LII  signal  at  a  given  detection  wavelength 
qualitatively  follows  the  soot  particle  tempera¬ 
ture  in  time,  with  the  exact  functional  relation¬ 
ship  determined  by  the  processes  described 
above. 

Computations  of  the  LII  in  response  to  an 
idealized  laser  pulse  based  on  the  above  black¬ 
body  radiation  laws  and  the  soot  particle  en¬ 
ergy  balance  have  been  performed  by  Melton 
[I]  and  Tait  and  Greenhalgh  [4],  In  particular, 
in  the  limit  of  high  laser  power  and  maximum 
soot  particle  temperature  near  its  vaporization 
point.  Melton  [1]  has  shown  that  the  intensity 
of  the  LII  signal  for  a  group  of  soot  particles 
has  a  dependence  on  mean  soot  particle  diam¬ 
eter  raised  to  the  power  of  (3  +  0.154Adc|), 
where  Adct  is  the  detection  wavelength  ex¬ 
pressed  in  microns.  For  At)c(  between  0.7  and 
0.4  fim,  for  example,  the  LII  signal  is  propor¬ 
tional  to  the  mean  soot  diameter  raised  to  the 
power  of  3.22  to  3.38,  or  approximately  to  the 
soot  volume  fraction.  This  forms  the  basis  for 
the  current  approach  of  using  LII  for  pointwise 
measurement  of  soot  volume  fraction. 
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EXPERIMENTAL  METHODS 

The  experimental  apparatus  involved  a  coan- 
nular  laminar  diffusion  flame  burner  that  was 
identical  to  the  burner  employed  in  this  labo¬ 
ratory  in  previous  studies  of  soot  properties 
[8-10];  thus,  only  a  brief  description  will  be 
given  here.  The  overventilated  laminar  flame 
burner  consisted  of  two  concentric  brass  tubes 
with  fuel  and  air  flowing  through  the  inner 
(11.1  mm  i.d.)  and  outer  (101.6  mm  i.d.)  tubes, 
respectively,  where  the  fuel  tube  extended  4 
mm  beyond  the  exit  plane  of  the  air  tube.  Flow 
conditioning  for  the  air  was  achieved  via  a 
layer  of  3.0-mm  glass  beads,  a  series  of  wire 
screens  and  ceramic  honeycomb  section,  while 
the  fuel  passage  contained  a  layer  of  3.0-mm 
glass  beads  and  a  single  wire  screen.  A  405-mm 
long  brass  cylinder  that  lit  over  the  outer  tube 
was  used  as  a  chimney  to  shield  the  flame  from 
laboratory  air  disturbances;  and  optical  access 
was  obtained  through  machined  slots  on  the 
brass  cylinder  which  traversed  with  the  burner 
assembly.  In  addition,  screens  and  a  flow  re¬ 
strictor  were  placed  at  the  chimney  exit  to 
achieve  a  highly  stable  flame  similar  to  previ¬ 
ous  studies  [8-10].  The  traverse  system  in¬ 
volved  a  stepper  motor  and  controller  (Daedal 
PC-410-288)  that  provided  positioning  capabil¬ 
ity  with  a  resolution  of  0.0127  mm. 

The  optical  setup  for  the  LII  included  an 
Nd:YAG  pulse  laser  (Continuum  Model 
NY61-10),  the  output  beam  of  which  was  ca.  7 
ns  in  duration  and  was  focused  to  approxi¬ 
mately  0.38  mm  in  diameter  using  a  bi-convcx 
lens  of  400  mm  focal  length.  The  probe  volume 
location  was  displaced  ~  25  mm  from  the  focal 
point  of  the  lens  in  order  to  achieve  this  0.38 
mm  diameter.  A  schematic  of  the  optical  ar¬ 
rangement  is  shown  in  Fig.  1.  Both  the  1064-nm 
and  frequency-doubled  532-nm  beams  from  the 
Nd:YAG  laser  were  used,  the  diameter  of 
which  prior  to  the  focusing  lens  was  approxi¬ 
mately  9  mm  with  a  nearly  Gaussian  profile.  In 
order  to  observe  the  effect  of  varying  the  laser 
fluence  for  the  532  nm  wavelength,  a  combina¬ 
tion  of  a  half-wave  plate,  mounted  in  a  rota¬ 
tion  stage,  in  series  with  a  polarizing  beam 
splitter  was  used  to  attenuate  the  laser  energy 
by  a  varying  amount.  A  neutral  density  filter 
(0.7  N.D.)  preceded  the  half-wave  plate  to  re¬ 


duce  the  laser  power  to  an  acceptable  value. 
By  rotating  the  half-wave  plate  to  vary  the 
ratio  of  vertically  polarized  to  horizontally  po¬ 
larized  light,  in  conjunction  with  aligning  the 
beamsplitter  to  transmit  only  the  vertically  po¬ 
larized  portion  of  the  beam,  laser  energies 
ranging  from  0.2  to  3.3  mJ  were  obtained  cor¬ 
responding  to  laser  fluences  between  2.5  X  1()7 
and  4.2  X  10s  W/cm;.  The  laser  energy  was 
monitored  during  the  experiment  using  a  py¬ 
rometer  (Molectron  J 1000),  and  was  main¬ 
tained  during  the  actual  measurements  of  the 
soot  volume  fraction  at  2.4  mJ  for  a  laser 
fluence  of  3.0  X  10*  W/cnr  in  order  to  mini¬ 
mize  the  effect  of  laser  beam  attenuation  across 
the  flame  (see  discussion  below).  The  LII  sig¬ 
nal  was  collected  at  a  90°  angle  using  a  focus¬ 
ing  lens  with  an  /-number  of  3  at  unit  magni¬ 
fication.  Except  for  the  spectral  scans,  an  infer¬ 
ence  filter  (400  +  10  nm)  was  placed  in  front 
of  the  monochromator  to  minimize  interfer¬ 
ence  from  light  scattering.  For  the  spectral 
scans,  the  interference  filter  was  removed  from 
in  front  of  the  monochromator  and  the  1064- 
nm  wavelength  laser  probe  beam  was  used  to 
produce  the  LII  signals.  For  these  measure¬ 
ments,  a  series  of  neutral  density  filters  were 
used  to  vary  the  laser  fluence  rather  than  the 
half-wave  platc/beam  splitter  arrangement 
used  with  the  532-nm  wavelength  laser  probe 
beam. 

Since  the  LII  signal  had  a  continuous  spec¬ 
trum  in  the  visible  wavelength  range,  while 


Fig.  I.  Optical  arrangement  for  laser-induced  incandes¬ 
cence  measurements. 
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interference  from  light-scattering  and  PAH 
fluorescence  was  expected  near  wavelengths  of 
532  nm  and  above  [10,  1 1 1,  the  detection  of  the 
Eli  signal  was  made  at  a  wavelength  of  400  nm 
lor  both  532-  and  1064-nm  wavelength  probe 
laser  beams.  Measurements  made  at  500-  and 
700-nm  detection  wavelengths  for  the  1064  nm 
wavelength  probe  laser  yielded  identical  results 
due  to  the  continuous  nature  of  the  LI  I  spec¬ 
trum  in  the  visible  wavelength  range.  The  de¬ 
tection  wavelength  was  set  by  using  a  0.25-m 
monochromator  (Instruments  SA  H20)  with  a 
grating  blazed  at  330  nm  with  1  mm  x  1  mm 
slits,  which  also  defined  the  length  on  the 
probe  volume.  The  bandpass  of  the  monochro¬ 
mator  was  estimated  to  be  4  nm  FWHM.  while 
the  spectral  response  of  the  monochromator 
was  calibrated  using  an  incandescent  lamp 
(Eppley  T24).  A  photomultiplier  tube 
(Hammamatsu  R928)  was  connected  to  the 
exit  slit  of  the  monochromator,  the  signal  from 
which  was  conditioned  using  a  boxcar  integra¬ 
tor  with  a  gate  width  of  10  ns  and  averaged 
over  100  laser  shots. 

Temporal  variations  of  the  Lll  signal  were 
observable  by  moving  the  boxcar  gate  in  2-1  ti¬ 
ns  increments.  LII  profiles  across  the  flame 
and  spectral  characteristics  of  the  LII  signal 
were  observed  by  using  the  burner  traverse 
system  described  above  and  by  the  scanning  of 
the  detection  wavelength  via  the  monochroma¬ 
tor,  respectively.  The  measurements  were  made 
in  nonsmoking  ethylene/air  diffusion  llames 
where  the  ethylene  and  air  flow  rates  were  3.85 
cmJ/s  and  1060  enr’/s,  respectively. 


RESULTS  AND  DISCUSSION 

Figure  2  shows  a  typical  temporal  variation  of 
the  LII  and  vertically  polarized  light-scattering 
signals  taken  at  a  height  of  40  mm  above  the 
fuel  tube  exit  for  an  ethylene  laminar  diffusion 
flame  at  the  radial  location  where  the  peak 
soot  volume  fraction  is  observed  (r  =  2.5  mm). 
The  variation  of  the  LII  signal  in  time  has 
been  obtained  by  increasing  the  boxcar  gate 
delay  in  2-10-ns  increments  with  respect  to  the 
laser  pulse  while  averaging  over  100  laser  shots, 
as  described  above.  It  can  be  observed  in  Fig.  2 
that  the  initial  phase  of  the  signal  involves  a 


Time  (na) 

f;ig.  2.  Temporal  response  of  the  laser-induced  incandes¬ 
cence. 


rapid  rise  in  the  LII  signal  intensity  due  to  the 
increase  in  the  soot  particle  temperature  dur¬ 
ing  the  laser  pulse  of  ca.  7  ns  in  duration. 
Subsequent  to  the  peak  in  the  LII  signal,  the 
soot  particles  undergo  conductive  and  radiative 
heat  loss  to  the  ambient  gas  and  the  LII  signal 
gradually  decreases,  although  a  sensible  LII 
signal  is  still  observed  at  approximately  600  ns 
after  the  laser  pulse.  The  temporal  variation  of 
the  LII  signal  shown  in  Fig.  2  is  qualitatively 
very  similar  to  the  LII  response  function  for  an 
idealized  laser  pulse  computed  by  Melton  [1], 
A  characteristic  time  constant  for  the  LII  pro¬ 
cess  for  soot  particles  has  been  shown  to  be 
linearly  proportional  to  the  soot  particle  diam¬ 
eter  [3],  and  is  estimated  to  be  approximately 
700  ns  for  a  diameter  of  100  nm.  The  decay 
time  observed  in  Fig.  2  is  approximately  600  ns 
for  the  signal  to  decrease  to  5%  of  the  peak 
value,  while  the  mean  soot  particle  diameter 
( )  at  this  location  is  approximately  135  nm 
(see  Fig.  7b).  Thus,  in  contrast  to  the  light- 
scattering  signal,  which  is  observable  in  princi¬ 
ple  only  during  the  duration  of  the  laser  pulse 
due  to  its  elastic  scattering  nature,  the  LII 
signal  exhibits  a  much  longer  temporal  charac¬ 
teristic.  as  shown  in  Fig.  2.  Melton  [1]  has 
discussed  the  potential  for  obtaining  particle 
size  distribution  parameters  from  the  temporal 
behavior  of  the  LII  signal.  Such  information  is 
better  obtained  using  a  shorter  probe  laser 
wavelength  than  that  used  in  the  present  study 
and  may  also  suffer  from  interference  from 
laser-influenced  fluorescence  from  PAH 
species  [I]. 
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A  comparison  between  the  soot  volume  frac¬ 
tion  measured  by  LI1  and  the  laser 
scattering/extinction  technique  is  shown  in 
Figs.  3a-3c,  where  the  soot  volume  fraction  is 
plotted  as  a  function  of  radial  location  at  se¬ 
lected  heights  ranging  from  10  to  70  mm  above 
the  fuel  tube  exit.  The  open  and  dark  symbols 
represents  laser  scattering/extinction  and  LII 
data,  respectively.  The  laser  scattering/extinc¬ 
tion  data  for  soot  volume  fraction  in  this  flame 
has  been  taken  from  Santoro  et  al.  [8,  9],  and 
involves  the  well-known  method  of  measuring 
the  line-of-sight  extinction  of  the  laser  beam 
followed  by  a  tomographic  inversion  in  order 
to  reconstruct  the  local  soot  volume  fraction. 
Further  details  of  this  technique  and  the  data 
can  be  found  in  Santoro  et  al.  [8,  9],  In  order  to 
calibrate  the  observed  LII  signal,  the  LII  signal 
at  a  single  spatial  point  corresponding  to  the 
radial  location  where  the  peak  soot  volume 
fraction  occurs  (r  =  2.5  mm)  at  the  40  mm 
height  has  been  equated  with  the  known  value 
of  soot  volume  fraction  at  this  location  from 
the  laser  scattering/extinction  measurements. 
All  other  LII  data  can  then  be  converted  to 
absolute  soot  volume  fraction  based  on  this 
single-point  calibration  procedure. 

The  radial  profiles  of  soot  volume  fraction 
obtained  in  this  manner,  as  shown  in  Figs. 
3a-3c,  exhibit  the  familiar  physics  of  soot 
growth  and  oxidization  in  this  flame.  At  low 
heights,  soot  particles  are  observed  in  the  an¬ 
nular  region  on  the  fuel-rich  side  of  the  flame. 
Soot  formed  in  this  region  undergoes  net 
growth  with  increasing  height  up  to  H  =  40 
mm  where  the  peak  soot  volume  fraction  is 
observed  at  a  radial  location  of  2.5  mm  from 
the  centerline,  while  soot  is  observable  at  the 
centerline  at  a  height  of  30  mm  in  Fig.  3a. 
Subsequent  development  involves  a  net  de¬ 
struction  of  soot  particles  through  soot  oxida¬ 
tion  with  soot  volume  fraction  in  the  annular 
region  diminishing  more  rapidly  than  in  the 
central  region. 

Figures  3a-3c  show  an  excellent  agreement 
between  the  LII  and  laser  extinction/scatter¬ 
ing  data  for  the  soot  volume  fraction  with  data 
being  within  5%-10%  of  one  another  at  most 
of  the  heights  where  measurements  have  been 
obtained.  However,  there  is  a  tendency  for  the 
LII  data  to  underestimate  the  soot  volume 


O 


Fig.  3.  Radial  profiles  of  the  soot  volume  fraction  obtained 
via  laser-induced  incandescence  and  laser  scattering/ex- 
tinction  at  several  heights  (//)  above  the  fuel  tube  exit  of 
the  burner  (a)  //  =  10,  20,  and  30  mm;  (b)  //  =  40  and  50 
mm;  and  (c)  H  =  60  and  70  mm. 


fraction  on  one  side  of  the  flame,  resulting  in 
slightly  asymmetric  soot  volume  fraction  pro¬ 
files  in  Fig.  3b.  This  effect  is  more  pronounced 
at  the  height  of  40  mm  than  elsewhere  and  is 
attributable  largely  to  the  fact  that  the  LII 
signal  from  the  far  soot  peak  traverses  the 
flame  in  order  to  arrive  at  the  signal  detection 
site  and,  thus,  is  subject  to  increased  path 
length  and  correspondingly  increased  absorp¬ 
tion  of  the  signal  by  the  soot  and  PAH  species 
in  the  flame.  This  effect  may  be  correctable  by 
estimating  and  integrating  the  local  extinction 
of  the  signal  across  the  flame.  Figure  4  simi¬ 
larly  shows  the  integrated  soot  volume  fraction 
as  a  function  of  height.  The  comparison  be¬ 
tween  the  LII  and  laser  extinction/scattering 
data  [8,  9]  yields  reasonably  good  agreement 
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descence  signal. 


with  the  maximum  discrepancy  being  approxi¬ 
mately  10%  at  a  height  of  30  mm. 

The  effect  of  laser  lluenee  on  the  LI  I  signal 
is  shown  in  Fig.  5.  The  laser  fluence  has  been 
varied  from  2.5  X  107  to  4.2  X  10s  W/cm2 
using  the  combination  of  the  half-wave  plate 
and  beamsplitter,  as  described  earlier.  For  a 
laser  fluence  from  2.5  X  107  to  1.2  x  1 0s 
W/cm2,  it  can  be  seen  in  Fig.  5  that  the  LI  I 
signal  increases  rapidly  as  the  laser  fluence 
increases.  This  effect  is  due  to  the  fact  that  the 
soot  particle  temperature  increases  as  a  func¬ 
tion  of  the  laser  fluence,  which  causes  a  corre¬ 
sponding  increase  in  the  LII  signal.  In  con¬ 
trast,  the  LII  signal  for  laser  fluence  beyond 
ca.  1.2  X  10*  W/cm2  exhibits  a  small  increase 
with  respect  to  an  increase  in  the  laser  fluence. 
The  influx  of  laser  energy  on  the  soot  particles 
can  cause  vaporization  of  small  carbon  frag¬ 
ments,  such  as  C2  and  C\,  from  the  soot  parti¬ 
cle  surface.  For  sufficiently  large  laser  tluence 
this  vaporization  mechanism  and  correspond¬ 
ing  mass  loss  can  become  the  dominant  effect 
which  limits  the  increase  in  soot  particle  tem¬ 
perature  and,  thus,  causes  a  leveling  of  the  LII 
signal  as  shown  in  Fig.  5.  However,  the  LII 
signal  intensity  in  this  “saturation  "  regime  in¬ 
creases  as  a  weak  function  of  the  laser  tluence 
in  Fig.  5  similar  to  the  results  of  Eekbrcth  [7] 
in  which  soot  surface  temperatures  as  high  as 
5000  K  were  observed  with  increasing  laser 
fluence.  From  the  data  shown  in  Fig.  5,  it  is 
estimated  that  the  “saturation  "  regime  for  the 
LII  signal  occurs  for  laser  fluences  above  1.2 
x  10H  W/cm2  in  the  present  studies.  During 
the  actual  measurements  of  soot  volume  frac¬ 


tion.  large  laser  tluence  levels  in  the  “satura¬ 
tion"  regime,  as  were  used  in  the  present  study, 
have  the  advantage  of  being  least  affected  by 
the  ef  fects  of  the  laser  beam  attenuation  across 
the  flame,  since  the  LII  signal  is  a  weak  func¬ 
tion  of  the  laser  tluence  in  this  regime.  For  the 
results  reported  here,  the  LII  response  varied 
by  less  than  5%  for  the  range  of  laser  fluence 
attenuations  ( ~  30%  maximum  attenuation) 
encountered  in  the  present  study.  The  ob¬ 
served  response  of  the  LII  signal  to  the  laser 
fluence  is  sensitive  to  the  laser  beam  intensity 
profile,  the  specific  focusing  arrangement  em¬ 
ployed  for  the  incident  beam  and  the  timing  of 
the  boxcar  gate  pulse  with  respect  to  the  laser 
pulse.  Other  researchers  have  observed  more 
significant  effects  of  the  laser  fluence  on  the 
LII  signal  [12]  and  further  research  is  needed 
in  this  area.  However,  for  the  conditions  de¬ 
scribed  in  the  present  work,  reproducible  re¬ 
sults  were  always  obtained  as  long  as  care  was 
taken  in  the  alignment  of  the  optics  and 
achievement  of  a  near-gaussian  beam  intensity 
profile. 

The  spectral  response  of  the  LII  signal  is 
shown  in  Fig.  6,  where  the  LII  signal  is  plotted 
as  a  function  of  the  detection  wavelength.  The 
spectral  scan  of  the  LII  signal  has  been  ob¬ 
tained  by  rotating  the  grating  in  the  monochro¬ 
mator  so  that  the  detection  wavelength  is  var¬ 
ied  in  l()-nm  increments.  Measurements  have 
been  taken  at  the  radial  location  where  the 
peak  soot  volume  fraction  occurs  at  the  40  mm 
height  for  a  probe  laser  beam  at  the  1064-nm 
wavelength  with  a  beam  diameter  of  0.5  mm 
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Fig.  6.  Spectral  response  of  the  laser-induced  incandes¬ 
cence  where  incident  laser  light  at  a  wavelength  of  1064 
nm  was  used  to  obtain  the  spectrum  shown. 


and  an  energy  of  1.5  mJ /pulse  or  laser  fluence 
of  9.5  x  1()7  W/cnr.  For  the  1064-nm  wave¬ 
length  laser  probe  beam,  interference  from 
light-scattering  and  PAH  fluorescence  is  mini¬ 
mal  in  the  visible  wavelength  range,  although  a 
fortuitous  peak  at  532  nm  is  seen  due  to  the 
leakage  of  the  532-nm  beam  from  the  laser 
and  corresponding  light-scattering  at  this  wave¬ 
length.  It  can  be  observed  that  as  expected  the 
LII  signal  exhibits  a  continuous  spectra  in  the 
visible  wavelength  range  with  the  signals  de¬ 
creasing  to  small  levels  below  300  nm.  The  LII 
spectrum  also  shows  no  distinct  peaks  and  con¬ 
tinues  up  to  750  nm  at  a  nearly  constant  level. 
Comparison  with  computed  blackbody  radia¬ 
tion  curves,  as  shown  in  Fig.  6,  indicates  that 
the  LII  spectrum  in  the  300-450-nm  range 
corresponds  to  a  soot  temperature  between 
5000  and  6000  K,  which  is  higher  than  the 
estimated  soot  vaporization  temperature  of  ca. 
4000  K.  That  the  LII  spectrum  continues  at  a 
nearly  constant  level  up  to  750  nm  and  beyond, 
while  the  computed  radiation  curves  begin  to 
decline,  is  attributed  to  the  fact  that  a  large 
soot  number  density,  and  corresponding  multi¬ 
plicity  of  soot  particle  size,  is  present  in  the 
probe  volume  with  different  soot  surface  tem¬ 
peratures  induced  by  the  laser  fluence.  Hence, 
a  more  continuous  and  level  spectral  response 
is  expected,  as  shown  in  Fig.  6  for  a  group  of 
soot  particles  in  comparison  to  a  radiation 
curve  computed  for  a  single  soot  particle  sur¬ 
face  temperature.  From  Fig.  6,  it  can  also  be 
observed  that  by  using  1064-nm  wavelength 
probe  laser,  LII  measurements  can  be  taken  at 


nearly  anv  wavelength  in  the  visible  range, 
except  near  532  nm  for  our  optical  arrange¬ 
ment,  due  to  the  absence  of  interference  from 
PAH  fluorescence  and  light-scattering  as  noted 
above,  and  measurements  taken  at  detection 
wavelengths  of  400  and  700  nm  have  yielded 
identical  soot  volume  fraction  results. 

The  heating  of  the  soot  particles  to  such 
high  temperatures  may  lead  one  to  be  con¬ 
cerned  about  altering  the  properties  of  the 
soot  particles.  In  the  present  experiments  we 
assume  that  the  soot  properties  are  not  seri¬ 
ously  affected  by  the  large  laser  fluence  before 
the  LII  measurements  are  made.  Furthermore, 
since  the  LII  signals  arc  calibrated  against 
independently  measured  soot  volume  fraction 
values,  any  changes  introduced,  as  long  as  they 
are  the  same  throughout  the  soot  field,  are 
compensated  for  in  the  calibration  procedure. 
These  assumptions  will  need  to  be  verified  in 
future  studies,  but  arc  not  anticipated  to  intro¬ 
duce  errors  any  more  significant  than  the  pre¬ 
sent  knowledge  of  soot  properties  necessarily 
include. 

Figure  7a  shows  the  soot  volume  fraction 
measured  by  LII  and  the  vertically  polarized 
component  of  the  light-scattering  signal 
while  Fig.  7b  contains  the  data  for  the  mean 
soot  particle  diameter  and  the  soot  number 
density  obtained  from  these  measurements.  To 
obtain  the  results  shown  in  Fig.  7b,  the  Qir 
and  soot  volume  fraction  measurements  were 
analyzed  in  a  manner  similar  to  that  described 
by  Santoro  et  al.  [8],  The  vertically  polarized 
light-scattering  signal  has  been  obtained  tor 
the  532-nm  wavelength  laser  probe  beam,  and 
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Fig.  7.  (a)  Radial  profiles  of  the  soot  volume  fraction  and 
vertically  polarized  light-scattering  signal,  Qrr,  at  a  height, 
//,  of  40  mm  above  the  fuel  tube  exit  of  the  burner. 
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Fig.  7.  (b)  Radial  profiles  of  the  mean  soot  particle  diame¬ 
ter.  and  soot  number  concentration,  /V,  at  a  height, 
//,  of  40  mm  above  the  fuel  tube  exit  of  the  burner. 


the  calibration  lor  the  absolute  value  of  Qrr 
has  been  obtained  by  matching  with  the  known 
value  of  this  signal  at  a  single  radial  position 
using  measurements  by  Santoro  et  al.  [8j  simi¬ 
lar  to  the  calibration  procedure  for  the  soot 
volume  fraction.  From  the  soot  volume  frac¬ 
tion  and  light-scattering  data,  which  is  propor¬ 
tional  to  the  sixth  moment  of  the  soot  particle 
diameter  distribution,  the  mean  soot  particle 
diameter  (DM)  and  the  soot  number  density 
( N )  can  be  obtained  as  follows  [8]: 


D„,  =  A4/-1 


'■Qrr 


1/3 


(1) 


where 


F(th) 


nr  -  1 
nr  +  2 


N  = 


12/r 


(2) 

(3) 


Here.  /,,  denotes  the  soot  volume  fraction  while 
the  complex  index  of  refraction.  /?/,  is  taken  as 
(1.57  -  0.56/)  following  Dalzell  and  Sarotim 
[13].  Similar  to  light  scattering/extinction,  the 
LI1  /light  scattering  technique  yields  the  ratio 
of  the  sixth  to  third  moment  (£>,,,)  of  the 
particle  diameter.  The  resultant  mean  soot 
particle  diameter  and  the  soot  number  density 
are  again  compared  with  the  previous  data 
obtained  for  this  flame  using  the  laser  scatter¬ 
ing/extinction  method  [8],  Figure  7b  shows 
that  the  mean  soot  particle  diameters  are  in 


very  good  agreement,  with  the  discrepancy  be¬ 
ing  mostly  limited  to  the  central  region  where 
due  to  the  relatively  low  signal  levels ’there  is 
the  largest  uncertainty  in  both  the  LII  and 
light-scattering  data.  The  mean  soot  particle 
diameter  from  LII  data  in  the  present  study 
ranges  from  75  to  130  nm,  while  diameters  of 
60-135  nm  have  been  observed  by  Santoro  et 
al.  [8], 

It  should  be  noted  that  the  particle  diame¬ 
ters  reported,  in  some  cases,  exceed  the  size 
normally  considered  to  be  appropriate  for  a 
Rayleigh  theory  analysis,  in  fact,  the  soot  par¬ 
ticles  in  this  llame  are  aggregates  composed  of 
primary  particles  that  do  satisfy  the  Rayleigh 
theory  limit  [14].  Thus,  strictly.  D  and  N  refer 
to  the  volume  equivalent  diameter  and  number 
density,  respectively,  of  these  aggregates.  Parti¬ 
cle  morphology  effects  are  not  included  explic¬ 
itly  in  the  current  study,  which  would  require 
multiple-angle  light  scattering  measurements. 
However,  such  measurements  are  feasible,  in 
principle,  and  represent  another  area  where 
the  LII  technique  can  contribute  to  research 
related  to  soot  particle  property  determination. 

The  soot  number  density  profiles  are  simi¬ 
larly  in  reasonably  good  agreement  with  the 
discrepancy  again  occurring  near  the  center- 
line.  Since  the  number  density  is  inversely  pro¬ 
portional  to  the  cube  of  the  mean  soot  particle 
diameter,  the  factor  of  2  difference  in  the  soot 
number  density  near  the  centerline  is  caused 
by  the  corresponding  difference  in  the  mean 
soot  particle  diameter  shown  in  Fig.  7b.  Thus, 
Fig.  7b  shows  that  the  single-point  or  two- 
dimensional  measurements  of  the  LII  and  light 
scattering,  which  can  be  set  up  with  relative 
ease,  a  complete  characterization  of  soot  parti¬ 
cle  properties  may  be  directly  obtained. 


CONCLUSIONS 

From  the  discussion  above,  the  following  con¬ 
clusions  concerning  the  LII  diagnostic  of  soot 
volume  fraction  are  made: 

1.  Laser-induced  incandescence  has  been  used 
to  obtain  spatially  resolved  measurements 
of  soot  volume  fraction  in  laminar  diffusion 
llames,  in  which  comparisons  with  laser 
scattcring/extinction  data  yield  excellent 
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agreement  (or  both  radial  profiles  and  inte¬ 
grated  volume  fraction.  Thus.  LI1  can  be 
used  as  an  instantaneous,  spatially  resolved 
diagnostic  of  soot  volume  fraction  without 
the  need  for  the  conventional  line-of-sight 
laser  extinction  method. 

2.  The  temporal  characteristics  of  the  Lll  sig¬ 
nal  are  observed  to  involve  a  rapid  rise  in 
intensity  followed  by  a  relatively  long  (ca. 
600  ns)  decay  period  subsequent  to  the  laser 
pulse,  while  the  effect  of  laser  fluence  is 
manifest  in  nonlinear  and  saturated  re¬ 
sponse  of  the  LII  signal  with  the  transition 
occurring  at  a  laser  fluencc  of  approxi¬ 
mately  1.2  X  10H  W/crrr  for  a  laser  pulse 
of  ca.  7  ns  in  duration. 

3.  Spectral  response  of  the  LII  involves  a  con¬ 
tinuous  spectrum  in  the  visible  wavelength 
range  due  to  the  blackbody  nature  of  the 
emission,  where  the  spectral  response  for 
300-450-nm  wavelength  range  indicates  a 
soot  surface  temperature  of  ca.  5000  K  with 
the  spectrum  continuing  at  a  nearly  level 
intensity  up  to  750-nm  wavelength  due  to 
the  multiplicity  of  the  soot  particle  sizes  in 
the  probe  volume. 

4.  Simultaneous  measurements  of  LII  and  the 
vertically  polarized  light-scattering  yield  en¬ 
couraging  results  concerning  the  mean  soot 
particle  diameter  and  number  concentra¬ 
tion;  thus,  significant  applications  exist  in 
two-dimensional  imaging  and  simultaneous 
measurements  of  LII  and  light  scattering  to 
generate  a  complete  soot  property  charac¬ 
terization. 
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Two-dimensional  imaging  of  soot  volume  fraction 
by  the  use  of  laser-induced  incandescence 

T.  Ni,  J.  A.  Pinson,  S.  Gupta,  and  R.  J.  Santoro 


recently  developed  laser-induced  incandescence  technique  is  used  to  make  novel  planar  measurements 
of  soot  volume  fraction  within  turbulent  diffusion  flames  and  droplet  flames.  The  two-dimensional 
imaging  technique  is  developed  and  assessed  by  systematic  experiments  in  a  coannular  laminar  diffusion 
name  m  which  the  soot  characteristics  have  been  well  established.  With  a  single  point  calibration 
procedure,  agreement  to  within  10%  was  found  between  the  values  of  soot  volume  fraction  measured  by 
is  technique  and  those  determined  by  conventional  laser  scattering-extinction  methods  in  the 
ame.  As  a  demonstration  of  the  wide  range  of  applicability  of  the  technique,  soot  volume  fraction 
images  are  also  obtained  from  both  turbulent  ethene  diffusion  flames  and  from  a  freely  failing  droplet 
flame  that  burns  the  mixture  of  75%  benzene  and  25%  methanol.  For  the  turbulent  diffusion  flames 
approximately  an  80%  reduction  m  soot  volume  fraction  was  found  when  the  Reynolds  number  of  the  fuel 
je  increased  from  4000  to  8000.  In  the  droplet  flame  case,  the  distribution  of  soot  field  was  found  to  be 
similar  to  that  observed  in  coannular  laminar  diffusion  flames. 

Key  words:  Laser-induced  incandescence,  soot,  flames.  ©  1995  Optical  Society  of  America 


1.  Introduction 

Recently  the  effectiveness  of  laser-induced  incandes¬ 
cence  (LII)  for  the  quantitative  determination  of  local 
soot  volume  fraction  in  laminar  diffusion  flames1-2 
was  demonstrated.  Because  this  diagnostic  tech¬ 
nique  provides  a  direct  measure  of  soot  volume 
fraction  without  the  need  for  a  tomographic  inversion 
process,  as  is  used  in  conventional  line-of-sight  laser 
extinction  methods,  it  was  suggested  that  LII  could 
be  applied  for  instantaneous  planar  measurements  of 
soot  volume  fraction  in  combustion  systems,  includ¬ 
ing  nonsymmetric  or  unsteady  flames.  The  applica¬ 
tions  of  LII  in  qualitative  two-dimensional  (2-D)  soot 
n6ld  imaging  of  diffusion  flame  and  diosol  engine 
combustion  have  been  previously  reported.3-5  With 
the  use  of  a  calibration  scheme  described  in  the 
earlier  work  for  point  measurements,1  it  is  possible  to 
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obtain  quantitative  2-D  imaging  of  soot  volume  frac¬ 
tion. 

LII  arises  from  near-blackbody  emissions  that  origi¬ 
nate  from  laser-heated  soot  particles.  Theoretical 
analysis  has  shown  that  when  the  laser  pulse  is 
sufficiently  intense,  LII  emission  in  the  visible  wave¬ 
length  range  is  approximately  proportional  to  local 
soot  volume  fraction.6  This  prediction  has  been 
verified  experimentally  in  the  previous  work  involv¬ 
ing  point  measurements  in  an  ethene-air  laminar 
diffusion  flame.1  Recent  work  that  address  pre¬ 
mixed  flat  flames  have  observed  similar  results  and 
have  examined  the  effects  of  interfering  radiation 
such  as  background  flame  luminosity  and  fluores¬ 
cence  from  polycyclic  aromatic  hydrocarbons7  (PAH’s). 
Studies  in  which  LII  is  used  to  measure  soot  volume 
fraction  in  flickering  ethene-air  diffusion  flames  have 
also  been  reported.8 

In  this  paper,  we  address  key  issues  associated  with 
quantitative  measurements  made  by  using  2-D  LII 
imaging  and  describe  the  novel  application  of  2-D  LII 
in  the  measurement  of  soot  volume  fraction  in  lami¬ 
nar  diffusion,  turbulent  diffusion,  and  free-falling 
droplet  flames.  To  gain  better  understanding  of  the 
behavior  of  the  LII  technique  and  to  provide  informa¬ 
tion  for  the  design  of  experiments  by  the  use  of  2-D 
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LII  imaging,  we  also  made  a  number  of  point  measure¬ 
ments, 

2.  Experimental  Approach 

A.  Experimental  Apparatus  For  Point  Measurements 

The  optical  setup  for  the  point  measurements  utilizes 
a  frequency-doubled  Nd:YAG  pulsed  laser  (Surelite 
II,  Continuum)  operating  at  10  Hz  with  a  pulse  width 
of  7  ns  FWHM.  For  studying  the  effects  of  the  laser 
beam  intensity  profile  on  the  LII  signal,  two  types  of 
intensity  profiles  were  used.  These  are  referred  to 
as  the  near-Gaussian  profile  and  the  rectangular 
profile  hereafter,  and  they  indicate  the  shape  of  the 
intensity  profile  across  the  laser  beam.  The  rectan¬ 
gular-profile  beam,  i.e.,  one  in  which  the  laser  inten¬ 
sity  is  uniform  across  the  beam,  was  produced  when  a 
1-mm-diameter  pinhole  was  placed  in  front  of  the 
flame,  which  allowed  only  the  center  portion  of  a 
6-mm-diameter  near-Gaussian  laser  beam  to  heat  the 
soot  particles.  Because  the  1-mm-diameter  pinhole 
is  only  6  cm  away  from  the  flame,  the  diffraction 
effects  from  the  pinhole  are  negligible.  The  inten¬ 
sity  of  the  laser  beam  was  varied  by  the  use  of  a 
combination  of  a  half-wave  plate  and  a  polarization 
beam  splitter.  The  LII  signal  at  a  90°  angle  was 
focused  onto  the  entrance  slit  of  a  0.25-m  monochrom¬ 
eter  (Instruments  SA  H20)  and  was  detected  by  a 
photomultiplier  tube  (PMT)  (Hammamatsu  R928) 
wired  for  fast  response.  The  combination  of  a 
0.25-mm  entrance  slit  width  for  the  monochrometer 
and  the  1-mm  beam  diameter  of  the  laser  defined  the 
probe  volume  for  the  point  measurements.  A  nar¬ 
row-band  interference  filter  ( 10  nm  FWHM,  transmit¬ 
tance  peaks  at  400  nm)  was  placed  in  the  front  of  the 
monochrometer  in  order  to  reject  scattered  laser 
light.  The  output  from  the  PMT  was  input  to  a 
100-MHz  digital  oscilloscope  i  Tektronix  TSA  2440) 
and  then  transferred  to  a  lab  computer  by  a  general- 
purpose  interface  bus.  A  coannular  laminar  diffu¬ 
sion  flame  that  burns  ethene  in  air  was  used  in  this 
study.  The  burner  that  was  identical  to  the  one  used 
in  previous  soot  studies9’10  consists  of  a  1.1-cm- 
diameter  fuel  tube  surrounded  by  a  10.2-cm-diameter 
air  tube.  The  flow  rates  of  ethene  and  air  were 
3.85  cm3/s  and  1060  cm3/s,  respectively.  At  this 
flow  condition,  the  distribution  of  soot  volume  frac¬ 
tion  and  soot  particle  size  have  been  measured  by  a 
light  scattering-extinction  method.9 

B.  Experimental  Apparatus  For  Two-Dimensional 
Laser-Induced  Incandescence  Imaging 

The  experimental  setup  for  2-D  LII  imaging  is  shown 
in  Fig.  1.  A  285-mJ/pulse,  532-nm  laser  beam 
(Gaussian  profile  at  far  field)  from  the  YAG  laser  was 
shaped  into  a  350-jxm-thick  (FWHM)  laser  sheet  by  a 
combination  of  a  diverging  cylindrical  lens  [focal 
length  (fl)  =  6  cm]  and  a  spherical  focusing  lens !  fl  =  2 
m).  The  thickness  of  the  laser  sheet  was  measured 
by  translation  of  a  razor  blade  across  the  laser  sheet 
and  measuring  the  transmitted  energy.  For  quanti- 


Fig.  1.  Experimental  setup  for  2-D  LII  imaging. 

tative  2-D  LII  measurements,  the  uniformity  of  the 
laser  sheet  is  particularly  important.  A  good  laser 
sheet  was  obtained  when  the  diverging  cylindrical 
lens  was  placed  close  the  laser  exit,  as  the  laser  beam 
profile  at  near  field  is  nearly  flat.  Although  the  laser 
sheet  at  the  focal  point  was  ~  32  cm  in  height,  only 
the  central  10  cm  were  used  for  2-D  LII  imaging  to 
ensure  the  uniformity  of  the  laser  sheet.  The  unifor¬ 
mity  of  the  center  10  cm  of  the  laser  sheet  was  better 
than  90 7c  as  determined  when  a  2-mm-wide  slit  was 
moved  vertically  along  the  laser  sheet  and  the  trans¬ 
mitted  laser  power  was  measured.  The  use  of  a  very 
long  fl  lens  (fl  =  2  m)  has  the  advantage  of  minimizing 
the  variation  of  the  thickness  of  the  laser  sheet  across 
a  flame,  which  is  especially  important  when  large- 
diameter  flames  are  being  studied.  For  an  imaging 
area  as  wide  as  5  cm,  the  variation  of  the  thickness  is 
estimated  to  be  less  than  57c.  Furthermore,  the 
sensitivity  of  LII  imaging  is  higher  when  a  longer  fl 
lens  is  used,  as  the  LII  signal  is  directly  proportional 
to  the  thickness  of  the  laser  sheet,  which  is  larger  for 
a  longer  fl  lens.  For  large-diameter  flames,  the  use 
of  a  thicker  sheet  is  particularly  important,  as  it  helps 
in  rejecting  background  luminosity,  which  is  nearly 
proportional  to  the  diameter  of  the  flame. 

The  LII  images  were  taken  at  a  90°  angle  to  the 
laser  sheet  by  a  gateable  intensified  CCD  camera 
(Princeton  Instruments,  Model  ICCD-576S/RB) 
equipped  with  a  105-mm  fl  UV  camera  lens  (Nikkon, 
fl 4.5).  The  spatial  resolution  for  these  2-D  imaging 
studies  is  determined  by  the  number  of  camera  pixels 
(578  x  384)  available  and  the  field  of  flew  imaged  in 
the  experiments.  The  resolutions  for  each  of  the 
studies  described  here  are  included  in  the  appropriate 
sections  below.  A  broadband  interference  filter 
(450  ±  25  nm)  was  placed  in  front  of  the  camera  to 
minimize  the  interference  from  scattered  light.  A 
small  amount  of  scattered  light  was  found  to  pass 
through  the  filter.  The  influence  of  the  scattered 
light  was  completely  eliminated  when  the  intensifier 
was  gated  on  after  the  laser  pulse.  Although  a 
stronger  signal  could  be  obtained  if  a  long  gate  time 
were  employed,  the  gate  time  for  the  intensifier  in 
this  experiment  was  set  at  18  ns,  which  is  much 
shorter  than  LII  decay  time,  to  minimize  the  effect  of 
particle  size,  which  is  discussed  below. 

The  interference  from  fluorescence  of  PAH’s  has 
been  examined.  For  ethene  coflow  diffusion  flames, 
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PAH’s  are  mostly  present  at  lower  flame  heights  and 
at  radial  locations  that  are  closer  to  the  center  line 
compared  with  soot  particles.2  In  these  PAH  forma¬ 
tion  regions  where  soot  particles  are  absent,  no 
significant  signal  was  observed  in  the  LII  imaging 
experiments,  indicating  that  fluorescence  from  PAH 
does  not  interfere  with  the  LII  signal.  This  may  be 
due  to  the  observation  that  PAH  species  emit  predomi¬ 
nantly  at  wavelengths  longer  than  the  excitation 
wavelength  (  532  nm),  whereas  the  LII  signal  was 
captured  at  a  shorter  wavelength  (450  nm).  To  pre¬ 
vent  flame  luminosity  from  leaking  through  the 
intensifier  when  the  gate  is  off,  a  mechanical  shutter 
that  opens  for  5  ms  during  image  acquisition  was 
placed  in  front  of  the  camera  lens.  With  a  modest 
intensifier  gain,  the  single-shot  signal  from  10  parts 
in  106  (ppm)  soot  in  an  ethene  diffusion  flame  is 
~  3000  counts  pixel  for  a  probe  volume  of  170  pm  x 
170  pm  x  350  pm.  The  sensitivity  can  be  enhanced 
by  increasing  the  bandwidth  of  the  interference  filter. 
In  the  present  imaging  experiment,  the  camera  was 
gated  on  immediately  after  the  laser  pulse.  When 
the  slightly  delayed  LII  signal  was  measured,  the 
intensity  of  the  acquired  images  was  found  to  remain 
nearly  constant  for  an  average  laser  fluence  beyond 
0.17  J/cm2.  Measurements  were  taken  at  a  laser 
fluence  of  0.32  J  cm2  or  an  intensity  density  of  4.6  x 
107  W/cm2.  This  laser  fluence  is  approximately  a 
factor  of  8  lower  than  that  used  in  our  previous  work 
for  point  measurements.1  The  use  of  low  laser  flu¬ 
ence  has  the  advantage  of  minimizing  intrusion  of  the 
LII  process  on  the  flame  that  is  due  to  laser  vaporiza¬ 
tion  of  soot  particles  and  thus  has  potential  for 
applications  in  systems  in  which  successive  images 
are  needed. 

After  dark  noise  was  subtracted,  the  LII  images 
were  corrected  for  the  combined  spatial  variations  in 
CCD  pixel  response  and  intensifier  gain  when  the 
ratio  of  the  LII  image  to  a  correction  image  was  taken. 
The  correction  image  was  obtained  by  the  uniform 
illumination  of  the  camera  by  a  lamp  located  behind  a 
diffusor.  The  corrected  LII  images  were  then  trans¬ 
lated  into  soot  volume  fraction  images  after  the 
imaging  system  was  calibrated  against  a  sooting 
source  whose  soot  volume  fraction  is  known.  In  this 
experiment,  an  ethene  diffusion  flame  served  as  the 
calibration  standard.1  The  LII  images  were  not 
corrected  for  the  absorption  of  LII  signal  by  soot 
particles  located  between  the  laser  sheet  and  the 
camera.  Although  using  longer  detection  wave¬ 
lengths  could  reduce  the  extinction  of  the  signal,  the 
background  luminosity  becomes  more  difficult  to 
reject. 

3.  Laser-induced  Incandescence  Signal 
Characteristics 

A.  Laser  Fluence  Effects 

The  intensity  of  a  LII  signal  depends  on  the  tempera¬ 
ture  of  the  laser-heated  soot  particles,  which  is  deter¬ 
mined  by  the  rate  of  laser  energy  absorption  and  heat 


loss  through  conduction,  radiation,  and  soot  vaporiza¬ 
tion.  It  has  been  observed1  that  initially  the  LII 
signal  increases  rapidly  as  laser  fluence  increases. 
Once  laser  fluence  reaches  a  saturation  threshold  the 
LII  signal  shows  a  small  or  no  increase  with  respect  to 
an  increase  in  laser  fluence.  At  high  laser  fluence  it 
is  believed  that  soot  vaporization  is  the  dominant 
effect  that  limits  the  increase  in  soot  particle  tempera¬ 
ture.  Using  a  laser  fluence  in  the  saturation  regime 
is  desirable  when  one  is  measuring  soot  volume 
fraction,  as  the  LII  signal  is  then  least  affected  by 
laser  beam  attenuation  across  the  flame.  However, 
there  is  a  significant  uncertainty  in  measuring  the 
local  laser  fluence  in  an  ordinary  laser  beam,  which  is 
not  uniform.  It  is  therefore  difficult  to  determine 
the  saturation  intensity  accurately  under  many  condi¬ 
tions.  In  order  to  study  the  effect  of  laser  fluence  on 
the  LII  signal,  a  rectangular-profile  laser  beam  was 
used  to  heat  the  soot  particles.  Theoretical  mod¬ 
els11’12  for  particle  heating  and  vaporization  have 
been  reported  previously  as  a  function  of  laser  fluence. 
The  use  of  the  rectangular-profile  beam  allows  the 
laser  fluence  to  be  easily  determined,  thereby  allow¬ 
ing  direct  comparison  with  theory. 

Figure  2  illustrates  the  effect  of  laser  fluence  on  the 
intensity  of  the  LII  signal  detected  immediately  after 
the  rectangular-profile  laser  pulse.  These  results 
were  obtained  in  an  ethene  laminar  diffusion  flame  at 
a  height  40  mm  above  the  fuel  tube  exit  and  at  a 
radial  location  where  peak  soot  volume  fraction  is 
observed  for  this  height.  Below  a  laser  fluence 
threshold  of  0.06  J/cm2,  no  observable  LII  signal  was 
detected.  Just  above  the  laser  fluence  threshold,  the 
LII  signal  exhibits  a  rapid  increase  with  increased 
laser  fluence.  After  the  LII  signal  reaches  a  peak  at 
a  laser  fluence  of  0.27  J/cm2,  further  increase  of  laser 
fluence  actually  reduces  the  LII  signal  intensity. 
Similar  results  have  been  observed  recently  by  Vander 
Wal  and  Weiland13,  who  used  a  Gaussian-profile  laser 
beam.  A  profile  that  exhibits  a  rapid  increase  in  the 


LASER  FLUENCE  (J/cm:) 

Fig.  2.  Effect  of  laser  fluence  on  LII  for  a  rectangular  laser 
intensity  profile  (rectangular  profile).  These  measurements  corre¬ 
spond  to  an  ethene-air  laminar  diffusion  flame  at  an  axial  position 
40  mm  above  the  fuel  tube  exit  and  a  radial  location  exhibiting  the 
peak  soot  volume  fraction  for  that  height. 
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(a)  (b) 

Fig.  3.  Effect  of  laser  fluence  on  the  spatially  resolved  LII  signal 
profile  at  a  height  40  mm  above  the  fuel  exit  of  an  ethene-air 
laminar  diffusion  flame:  (a)  laser  fluence  0.15  J/cm2,  (b)  laser 
fluence  0.41  J/cm2. 

LII  signal  followed  by  a  slow  decrease  as  a  function  of 
laser  fluence  has  been  predicted  by  Tait  and  Green- 
halgh3  from  numerical  calculations.  They  predicted 
that  a  decrease  in  particle  mean  volume  that  is  due  to 
vaporization  at  high  laser  fluence  reduces  the  LII 
signal  intensity.  The  spatially  resolved  LII  signal 
profiles  across  an  axisymmetric  ethene  diffusion  flame 
at  the  40-mm  height  are  shown  in  Fig.  3.  Because  of 
the  extinction  of  the  laser  beam  by  soot  particles, 
which  results  in  a  20%-25%  reduction  in  the  laser 
intensity  across  the  flame,  the  LII  profiles  are  asym¬ 
metric.  For  a  laser  fluence  below  saturation,  the  LII 
signal  from  the  peak  that  is  closer  to  the  laser  source 
is  stronger  [shown  in  Fig.  3(a)],  as  the  effect  of  laser 
attenuation  dominates  over  soot  vaporization.  In 
contrast,  for  a  laser  fluence  beyond  the  saturation 
fluence,  the  signal  from  the  far  peak  is  stronger 
[shown  in  Fig.  3(b)],  as  soot  vaporization  is  now  the 
dominant  perturbation  that  affects  the  LII  signal. 

A  temporal  profile  of  the  vertically  polarized  light¬ 
scattering  signal  and  the  influence  of  laser  intensity 
on  the  temporal  profile  of  the  LII  signal  are  shown  in 
Fig.  4.  For  a  laser  intensity  slightly  above  the 


Fig.  4.  Effect  of  laser  fluence  on  the  temporal  profile  of  the  LII 
signal:  dashed  curve,  laser  scattering  signal;  trace  1,  laser  fluence 
0.072  J/cm2;  trace  2,  laser  fluence  0.14  J/cm2;  trace  3,  laser 
fluence  0.27  J/cm2;  trace  4,  laser  fluence  0.72  J/cm2. 


threshold,  the  LII  signal  rises  slowly  and  reaches  its 
maximum  just  after  the  laser  pulse,  which  indicates 
that  the  temperature  of  the  soot  particles  is  con¬ 
stantly  increasing  during  the  laser  pulse.  As  laser 
fluence  increases,  the  rise  time  and  the  decay  time  of 
the  LII  signal  decrease,  while  the  maximum  in  the  LII 
signal  is  observed  to  increase  at  lower  laser  fluence 
(<0.27  J/cm2)  and  decrease  at  higher  laser  fluence 
(>0.27  J/cm2).  The  initial  increase  in  the  maxi¬ 
mum  of  the  LII  signal  and  the  decrease  of  the  LII 
signal  rise  time  with  increased  laser  fluence  are 
consistent  with  the  observation  that  soot  particles  are 
heated  to  a  higher  temperature  and  that  the  tempera¬ 
ture  rises  faster  at  a  higher  laser  intensity.11  The 
dramatic  decrease  of  the  LII  decay  time  with  in¬ 
creased  laser  fluence,  as  shown  in  Fig.  5,  is  not  fully 
understood.  One  explanation  is  that  the  vaporiza¬ 
tion  of  soot  particles  under  intense  laser  heating 
reduces  their  size,  and  the  smaller  particles  conse¬ 
quently  cool  faster.  Dasch12  showed  that  vaporiza¬ 
tion  reduces  the  diameter  of  soot  particles,  D ,  exponen¬ 
tially  with  laser  fluence,  F,  as 

D=  D0  exp[— (jF  -  F*)/F0\9  (1) 

where  D0  is  the  initial  particle  diameter  and  F*  and  FQ 
are  empirical  parameters  that  can  be  determined 
experimentally.  For  a  7-ns  laser  pulse,  the  thresh¬ 
old  fluence  F*  and  decay  fluence  F0  have  been  re¬ 
ported  to  be  0.23  J/cm2  and  1.9  J/cm2,  respectively.12 

For  a  laser  fluence  of  0.54  J/cm2,  the  calculated 
reduction  in  soot  particle  size  from  Eq.  (1)  is  15%. 
Under  the  same  laser  fluence,  the  LII  decay  time  is 
approximately  a  factor  of  10  times  shorter  than  the 
decay  time  under  a  laser  fluence  of  0.11  J/cm2. 
Such  a  significant  decrease  in  the  decay  time  cannot 
be  explained  by  the  vaporization  mechanism  alone. 
It  is  likely  that,  in  addition  to  the  vaporization 
mechanism,  other  physical  or  chemical  properties  of 
the  soot  particles  change  under  the  intense  laser 
fluences  involved  in  the  particle-heating  process. 
Dasch12  measured  light  scattering  and  absorption  of  a 


Fig.  5.  Effect  of  laser  fluence  on  the  LII  signal  decay  time  (the 
time  for  the  signal  to  decay  to  1  /e  of  peak  intensity). 
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cw  probe  laser  along  a  high-energy  pulsed  laser  beam 
and  concluded  that  soot  particles  vaporize  rather 
than  fragment.  However,  by  a  careful  examination 
of  the  time-resolved  light-scattering  signal  obtained 
by  Dasch,12  it  can  be  found  that  there  is  a  transient, 
but  significant,  reduction  in  light  scattering  during  or 
immediately  after  the  laser  pulse,  which  Dasch  attrib¬ 
uted  to  scattered  light  from  the  laser  pulse  or  induced 
incandescence.  This  transient  signal  could  also  be 
interpreted  as  indicating  that  soot  aggregates  un¬ 
dergo  a  rapid  change  in  shape,  resulting  in  an  in¬ 
creased  surface  area  and  a  faster  heat  loss. 

From  these  results,  it  is  clear  that  at  the  laser 
fluences  required  for  LII  measurements,  significant 
effects  on  the  LII  signal  decay  time  can  be  observed  as 
the  laser  fluence  is  varied.  Thus  the  observed  re¬ 
sponse  of  LII  signal  intensity  to  the  laser  fluence  is 
sensitive  to  the  delay  time  for  detector  gate  pulse  with 
respect  to  the  laser  pulse.  The  LII  signal  during  the 
initial  laser-heating  period  is  particularly  sensitive  to 
the  laser  fluence.  F or  quantitative  soot  volume  frac¬ 
tion  measurements,  the  detector  gate  should  be  de¬ 
layed  by  few  nanoseconds  with  respect  to  the  laser 
pulse. 

B.  Particle-Size  Effects 

Figure  6  shows  the  temporal  profiles  of  a  LII  signal 
obtained  in  the  ethene-air  laminar  diffusion  flame  at 
heights  of  10  and  30  mm  above  the  fuel  tube  exit  and 
at  the  radial  locations  corresponding  to  peak  soot 
volume  fraction  for  these  heights.  These  observa¬ 
tions  were  made  with  the  rectangular-profile  laser 
beam  with  a  laser  fluence  of  0.27  J/cm2.  Spatially 
resolved  particle-size  measurements  on  this  flame 
have  been  performed  by  Santoro  et  al9  The  mean 
particle  diameter  (D63)  of  soot  particles  at  the  30-mm 
height  is  160  nm,  which  is  significantly  larger  than 
that  observed  at  the  10-mm  height,  where  the  mean 
size  is  64  nm.  As  can  be  seen  in  Fig.  6,  the  leading 
edges  of  LII  temporal  profiles  at  the  two  locations  in 
the  flame  closely  match,  which  is  consistent  with  the 


Fig.  6.  Temporal  profile  of  a  LII  signal  obtained  in  the  ethene-air 
laminar  diffusion  flame  at  heights  of  10  and  30  mm  above  the  fuel 
tube  exit  and  at  the  radial  locations  corresponding  to  peak  soot 
volume  fraction  for  these  heights. 


theoretical  prediction,3  indicating  that  in  the  Ray¬ 
leigh  regime  all  particles,  regardless  of  size,  increase 
in  temperature  at  the  same  rate.  In  contrast,  the 
LII  signal  at  the  10-mm  flame  height  undergoes  a 
faster  decay  than  that  observed  at  the  30-mm  height. 
This  behavior  is  due  to  the  difference  of  the  soot 
particle  size,  as  smaller  particles  have  a  larger  surface 
area  per  unit  mass  and  are  expected  to  cool  faster. 

To  demonstrate  this  particle-size  effect  and  esti¬ 
mate  the  error  associated  with  the  detection  delay 
time,  a  one-dimensional  (1-D)  LII  imaging  experi¬ 
ment  was  performed.  The  LII  signal  generated  by 
the  rectangular-profile  laser  beam  was  imaged  onto 
the  gated  intensified  CCD  camera.  A  narrow-band 
interference  filter  with  transmittance  centered  at 
400  nm  was  placed  in  the  front  of  the  camera  to  reject 
scattered  laser  light  and  PAH  fluorescence.  The 
camera  was  gated  on  for  5  ns.  Figure  7(a)  shows  the 
normalized  1-D  LII  signal  profiles  from  the  ethene- 
air  laminar  diffusion  flame  at  the  40  mm  height, 
obtained  by  the  varying  of  the  delay  time  of  the 
intensifier  gate  with  respect  to  the  laser  pulse. 
Figure  7(a)  contrasts  the  results  obtained  when  the 
gate  is  placed  immediately  after  the  laser  pulse 
compared  with  the  case  in  which  the  gate  is  delayed 
180  ns.  It  can  be  seen  that  the  LII  signal  ratio  of 
center  line  to  the  peaks  decreases  as  the  intensifier 
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Fig.  7.  Effect  of  the  delay  time  of  the  detector  gate  on  the  spatial 
profile  of  the  LII  signal  at  the  40-mm  height  in  the  ethene-air 
laminar  diffusion  flame:  (a)  solid  curve,  LII  signal  observed 
immediately  after  laser  pulse;  dashed  curve,  LII  signal  observed 
180  ns  after  the  laser  pulse;  (b>  laser  scattering-extinction  data,9 
where  the  solid  curve  is  the  soot  mean  particle  diameter,  and  the 
dashed  curve  is  the  soot  volume  fraction. 
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gate  pulse  was  delayed  for  180  ns  with  respect  to  the 
laser  pulse.  This  observation  is  consistent  with  the 
particle-size  measurements  at  that  height,  which  are 
shown  in  Fig.  7(b).  Compared  with  soot  volume 
fraction  data  at  the  same  flame  height  [shown  in  Fig. 
7(b)]  measured  by  the  laser  scattering-extinction 
method,9  the  LII  profile  taken  immediately  after  the 
laser  pulse  matches  closely  with  the  laser  scattering- 
extinction  data.  Similar  results  were  also  observed 
for  the  Gaussian-profile  laser  beam.  Because  of  this 
size  dependence,  the  delay  time  between  laser  pulse 
and  detector  gate  for  an  imaging  experiment  should 
be  made  as  short  as  possible  to  minimize  this  particle- 
size  effect.  Even  though  the  LII  signal  can  last  as 
long  as  600  ns,  a  short  detector  gate  width  should  be 
used  so  that  the  contribution  from  the  size-sensitive 
portion  of  LII  signal  is  minimized.  For  a  gate  width 
of  50  ns,  the  error  that  is  due  to  the  size  effect  could 
be  as  large  as  10%,  as  estimated  by  the  integration  of 
the  LII  signals,  shown  in  Fig.  6,  over  the  time  period 
of  the  gate  width. 

C.  Laser  Beam  Shape  Effects 

The  results  describe  above  were  obtained  with  the  use 
of  a  laser  pulse  whose  spatial  profile  is  rectangular 
compared  with  a  Gaussian  beam.  However,  typi¬ 
cally  lasers  do  not  generate  a  rectangular-profile 
beam.  Although  the  use  of  a  pinhole  can  reshape  the 
Gaussian  beam  into  a  rectangular-profile  beam,  the 
spatial  resolution  is  relatively  low  because  of  the 
diffraction  effect  from  the  pinhole.  The  use  of  a 
pinhole  also  reduces  the  efficiency  of  laser  power 
utilization.  Thus  it  is  not  feasible  to  use  a  rectangu¬ 
lar-profile  beam  in  a  2-D  LII  imaging  measurement 
that  requires  very  high  laser  energies  to  form  a  laser 
sheet  intense  enough  to  saturate  the  LII  signal. 

For  a  near-Gaussian-profile  beam,  the  laser  inten¬ 
sity  is  not  uniformly  distributed.  Increasing  laser 
energy  will  result  in  an  increase  in  the  effective  LII 
probe  volume,  in  which  the  local  laser  fluence  exceeds 
the  threshold  for  the  LII  signal  to  be  observed. 
Figure  8  shows  a  plot  of  the  LII  signal  taken  immedi¬ 
ately  after  the  laser  pulse  as  a  function  of  average 
laser  fluence  of  the  laser  sheet.  The  LII  signal  was 
obtained  in  the  ethene-air  laminar  diffusion  flame  at 
a  height  of  40  mm  above  the  fuel  tube  exit  and  at  the 
radial  location  that  corresponds  to  peak  soot  volume 
fraction  for  this  height.  The  intensified  CCD  detec¬ 
tor  was  gated  on  for  18  ns.  The  laser  intensity 
across  the  laser  sheet  has  a  near-Gaussian  profile. 
The  thickness  of  the  laser  sheet  is  -350  pm(FWHM). 
The  LII  signal  intensity  shown  in  Fig.  8  is  seen  to  be 
nearly  constant  for  a  laser  fluence  beyond  0.17  J/cm2. 
This  behavior  is  different  from  that  observed  for  the 
rectangular-profile  beam,  in  which  the  LII  signal 
intensity  decreases  after  saturation. 

Previous  point  measurements1  in  which  a  near- 
Gaussian-profile  beam  was  used  have  shown  that  for 
a  sufficiently  intense  laser  fluence  the  LII  signal 
exhibits  a  small  increase  with  an  increase  in  laser 


Fig.  8.  Effect  of  average  laser  fluence  on  LII  for  a  Gaussian- 
shaped  laser  beam.  These  measurements  correspond  to  an  ethene- 
air  laminar  diffusion  flame  at  an  axial  position  40  mm  above  the 
fuel  tube  exit  and  a  radial  location  that  exhibits  the  peak  soot 
volume  fraction  for  that  height. 


fluence.  A  similar  result  was  found  by  Tait  and 
Greenhalgh,3  who  used  a  laser  sheet.  It  was  sug¬ 
gested  that,  as  laser  energy  increases,  the  decrease  in 
the  LII  signal  that  is  due  to  vaporization  of  the  soot 
particles  was  compensated  for  by  an  increase  in  the 
effective  probe  volume,  resulting  in  a  near-constant 
LII  signal  intensity.  We  have  tried  to  generate  LII 
signals  by  using  a  doughnut-shaped  beam.  How¬ 
ever,  the  LII  signal  intensity  observed  continued  to 
increase  with  higher  laser  fluence,  suggesting  that 
the  response  of  the  LII  signal  to  average  laser  fluence 
is  sensitive  to  the  laser  beam  intensity  profile.  This 
result  is  likely  due  to  differences  in  the  expansion  of 
the  probe  volume  with  respect  to  the  degree  of  soot 
particle  vaporization. 

The  interactions  between  laser  beam  and  soot 
particles  may  also  have  significant  effects  on  the  beam 
profile.  For  a  Gaussian  beam,  more  soot  content  is 
vaporized  in  the  center  of  the  beam  than  in  the  outer 
portion,  because  the  laser  intensity  in  the  center  of 
near-Gaussian-profile  beam  is  larger.  The  vaporiza¬ 
tion  effect  decreases  the  absorption  of  laser  energy  by 
soot  particles,  particularly  in  the  center  of  the  laser 
beam.  Consequently,  as  the  laser  beam  propagates 
through  the  soot  field,  the  beam  diameter  becomes 
smaller,  because  the  wings  are  preferentially  absorbed. 
The  vaporization  of  soot  particles  also  affects  the 
temporal  profile  of  the  laser  pulse  as  the  laser  beam 
propagates  through  the  soot  field.  This  effect  is  due 
to  the  fact  that  the  soot  volume  fraction  in  the  laser 
beam  decreases  gradually  during  the  laser  pulse,  the 
less  energy  in  the  later  part  of  the  laser  pulse  is 
absorbed.  The  change  in  the  temporal  profile  of  the 
laser  pulse  may  affect  the  response  of  the  LII  signal  to 
laser  fluence.  To  minimize  the  experimental  error 
caused  by  these  effects,  the  laser  fluence  should  be 
kept  as  low  as  possible.  For  these  reasons,  it  is  not 
recommended  that  the  technique  be  used  in  large- 
diameter  and  heavily  sooting  flames  unless  care  has 
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been  exercised  to  examine  the  interaction  between 
the  laser  beam  and  the  soot  particle  field. 

4.  Imaging  Results  and  Discussion 

In  order  to  examine  and  demonstrate  various  aspects 
of  2-D  LII  imaging  for  soot  volume  fraction  measure- 
ments,  a  series  of  experiments  were  conducted  that 
involved  an  ethene-air  laminar  diffusion  flame,  a 
turbulent  ethene  diffusion  flame  and  a  benzene- 
methanol  droplet  flame.  These  experiments  were 
selected  to  demonstrate 

.  (1)  the  comparative  accuracy  of  the  2-D  LII  tech¬ 
nique  compared  with  previous  laser  scattering- 
extinction  measurements  in  laminar  diffusion  flames’ 

(2)  the  capabilities  of  the  2-D  LII  technique  in 
gaseous,  unsteady  turbulent  diffusion  flames; 

(3)  that  2-D  LII  imaging  can  be  implemented  in 
the  presence  of  a  liquid  fuel  droplet  that  introduces 
additional  light-scattering  interferences. 

A.  Laminar  Diffusion  Flame  Measurements 

The  soot  volume  fraction  image  for  an  ethene— air 
laminar  diffusion  flame  measured  with  the  2-D  LII 
imaging  technique  is  shown  in  Fig.  9.  Because  soot 
particle  size  and  concentration  information  through¬ 
out  the  ethene  diffusion  flame  has  been  deter¬ 
mined,9-10  it  is  possible  to  use  this  flame  for  the 
examination  of  the  accuracy  of  the  2-D  imaging 
technique. 

To  provide  access  for  the  laser  sheet,  the  brass 
chimney  for  the  burner  used  in  previous  work9-10  was 
replaced  with  a  glass  chimney.  Part  of  the  glass 
chimney  was  painted  black  to  minimize  the  reflection 
of  LII  from  the  chimney  wall  onto  the  detector.  The 
flow  conditions  for  the  burner  were  duplicated  from 
previous  work.9  The  laser  sheet  was  directed  to  pass 
through  the  center  line  of  the  flame.  The  flame  has 
a  visible  flame  height  of  88  mm  and  does  not  emit 
smoke  from  the  flame  tip.  Thus  the  10-cm-wide 
laser  sheet  is  large  enough  to  cover  the  entire  soot 
field  in  the  flame.  The  spatial  resolution  for  the  LII 
images  is  170  jrm  x  170  jxm  x  350  jim. 

To  obtain  the  soot  volume  fraction  image,  the 
observed  LII  image  was  scaled  so  that  the  LII  signal 
from  a  pixel  corresponding  to  the  radial  location 
where  the  peak  soot  volume  fraction  occurs  at  the 
40-mm  height  is  equal  to  the  soot  volume  fraction  at 
this  location  from  laser  scattering— extinction  measure- 


Fig.  9.  Soot  volume  fraction  (fv)  image  of  an  ethene-air  laminar 
diffusion  flame.  The  maximum  soot  volume  fraction  in  the  imafre 
is  9.8  x  10"6, 


ments.  A  comparison  between  the  soot  volume  frac¬ 
tion  measured  by  LII  imaging  and  the  laser  scattering- 
extinction  method  for  three  heights  in  the  flame  is 
shown  in  Fig.  10.  It  can  be  seen  that  the  data 
measured  by  the  LII  technique  agrees  very  well  with 
that  obtained  by  the  laser  extinction-scattering 
method.9  The  small  different  (~5%-10%)  between 
LII  data  and  laser  scattering-extinction  data  is  prob¬ 
ably  due  to  the  effect  of  particle  size.  As  can  be  seen, 
the  deviation  usually  occurs  at  lower  heights  or  along 
the  flame  center  line  where  soot  particles  are  smaller. 
Numerical  analysis  performed  by  Melton6  shows  that 
the  intensity  of  LII  signal  at  the  maximum  tempera¬ 
ture  is  proportional  to  the  mean  particle  diameter 
raised  to  the  power  of  (3  +  0.154\det_1),  where  Adet  is 
the  detection  wavelength  in  micrometers.  For  Adet  = 
450  nm,  the  exponent  is  3.34.  The  deviation  of  the 
exponent  from  3  could  result  in  an  underestimation 
of  the  soot  volume  fraction  for  smaller  particles  if  the 
detection  system  is  calibrated  against  larger  particles. 
The  mean  soot  particle  diameter  (D63)  in  the  ethene 
flame  has  been  measured  by  Santoro  et  al9  to  range 
from  15  to  160  nm.  However,  the  difference  be¬ 
tween  the  LII  data  and  the  laser  scattering-extinc¬ 
tion  data  is  significantly  smaller  than  the  theoretical 
prediction. 

B.  Turbulent  Diffusion  Flame  Measurements 

The  turbulent  diffusion  flames  were  established  by 
the  injection  of  ethene  into  still  air  through  a  1-mm- 
diameter  250-mm-long  stainless  steel  tube  mounted 
on  a  traversing  stage.  The  flow  rate  of  ethene  was 
metered  by  a  rotameter.  The  2-D  LII  measurements 
were  made  for  two  different  burner  Reynolds  num¬ 
bers  (4000  and  8000).  The  visible  flame  heights  for 
Re  =  4000  and  Re  —  8000  flames  are  ^  300  mm. 

The  quantitative  single-shot  images  for  soot  vol¬ 
ume  fraction  distribution  were  taken  when  the  10-cm- 
wide  laser  sheet  was  passed  through  the  center  line  of 
the  turbulent  flames,  resulting  in  a  spatial  resolution 
identical  to  the  laminar  flame  case  discussed  above. 


RADIAL  POSITION  (mm) 

Fig.  10.  Comparison  of  the  soot  volume  fraction  profiles  obtained 
by  LII  and  laser  scattering-extinction.9  Solid  curve,  LII  data; 
dashed  curve,  laser  scattering-extinction  data. 
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The  error  involving  the  contribution  from  back¬ 
ground  luminosity  to  the  signal  intensity  in  the 
flames  was  found  to  be  less  than  2%.  However,  it 
must  noted  that  the  extinction  of  the  LII  signal  by 
soot  particles  located  between  the  laser  sheet  and  the 
detector  may  reduce  the  signal  intensity  by  as  much 
as  10%.  Typical  soot  volume  fraction  images,  extend¬ 
ing  over  the  region  from  10  to  20  cm  above  the  fuel 
tube  exit,  are  shown  in  Fig.  11.  The  LII  images, 
which  show  soot  particles  being  confined  in  turbulent 
eddies,  reveal  that  the  soot  distribution  is  highly 
intermittent.  This  soot  field  structure  is  similar  to 
the  observations  made  by  light-scattering  studies  for 
turbulent  acetylene14*15  and  natural  gas16  flames. 
Figure  11  also  shows  that  the  soot  volume  fractions  in 
the  Re  =  8000  flame  are  generally  lower  than  those  in 
the  Re  =  4000  flame  by  at  least  a  factor  of  5.  It  is 
believed  that  the  higher  Reynolds  number  associated 
with  higher  flow  velocity  enhances  the  fuel-air  mix¬ 
ing,  thereby  reducing  the  time  for  soot  growth.  A 
similar  fuel  velocity  effect  on  soot  formation  w^as 
reported  in  the  study  of  turbulent  acetylene  flames  by 
Kent  and  Bastin.17 

C.  Free-Falling  Droplet  Flame  Measurements 

Single  benzene-methanol  droplets  generated  from  an 
aerodynamic  droplet  generator  wrere  injected  down¬ 
ward  into  the  postcombustion  gas  of  a  flat  premixed 
flame  burning  methane.  The  temperature  of  the 
postcombustion  gas,  which  contains  21%  oxygen,  is 
-  670  °C.  To  increase  spatial  resolution  for  LII  im¬ 
aging  in  the  droplet  flame,  the  2-m  spherical  lens  was 
replaced  with  aim  lens  to  shape  the  laser  beam  into 


Fig.  11.  Soot  volume  fraction  i  fL.)  images  of  turbulent  ethene 
diffusion  flames:  a  Reynolds  number  4000,  ib)  Reynolds  number 
8000. 


Fig.  12.  Soot  volume  fraction  l/J  image  of  a  freely  falling  droplet 
flame. 


a  200-  pm-thick  sheet.  A  close-up  bellows  was  placed 
between  the  intensified  CCD  camera  and  the  105-mm 
camera  lens  to  increase  the  magnification  of  the  LII 
images.  For  these  conditions,  the  spatial  resolution 
is  30  gm  x  30  pm  x  200  pm. 

Figure  12  shows  the  distribution  of  soot  formed 
from  a  680-pm-diameter  burning  droplet  that  con¬ 
tains  75  volume  %  benzene  and  25  volume  %  metha¬ 
nol.  The  image  was  taken  50  ms  after  the  droplet 
was  injected  into  the  postflame  gas.  Similar  to  the 
structure  of  coannular  laminar  diffusion  flames,  the 
maximum  soot  volume  fractions  in  this  droplet  flame 
are  found  in  an  annular  region,  not  along  the  center 
line.  Because  benzene  (bp  80  °C)  and  methanol 
(bp  65  °C)  are  highly  volatile,  the  gaseous  fuel  vapor¬ 
ized  from  the  droplet  is  concentrated  in  the  wake 
region  behind  the  droplet,  forming  a  diffusion  flame 
similar  to  the  coannular  gaseous  flame  geometry. 
The  higher  soot  volume  fractions  observed  in  the 
droplet  flame  compared  wdth  the  ethene— air  laminar 
diffusion  flame  is  consistent  with  the  higher  sooting 
tendency  of  benzene.  The  soot  volume  fraction  ob¬ 
served  in  the  droplet  flame  is  sensitive  to  droplet  size 
and  methanol  concentration.  Significant  reductions 
in  soot  volume  fraction  were  observed  by  a  decrease  in 
the  droplet  size  or  an  increase  in  the  concentration  of 
methanol. 

5.  Conclusions 

The  present  investigation  extends  our  previous  work 
on  measurements  of  soot  volume  fraction  by  the  use 
of  LII.  The  results  of  this  investigation  indicate  the 
following: 

(1)  The  application  of  LII  for  instantaneous,  quan¬ 
titative  2-D  measurement  of  soot  volume  fraction  was 
confirmed  by  this  investigation.  There  is  good  agree¬ 
ment  between  the  LII  results  and  laser  scattering- 
extinction  data,  although  the  soot  volume  fractions 
measured  by  LII  are  a  little  lower  in  the  regions 
where  the  mean  soot  particle  diameter  is  small. 

(2)  Time-resolved  measurements  show  that  the 
decay  time  of  the  LII  signal  is  as  long  as  several 
hundred  nanoseconds  and  is  sensitive  to  laser  fluence. 
This  significant  effect  of  laser  fluence  on  the  LII  decay 
time  may  be  related  to  fragmentation  or  shape  changes 
of  the  soot  particles  under  intense  laser  heating. 
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(3)  The  LII  signal  was  found  to  decay  faster  in 
regions  where  the  mean  soot  particle  diameter  is 
smaller.  To  avoid  the  interference  of  elastic  scatter¬ 
ing  and  minimize  particle-size  effects,  detection  should 
be  implemented  immediately  after  the  laser  pulse, 
and  the  sampling  (gate)  time  should  be  short  (less 
than  50  ns). 

(4)  The  technique  has  been  used  to  obtain  2-D 
images  of  the  soot  volume  fraction  field  in  turbulent 
diffusion  flames  and  a  free-falling  droplet  flame, 
demonstrating  that  this  technique  can  provide  quan¬ 
titative,  instantaneous,  and  high  spatially  resolved 
measurements  of  soot  volume  fraction. 
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